Tổng quan về luận án

Sự biến đổi sâu sắc của chế độ thủy văn toàn cầu và các hình thái thời tiết cực đoan đang đặt hệ thống nông nghiệp vùng đồng bằng châu thổ trước những thách thức sinh tồn chưa từng có. Tại Việt Nam, Đồng bằng sông Cửu Long (ĐBSCL) giữ vị thế chiến lược khi đóng góp 50% sản lượng lúa, 65% sản lượng nuôi trồng thủy sản, 60% sản lượng cá xuất khẩu, 70% sản lượng trái cây và 95% lượng gạo xuất khẩu của cả nước. Tuy nhiên, toàn vùng có tổng diện tích tự nhiên khoảng 4,0 triệu ha thì có tới xấp xỉ 2,0 triệu ha chịu tác động thường niên của ngập lũ và xâm nhập mặn (XNM). Hiện nay, khoảng 45% diện tích đất tại ĐBSCL đang đối mặt với nguy cơ thoái hóa do nhiễm mặn nghiêm trọng dưới tác động kép của biến đổi khí hậu, nước biển dâng và việc vận hành chuỗi thủy điện bậc thang ở thượng nguồn sông Mê Công.

Khoảng trống nghiên cứu (research gap) cốt lõi tồn tại trong y văn học thuật trước đây là sự thiếu hụt các mô hình lượng hóa động thái chuyển dịch muối trong môi trường đất chưa bão hòa (vadose zone) gắn liền với nguồn nước tưới bề mặt vùng triều, đồng thời chưa tích hợp đồng bộ các mô phỏng vật lý này vào khung đánh giá thích nghi đất đai đa tiêu chí phục vụ quy hoạch không gian. Phần lớn các công trình trước chỉ tiếp cận XNM ở quy mô thủy lực mạng lưới sông kênh vĩ mô hoặc đánh giá tĩnh hiện trạng đất đai, thiếu tính dự báo định lượng theo các kịch bản sử dụng đất nông nghiệp (Land Use Types - LUTs).

Luận án Tiến sĩ Kỹ thuật của Nghiên cứu sinh Nguyễn Ngọc Thy với đề tài "Nghiên cứu các đặc điểm nguồn nước và đề xuất các giải pháp sử dụng đất bị nhiễm mặn do nguồn nước vùng ven biển Đồng bằng sông Cửu Long (Trường hợp nghiên cứu điển hình cho huyện Cần Giuộc, Long An)" (Chuyên ngành: Kỹ thuật Tài nguyên nước, Mã số: 9 58 02 12, thực hiện tại Viện Khoa học Thủy lợi Miền Nam dưới sự hướng dẫn của PGS.TS. Võ Khắc Trí và TS. Hoàng Quang Huy, bảo vệ năm 2024) đã giải quyết trọn vẹn khoảng trống trên.

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học tường minh:

  1. Câu hỏi nghiên cứu 1 (RQ1): Động thái xâm nhập mặn trên hệ thống sông rạch và sự lan truyền mặn vào phẫu diện đất canh tác diễn ra theo quy luật không - thời gian nào trong bối cảnh thủy văn bị điều tiết?
  2. Câu hỏi nghiên cứu 2 (RQ2): Mức độ tích tụ mặn trong tầng đất canh tác ở các mô hình sản xuất nông nghiệp điển hình chịu tác động ra sao khi sử dụng nguồn nước mặt trong mùa khô hạn?
  3. Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để lượng hóa khả năng thích nghi đất đai và tối ưu hóa cơ cấu cây trồng, vật nuôi bền vững dựa trên sự tích hợp giữa mô hình mô phỏng thủy lực - truyền chất và công nghệ GIS?

Các giả thuyết nghiên cứu tương ứng:

  • Giả thuyết 1 (H1): Sự gia tăng độ mặn trong đất canh tác vùng cửa sông ven biển chủ yếu chịu sự chi phối trực tiếp của chế độ sử dụng nước mặt tưới tiêu trong các tháng kiệt đỉnh điểm (tháng 2 - tháng 3), vượt trội so với quá trình mao dẫn mặn từ mực nước ngầm sâu.
  • Giả thuyết 2 (H2): Mô hình lan truyền chất 1D kết hợp hàm truyền thủy lực van Genuchten có khả năng mô phỏng chính xác biến thiên nồng độ muối theo độ sâu phẫu diện đất dưới các kịch bản canh tác khác nhau.
  • Giả thuyết 3 (H3): Tái cấu trúc sử dụng đất dựa trên tích hợp mô hình HYDRUS-1D, khung đánh giá thích nghi FAO và phương pháp phân tích thứ bậc AHP sẽ giảm thiểu tối đa diện tích đất bỏ hoang do nhiễm mặn và tối đa hóa hiệu quả kinh tế - sinh thái.

Phạm vi nghiên cứu tập trung vào vùng cửa sông ven biển chịu ảnh hưởng bán nhật triều Biển Đông tại huyện Cần Giuộc, tỉnh Long An – cửa ngõ chuyển tiếp sinh thái giữa Đông Nam Bộ và Tây Nam Bộ, nơi biên mặn 4 g/l và 1 g/l xâm nhập sâu lần lượt từ 85 - 130 km và 108 - 164 km trên hệ thống sông Vàm Cỏ Đông và Vàm Cỏ Tây.

Literature Review và Positioning

Tổng quan tài liệu cho thấy nghiên cứu về xâm nhập mặn và thoái hóa đất mặn tại ĐBSCL được định hình bởi nhiều dòng nghiên cứu lớn. Về động lực học xâm nhập mặn vùng cửa sông, Lê Sâm (2004) đã phân vùng xâm nhập mặn ĐBSCL thành 4 vùng đặc thù và chỉ ra phạm vi ảnh hưởng mặn vượt trên 50% diện tích toàn vùng. Tăng Đức Thắng và cộng sự (2019) phát triển lý thuyết thành phần nguồn nước trên nền tảng mô hình MIKE 11 để tính toán thủy động lực và lan truyền chất ô nhiễm tại Bán đảo Cà Mau. Viện Khoa học Thủy lợi Miền Nam (2009) ghi nhận độ mặn 10 - 40 g/l xâm nhập sâu từ 50 - 70 km vào nội đồng trong các đợt khô hạn kỷ lục.

Về hóa lý đất mặn và quá trình thoái hóa, Hồ Quang Đức, Nguyễn Văn Đạo, Trương Xuân Cường và Lê Thị Mỹ Hảo (2009) khi đánh giá biến động đất mặn ĐBSCL sau 30 năm đã chỉ rõ diện tích đất tái nhiễm mặn chiếm tới 46% tổng diện tích đất mặn toàn vùng, trong đó nhóm đất mặn ít và trung bình (Molli Salic Fluvisols) chiếm đa số với diện tích khoảng 700.000 ha. Nguyễn Mỹ Hoa, Trần Sơn Tùng, Nguyễn Hồng Giang và Võ Thị Gương (2016) khảo sát tại Bến Tre đã chứng minh hiện tượng mặn hóa đất vào mùa khô đạt ngưỡng sodic hóa (tỷ số hấp phụ natri SAR và tỷ lệ natri trao đổi ESP tăng vọt), tuy nhiên có thể được rửa mặn một phần trong mùa mưa nhờ canh tác lúa. Ngược lại, Nguyễn Đình Vượng, Trần Minh Tuấn, Vũ Thị Thược (2016) tại Bạc Liêu khẳng định mô hình nuôi tôm thâm canh tự phát thiếu quy hoạch thủy lợi là nguyên nhân nhân sinh thúc đẩy nhanh quá trình mặn hóa không thể đảo nghịch, dẫn đến suy thoái đất nghiêm trọng.

Trong y văn tồn tại cuộc tranh luận học thuật sâu sắc giữa hai trường phái:

  • Trường phái kiểm soát công trình cứng (Hard-engineering approach): Ưu tiên xây dựng hệ thống đê bao, cống ngăn mặn khép kín để "ngọt hóa" triệt để phục vụ thâm canh lúa (điển hình tại các dự án Gò Công, Ba Lai, Nam Măng Thít).
  • Trường phái thích ứng sinh thái mềm (Soft-adaptation/Living with salinity): Điển hình là quan điểm của Võ Tòng Xuân (2013), coi nước mặn là tài nguyên kinh tế cho nuôi trồng thủy sản nước lợ và đề xuất nhân rộng mô hình tôm - lúa luân canh kết hợp giám sát mặn dựa vào cộng đồng.

Trên bình diện quốc tế, mô hình hóa sự chuyển dịch muối trong đất được phát triển mạnh mẽ:

  • Madyaka (2008) tại tỉnh Nakhon Ratchasima (Thái Lan) đã tích hợp mô hình Saltmod với kỹ thuật viễn thám và GIS để dự báo không gian đất mặn qua chu kỳ 20 năm canh tác.
  • Daghari (2013) tại vùng Kalâat El Andalous (Tunisia) ứng dụng thành công Saltmod hiệu chỉnh cân bằng nước - muối cho hệ thống tưới tiêu ngầm tầng canh tác các cây trồng cạn và lúa mì.
  • Ting-ting Zhang, Sheng-Lan Zeng, Yu Gao (2011) tại đồng bằng sông Hoàng Hà (Trung Quốc) ứng dụng mô hình tự hồi quy không gian (GIS_SAR) phân tích sự tương quan giữa hoạt động nhân sinh và nguy cơ nhiễm mặn thứ cấp.
  • Mohamed Saidul Islam và Mohammad Atiqur Rahman (2012) cùng Monirul Islam (2014) tại Bangladesh sử dụng GIS và viễn thám lượng hóa tốc độ suy giảm diện tích đất nông nghiệp do XNM, chứng minh xu thế bắt buộc chuyển dịch sang mô hình tôm - cua nước lợ.

Luận án của Nguyễn Ngọc Thy định vị chính xác điểm giao thoa giữa các dòng nghiên cứu trên: không dừng lại ở mô phỏng cân bằng khối vĩ mô (như Saltmod) hay mô hình thủy văn thuần túy, luận án đi sâu vào cơ chế vi mô của dòng thấm chưa bão hòa bằng mô hình toán HYDRUS-1D, kết hợp quan trắc thực nghiệm đa tầng đất và tích hợp đánh giá thích nghi đất đai đa tiêu chí (FAO - MCE - AHP) trên nền tảng GIS cho một tiểu vùng sinh thái cửa sông ven biển đặc thù.

Đó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 lý thuyết chuyển dịch ẩm và chất tan trong môi trường xốp chưa bão hòa được khởi xướng bởi Richards (1931) và hoàn thiện bởi van Genuchten (1980) cùng Mualem (1976), ứng dụng đặc thù vào loại đất phù sa nhiễm mặn ven biển nhiệt đới ẩm gió mùa (Gleyic Salic Fluvisols và Molli Salic Fluvisols). Công trình xây dựng luận chứng lý thuyết về "chu trình tái nhiễm mặn thứ cấp cưỡng bức do nguồn nước mặt": sự suy giảm lưu lượng kiệt thượng nguồn kết hợp với áp lực chuyển đổi cơ cấu sản xuất nông nghiệp (tự phát lấy nước mặn nuôi tôm hoặc tưới lúa muộn) đã làm biến đổi cân bằng áp suất thẩm thấu và phá vỡ cấu trúc trao đổi cation của đất, dẫn đến quá trình sodic hóa tầng rễ cây.

Khung lý thuyết của luận án thiết lập 3 mệnh đề khoa học cơ bản:

  • Mệnh đề 1: Tốc độ tích tụ muối tại tầng canh tác (0 - 30 cm) tỷ lệ thuận với tần suất và nồng độ mặn của nguồn nước mặt được đưa vào đồng ruộng trong các tháng 2 và 3, vượt qua khả năng đệm tự nhiên của dung dịch đất.
  • Mệnh đề 2: Sự suy thoái sức sản xuất của đất mặn là một hàm số phức hợp giữa áp suất thẩm thấu, tỷ số SAR và thời gian duy trì nồng độ Cl- trong mùa kiệt, quyết định ngưỡng chuyển hóa sinh thái của cây trồng.
  • Mệnh đề 3: Đánh giá thích nghi đất đai bền vững đòi hỏi phải đồng quy 3 trụ cột (Thổ nhưỡng - Thủy văn - Hiệu quả kinh tế xã hội) thông qua ma trận trọng số thứ bậc động, thay thế cho các đánh giá tĩnh truyền thống.

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

Khung phân tích của luận án tích hợp liên ngành 3 cấu phần cốt lõi:

[Thủy động lực & Chất lượng nước sông] -> [Mô hình Toán HYDRUS-1D (Phẫu diện đất)]
                                                       |
[Dữ liệu Điều tra RRA & Thổ nhưỡng GIS] -------> [Đánh giá Đa tiêu chí MCE-AHP]
                                                       |
                                            [Quy hoạch Tái cơ cấu LUTs Bền vững]

Phương pháp luận này kết hợp chặt chẽ:

  1. Lý thuyết lan truyền chất theo phương trình bình lưu - phân tán (Advection-Dispersion Equation - ADE).
  2. Quy trình đánh giá đất đai của FAO (1976, 2006) với cấu trúc phân cấp: Loại hình sử dụng đất (LUTs), Hệ thống sử dụng đất (LUS), Yêu cầu sử dụng đất (LUR), Chất lượng đất đai (LQ) và Tính chất đất đai (LC).
  3. Đánh giá đa tiêu chí (Multi-Criteria Evaluation - MCE) phối hợp kỹ thuật liên kết trọng số tuyến tính (Weighted Linear Combination - WLC) và phân tích thứ bậc (Analytic Hierarchy Process - AHP) của Saaty.

Điều kiện biên của khung phân tích được xác định rõ: áp dụng cho vùng cửa sông ven biển chịu chi phối trực tiếp của chế độ bán nhật triều, có mùa khô kéo dài từ tháng 12 đến tháng 5 năm sau, độ mặn nguồn nước mặt dao động từ 1 g/l đến trên 15 g/l.

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

Thiết kế nghiên cứu

Nghiên cứu vận dụng thế giới quan thực chứng kết hợp thuyết hiện thực phê phán (Critical Realism). Thiết kế phương pháp hỗn hợp (Mixed-Methods Research) đa cấp độ được triển khai công phu:

  • Cấp độ vĩ mô (Macro-level): Đánh giá mạng lưới sông rạch chính (sông Vàm Cỏ Đông, sông Vàm Cỏ Tây, sông Soài Rạp, sông Cần Giuộc) qua chuỗi số liệu thủy văn nhiều năm (2002 - 2021).
  • Cấp độ trung mô (Meso-level): Toàn bộ địa bàn huyện Cần Giuộc với phân tích không gian GIS, chồng xếp các lớp dữ liệu thổ nhưỡng, địa hình, hiện trạng đê bao, kênh rạch và phân bố dân cư.
  • Cấp độ vi mô (Micro-level): Thiết lập 14 điểm quan trắc thực nghiệm bằng ống đo chuyên dụng (ống PN), đào và lấy mẫu phân tích 5 phẫu diện đất điển hình đại diện cho các hệ sinh thái sản xuất (ký hiệu CG02 tại xã Phước Vĩnh Tây - chuyên lúa; CG02 tại Phước Vĩnh Đông - chuyên tôm; CG10 tại Tân Tập; CG15 tại Đông Thạnh; CG18 tại Long Phụng).

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm và thu thập dữ liệu tuân thủ nghiêm ngặt các tiêu chuẩn khoa học:

  • Điều tra thực địa: Tiến hành phương pháp đánh giá nhanh nông thôn có sự tham gia (Rapid Rural Appraisal - RRA), phỏng vấn nông hộ theo bảng hỏi cấu trúc tại các xã trọng điểm để thu thập lịch thời vụ, chế độ lấy nước, chi phí đầu vào và năng suất cây trồng, thủy sản.
  • Quan trắc mặn đất và nước: Lắp đặt mạng lưới 14 ống quan trắc PN trực tiếp trên kênh nội đồng và trên mặt ruộng. Thu mẫu đất định kỳ ở các tầng độ sâu 0 - 20 cm, 20 - 40 cm, 40 - 60 cm, 60 - 80 cm và 80 - 100 cm trong các chu kỳ mùa khô (từ tháng 7/2016 đến tháng 7/2018 và mở rộng đến 2021).
  • Phân tích phòng thí nghiệm: Xác định thành phần cơ giới đất, độ dẫn điện của dung dịch đất trích bão hòa ($EC_e$), pH ($pH_{H_2O}$, $pH_{KCl}$), hàm lượng ion hòa tan ($Cl^-$, $SO_4^{2-}$, $Na^+$, $Ca^{2+}$, $Mg^{2+}$, $K^+$), tỷ số hấp phụ natri ($SAR$), tỷ lệ natri trao đổi ($ESP$), hàm lượng chất hữu cơ (mùn), đạm tổng số ($N%$), lân tổng số ($P_2O_5%$) và kali tổng số ($K_2O%$).
  • Kiểm định độ tin cậy và giá trị: Dữ liệu điều tra kinh tế - xã hội được kiểm định phân phối chuẩn thông qua biểu đồ phân phối lũy kế chuẩn P-P Plot và kiểm định dư lượng chuẩn hóa. Độ nhất quán của ma trận trọng số AHP được kiểm soát chặt chẽ với Tỷ số nhất quán $CR < 0,10$ ($CR = CI / RI$).

Data và phân tích

Mô hình toán HYDRUS-1D được ứng dụng để mô phỏng dòng thấm nước và lan truyền muối trong môi trường đất chưa bão hòa 1 chiều thẳng đứng. Phương trình cơ bản mô tả dòng chảy là phương trình Richards phi tuyến:

$$\frac{\partial \theta}{\partial t} = \frac{\partial}{\partial z} \left[ K(h) \left( \frac{\partial h}{\partial z} + 1 \right) \right] - S$$

Trong đó: $\theta$ là độ ẩm thể tích đất ($cm^3/cm^3$); $h$ là thế áp suất nước trong đất ($cm$); $K(h)$ là hàm độ dẫn thủy lực ($cm/ngày$); $z$ là tọa độ không gian thẳng đứng hướng lên ($cm$); $t$ là thời gian ($ngày$); $S$ là số hạng tiêu thoát do rễ cây hút nước ($ngày^{-1}$).

Đặc tính giữ nước của đất được mô tả qua hàm van Genuchten (1980):

$$\theta(h) = \theta_r + \frac{\theta_s - \theta_r}{\left[ 1 + |\alpha h|^n \right]^m} \quad (m = 1 - 1/n)$$

$$K(h) = K_s S_e^l \left[ 1 - (1 - S_e^{1/m})^m \right]^2 \quad \text{với } S_e = \frac{\theta - \theta_r}{\theta_s - \theta_r}$$

Các thông số thủy lực đầu vào ($\theta_r$, $\theta_s$, $\alpha$, $n$, $K_s$, $l$) được xác định thông qua phần mềm RETC dựa trên cơ sở dữ liệu cấp phối hạt USDA và được hiệu chỉnh bằng bài toán nghịch (inverse problem) tối ưu hóa tham số từ dữ liệu quan trắc phẫu diện thực tế (Bảng thông số phẫu diện CG02, CG10, CG15, CG18). Quá trình vận chuyển chất tan ($Cl^-$) được giải theo phương trình Advection-Dispersion phi tuyến có xét đến cơ chế hấp phụ và trao đổi cation.

Đánh giá độ chính xác của mô hình mô phỏng dựa trên chỉ số $R^2$, sai số căn bậc hai trung bình chuẩn hóa ($NRMSE$) và chỉ số hiệu quả Nash-Sutcliffe ($EF$). Kết quả mô phỏng độ ẩm và độ mặn đạt mức tương thích cao giữa số liệu thực đo và mô phỏng qua các điểm nút phẫu diện (Node 1 đến Node 5).

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

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

Luận án đã công bố 5 phát hiện mang tính đột phá về cơ chế nhiễm mặn và sử dụng đất:

  1. Phát hiện 1: Quy luật nghịch đảo thời gian giữa đỉnh mặn và nhu cầu nước sản xuất. Phân tích dữ liệu thủy văn sông Vàm Cỏ Đông và Vàm Cỏ Tây giai đoạn 2002 - 2021 chỉ ra rằng đỉnh xâm nhập mặn hàng năm luôn tập trung cao độ từ cuối tháng 2 đến tháng 3 (trùng khớp với thời kỳ trạm Tân An ghi nhận mặn đạt đỉnh 7,4 - 11,7 g/l, trạm Bến Lức đạt 11,1 - 16,0 g/l). Thời điểm này trùng khớp hoàn toàn với giai đoạn sinh trưởng nhạy cảm nhất và đòi hỏi lượng nước tưới tối đa của cây trồng vụ Đông Xuân, tạo nên "điểm nghẽn sinh thái" khốc liệt khiến năng suất lúa suy giảm nghiêm trọng hoặc mất trắng khi người dân bơm nước mặt trực tiếp vào ruộng.

  2. Phát hiện 2: Diễn biến mặn đất theo độ sâu và vai trò quyết định của nguồn nước tưới bề mặt. Số liệu thực nghiệm từ 14 điểm ống đo PN và 5 phẫu diện chứng minh rằng nguồn gây mặn chính trong đê bao không phải là nước ngầm dâng mao dẫn mà là do quá trình nông dân lấy nước tưới từ hệ thống kênh rạch đã bị nhiễm mặn trong mùa khô. Phân tích các tầng đất cho thấy nồng độ tổng muối hòa tan và $Cl^-$ tích tụ cao nhất ở tầng mặt 0 - 20 cm (đạt mức $>0,65$ g/l, thuộc cấp đất mặn nhiều), sau đó giảm dần ở tầng chuyển tiếp 20 - 40 cm và ổn định ở tầng sâu $>60$ cm.

  3. Phát hiện 3: Nguy cơ mặn hóa và sodic hóa vĩnh viễn từ các mô hình nuôi tôm tự phát. Tại các xã vùng hạ như Phước Vĩnh Đông, Tân Tập, Đông Thạnh, việc người dân tự ý đưa nước mặn ngoài sông vào nội đồng để nuôi tôm thẻ/tôm sú thâm canh sau nhiều vụ thất bại đã làm cho phức hệ trao đổi đất bị bão hòa bởi $Na^+$. Tỷ số SAR tăng cao, kết cấu đất bị phá hủy cấu trúc (chai cứng, phân tán hạt sét), khiến đất bị hoang hóa và mất hoàn toàn khả năng quay trở lại trồng lúa hoặc hoa màu nếu không có giải pháp cải tạo thủy lợi - hóa học đồng bộ trong nhiều năm.

  4. Phát hiện 4: Sự phân hóa không gian thích nghi đất đai của 6 loại hình sử dụng đất (LUTs). Kết quả đánh giá đa tiêu chí MCE-AHP tích hợp GIS trên địa bàn Cần Giuộc cho thấy:

    • Mô hình chuyên lúa 2 vụ ($LUT_1$) bị thu hẹp mạnh diện tích thích nghi cao ($S_1$) do thiếu nước ngọt mùa kiệt;
    • Mô hình lúa - rau màu ($LUT_2$) và chuyên rau ứng dụng công nghệ cao/GAP ($LUT_3$) đạt hiệu quả kinh tế và thích nghi sinh thái cao nhất tại các vùng trung và thượng huyện (nơi có hệ thống đê bao hoàn chỉnh, mặn $<2$ g/l);
    • Mô hình tôm - lúa luân canh ($LUT_4$) và chuyên thủy sản nước lợ ($LUT_5$) thích nghi tối ưu tại các xã ven sông Soài Rạp và sông Cần Giuộc (vùng mặn dao động 4 - 15 g/l).

  5. Phát hiện 5: Độ chính xác cao của mô hình HYDRUS-1D trong dự báo kịch bản mặn tầng canh tác. So sánh kết quả mô phỏng các kịch bản sử dụng đất cho thấy sai số giữa mặn mô phỏng và mặn thực đo tại các phẫu diện thực nghiệm CG02, CG15 năm 2017 duy trì ở mức rất nhỏ ($R^2 > 0,88$, sai số mô hình lúa và tôm nằm trong khoảng tin cậy 95%), xác nhận tính khả thi của việc sử dụng HYDRUS-1D như một công cụ dự báo định lượng phục vụ lập quy hoạch sử dụng đất cấp huyện.

Implications đa chiều

  • Ý nghĩa lý thuyết: Xác lập cơ sở khoa học định lượng về động thái nước - muối trong đất phù sa ven biển ĐBSCL, bắc cầu thành công giữa lý thuyết dòng thấm chưa bão hòa và lý thuyết quy hoạch không gian sử dụng đất bền vững.
  • Ý nghĩa phương pháp luận: Chuẩn hóa quy trình tích hợp liên thông: Khảo sát thực địa (RRA/PN-tube) $\rightarrow$ Thí nghiệm hóa lý $\rightarrow$ Mô hình toán (HYDRUS-1D) $\rightarrow$ Đánh giá đất đai (FAO-GIS-AHP). Quy trình này có tính tổng quát hóa cao, dễ dàng chuyển giao cho các tỉnh ven biển khác trong vùng ĐBSCL như Bến Tre, Trà Vinh, Sóc Trăng, Bạc Liêu, Cà Mau.
  • Ý nghĩa thực tiễn & chính sách: Cung cấp luận cứ kỹ thuật trực tiếp phục vụ triển khai Nghị quyết 120/NQ-CP của Chính phủ về phát triển bền vững ĐBSCL thích ứng với biến đổi khí hậu. Cung cấp bản đồ phân vùng sử dụng đất nông nghiệp chi tiết và giải pháp tái cơ cấu mùa vụ cho Ủy ban nhân dân huyện Cần Giuộc và Sở Nông nghiệp & PTNT tỉnh Long An, ngăn chặn tình trạng chuyển đổi đất tự phát gây tổn thất hàng trăm tỷ đồng mỗi năm.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  • Giới hạn không gian của mô hình thấm 1D: HYDRUS-1D chỉ mô tả dòng thấm và lan truyền chất theo phương thẳng đứng, chưa mô phỏng được dòng chảy ngang (lateral flow) phức tạp qua bờ ruộng, dòng thấm qua khe nứt nẻ của đất sét co ngót trong mùa khô hạn gắt gao.
  • Dữ liệu khí tượng thủy văn ngoại biên: Chuỗi số liệu bốc thoát hơi nước và mưa phụ thuộc chủ yếu vào trạm khí tượng thủy văn Tân An và Bến Lức, chưa bao quát hết sự biến thiên vi khí hậu cục bộ tại các dải đất giáp cửa biển Soài Rạp.
  • Quy mô lấy mẫu thực nghiệm: Mặc dù 14 trạm quan trắc ống PN và 5 phẫu diện đất được lựa chọn đại diện kỹ lưỡng, mật độ này vẫn cần được gia tăng khi áp dụng quy hoạch chi tiết đến từng cụm thửa đất nông hộ nhỏ lẻ.

Định hướng nghiên cứu tiếp theo (Future Research Agenda):

  1. Mở rộng ứng dụng mô hình 2 chiều/3 chiều (HYDRUS-2D/3D hoặc SWAT-MODFLOW) để giải bài toán tương tác thủy lực giữa kênh mương - đê bao - mặt ruộng và dòng chảy ngầm ngang.
  2. Tích hợp các kịch bản biến đổi khí hậu toàn cầu mới nhất (IPCC AR6 - CMIP6) với các kịch bản khai thác thủy điện Mê Công thượng nguồn để dự báo động thái mặn đất đến năm 2050.
  3. Nghiên cứu sâu cơ chế biến tính keo đất và ứng dụng các chế phẩm sinh học, vật liệu cải tạo đất (như biochar, thạch cao phospho-gypsum) kết hợp kỹ thuật tưới tiết kiệm nước nhằm đẩy nhanh tốc độ phục hồi đất hoang hóa sau nuôi tôm.
  4. Xây dựng hệ thống hỗ trợ ra quyết định không gian (Spatial Decision Support System - SDSS) thời gian thực dựa trên nền tảng WebGIS và cảm biến IoT giám sát mặn tự động.

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

Luận án tạo ra những tác động đa tầng đối với khoa học và xã hội:

  • Tác động học thuật: Đặt nền móng cho việc chuẩn hóa phương pháp mô phỏng số động thái mặn tầng đất canh tác trong hệ thống các viện nghiên cứu thủy lợi và trường đại học chuyên ngành Tài nguyên nước, Thổ nhưỡng, Quản lý đất đai.
  • Chuyển đổi ngành nông nghiệp: Thúc đẩy chuyển đổi phương thức sản xuất từ "chống mặn cực đoan" sang "chủ động thích ứng", chuyển đổi hơn 1.000 ha đất canh tác lúa kém hiệu quả vùng hạ Cần Giuộc sang mô hình luân canh tôm - lúa, chuyên tôm công nghệ cao và rau an toàn đạt chuẩn GAP, nâng cao giá trị gia tăng trên một đơn vị diện tích.
  • Tác động chính sách: Cung cấp cơ sở khoa học có độ tin cậy cao cho Sở Tài nguyên & Môi trường, Sở Nông nghiệp & PTNT tỉnh Long An trong việc lập Quy hoạch sử dụng đất thời kỳ 2021 - 2030, tầm nhìn đến năm 2050 và kế hoạch điều tiết vận hành hệ thống công trình thủy lợi ngăn mặn, cấp ngọt vùng 2 sông Vàm Cỏ.
  • Lợi ích kinh tế - xã hội: Giúp cộng đồng nông dân địa phương giảm thiểu rủi ro mất mùa do lấy nhầm nước mặn tưới lúa trong mùa kiệt, ổn định sinh kế nông thôn, bảo vệ tài nguyên đất và nguồn nước ngầm quý giá khỏi nguy cơ ô nhiễm mặn vĩnh viễn.

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

  • Nghiên cứu sinh, Học viên cao học & Nhà nghiên cứu: Tiếp cận quy trình phương pháp luận mẫu mực về mô hình hóa HYDRUS-1D tích hợp viễn thám/GIS và đánh giá thích nghi đất đai FAO.
  • Giảng viên và các cơ sở đào tạo: Bổ sung tư liệu giảng dạy thực tế, cập nhật số liệu chuyên sâu về tài nguyên nước và quản lý đất mặn vùng ĐBSCL.
  • Cơ quan quản lý nhà nước (Bộ Nông nghiệp & PTNT, Bộ TN&MT, UBND các tỉnh, huyện ven biển ĐBSCL): Sử dụng các kết quả đánh giá thích nghi và kịch bản mô phỏng làm công cụ chỉ đạo sản xuất, phân bổ nguồn vốn đầu tư hạ tầng thủy lợi và hoạch định chính sách thích ứng biến đổi khí hậu.
  • Doanh nghiệp nông nghiệp và Nông dân vùng cửa sông ven biển: Ánh xạ lịch thời vụ khoa học, lựa chọn đúng cây trồng, vật nuôi thích nghi với từng ngưỡng mặn, tiếp cận các giải pháp cải tạo đất bị thoái hóa để gia tăng năng suất và thu nhập 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? Trả lời: Đóng góp độc đáo nhất là việc mô hình hóa thành công cơ chế "tích tụ mặn thứ cấp do sử dụng nước tưới bề mặt" trong đất phù sa nhiệt đới bằng phương trình Richards và hàm van Genuchten - Mualem trong mô hình HYDRUS-1D, mở rộng lý thuyết chuyển dịch ẩm và chất tan từ môi trường đồng nhất sang môi trường đất cửa sông phân tầng phức tạp chịu ảnh hưởng của chế độ thủy văn triều.

2. Điểm cải tiến phương pháp luận so với ít nhất 2 nghiên cứu quốc tế tương đồng là gì? Trả lời: So với nghiên cứu của Madyaka (2008) tại Thái Lan và Daghari (2013) tại Tunisia (chỉ sử dụng mô hình cân bằng khối hộp Saltmod với các thông số phân bố trung bình theo tháng), luận án đã cải tiến vượt bậc khi mô phỏng động thái mặn theo từng điểm nút thẳng đứng phẫu diện đất với bước thời gian tính toán linh hoạt theo ngày, đồng thời ghép nối trực tiếp kết quả mô phỏng vào ma trận đánh giá thích nghi đa tiêu chí AHP-GIS, điều mà các công trình trước chưa thực hiện.

3. Phát hiện nào bất ngờ nhất từ dữ liệu thực nghiệm và giải thích cơ chế khoa học? Trả lời: Phát hiện bất ngờ nhất là độ mặn của đất tại tầng rễ cây (0 - 30 cm) trong đê bao tăng vọt trong mùa khô không bắt nguồn từ hiện tượng mao dẫn nước ngầm mặn từ tầng sâu lên (như giả định phổ biến trước đây), mà nguyên nhân chính là do người dân tranh thủ lấy nước mặt ngoài sông kênh vào đồng ruộng trong các đợt triều cường tháng 2 và tháng 3 – đúng thời điểm sông đạt nồng độ mặn đỉnh ($>4 - 11,7$ g/l) nhưng trực quan cảm quan khó nhận biết.

4. Luận án có cung cấp quy trình nhân rộng (Replication Protocol) rõ ràng không? Trả lời: Có. Luận án xây dựng sơ đồ quy trình 5 bước tường minh từ: (1) Thiết lập mạng lưới quan trắc ống PN và lấy mẫu phẫu diện; (2) Hiệu chuẩn thông số thủy lực đất qua phần mềm RETC; (3) Thiết lập điều kiện biên và chạy mô phỏng HYDRUS-1D; (4) Xây dựng bản đồ đơn vị đất đai LUM trên ArcGIS; (5) Tính toán ma trận trọng số AHP theo tiêu chí FAO để phân hạng thích nghi cho từng loại hình sử dụng đất.

5. Chương trình nghiên cứu 10 năm (10-year Research Agenda) được vạch ra như thế nào? Trả lời: Trọng tâm 10 năm tới gồm: Chuyển dịch từ mô hình 1D sang hệ thống mô phỏng ghép nối 3D toàn lưu vực (Surface water - Unsaturated soil - Saturated groundwater); Tích hợp trí tuệ nhân tạo (AI/Machine Learning) dự báo sớm độ mặn đất dựa trên dữ liệu viễn thám radar vệ tinh (Sentinel-1/2); Nghiên cứu công nghệ sinh học nano cải tạo đất sodic hóa sau nuôi tôm; Thiết lập mạng lưới quan trắc thông minh IoT kết nối trực tiếp với ứng dụng chỉ đạo mùa vụ trên thiết bị di động của nông dân.

Kết luận

  1. Luận án đã làm sáng tỏ một cách toàn diện, khoa học và có hệ thống đặc điểm nguồn nước mặt, quy luật xâm nhập mặn trên hệ thống sông Vàm Cỏ và mối tương quan nhân quả trực tiếp đến quá trình mặn hóa đất canh tác nông nghiệp vùng ven biển huyện Cần Giuộc, tỉnh Long An.
  2. Ứng dụng thành công và chuẩn hóa mô hình toán HYDRUS-1D trong việc mô phỏng chính xác diễn biến độ mặn theo thời gian và không gian phẫu diện đất cho các mô hình sản xuất chuyên lúa, luân canh rau màu và nuôi trồng thủy sản.
  3. Đột phá trong phương pháp luận khi tích hợp thành công dữ liệu mô phỏng thủy lực ngầm với hệ thống thông tin địa lý GIS và đánh giá đa tiêu chí MCE-AHP theo tiêu chuẩn của FAO, thiết lập cơ sở khoa học vững chắc cho công tác phân vùng thích nghi đất đai.
  4. Đề xuất mô hình tái cơ cấu sử dụng đất nông nghiệp bền vững với cơ cấu mùa vụ, cây trồng, vật nuôi tối ưu cho từng tiểu vùng sinh thái (thượng - trung - hạ) của huyện Cần Giuộc, hóa giải xung đột mặn - ngọt trong sản xuất.
  5. Tạo bước chuyển dịch quan điểm (paradigm shift) từ tư duy công trình "ngăn mặn, ngọt hóa triệt để" sang tư duy "chủ động thích ứng sinh thái, quản trị rủi ro mặn", góp phần hiện thực hóa mục tiêu phát triển bền vững ĐBSCL theo tinh thần Nghị quyết 120/NQ-CP của Chính phủ.
  6. Mở ra 3 hướng nghiên cứu chuyên sâu mới cho ngành kỹ thuật tài nguyên nước và khoa học đất: mô hình hóa 3 chiều tương tác nước mặt - đất - nước ngầm vùng cửa sông; dự báo mặn hóa dưới kịch bản đa hồ chứa thượng nguồn Mê Công; và công nghệ tích hợp IoT - GIS phục vụ quản lý thích ứng đất nhiễm mặn trong thời gian thực.