Tổng quan nghiên cứu

Biến đổi khí hậu và nước biển dâng đang đặt ra thách thức sống còn đối với nền nông nghiệp Đồng bằng sông Cửu Long, đặc biệt là tại các vùng hạ lưu ven biển. Theo các kịch bản biến đổi khí hậu, nếu mực nước biển dâng cao thêm 1 m, khoảng 43,7% diện tích tự nhiên của tỉnh Sóc Trăng sẽ bị ngập, tác động trực tiếp đến hơn 450.000 người dân, tương đương 35% tổng dân số toàn tỉnh. Huyện Long Phú với diện tích 691 km2 và địa hình thấp trũng ven sông Hậu (cao trình phổ biến từ 1,2 m đến 2,0 m) là một trong những điểm nóng chịu tác động gay gắt nhất của hiện tượng xâm nhập mặn trong mùa khô.

Nghiên cứu tập trung giải quyết bài toán định lượng tác động của độ mặn trong nước tưới và nước ngầm đến sinh trưởng, năng suất của giống lúa OM5451 vụ Đông Xuân tại xã Tân Quy A, huyện Long Phú. Mục tiêu trọng tâm của luận văn là ứng dụng mô hình toán nông học AquaCrop do Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc (FAO) phát triển để mô phỏng chính xác phản ứng của cây lúa dưới các mức độ nhiễm mặn khác nhau; đồng thời dự báo biến động năng suất lúa theo các kịch bản phát thải khí nhà kính B1, B2 và A2 trong các giai đoạn từ năm 2020 đến năm 2050.

Kết quả nghiên cứu cung cấp cơ sở khoa học định lượng với độ tin cậy cao, xác lập các ngưỡng chịu mặn tới hạn và xây dựng gói giải pháp kỹ thuật canh tác tiết kiệm nước thích ứng hiệu quả với biến đổi khí hậu, phục vụ đắc lực cho chiến lược tái cơ cấu ngành nông nghiệp tỉnh Sóc Trăng.

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

Khung lý thuyết áp dụng

Luận văn vận dụng khung lý thuyết cân bằng nước và bốc thoát hơi nước theo tiêu chuẩn FAO kết hợp lý thuyết sinh lý cây trồng chịu mặn. Mô hình AquaCrop phát triển dựa trên nguyên lý tách quá trình bốc thoát hơi nước (ET) thành bốc hơi từ mặt đất (E) và thoát hơi nước qua tán lá cây (Tr), đồng thời chuyển đổi sinh khối tích lũy (B) thành năng suất kinh tế (Y) thông qua chỉ số thu hoạch (Harvest Index - HI) theo phương trình hàm số chuẩn: Y = HI x B.

Các khái niệm chuyên ngành then chốt được chuẩn hóa trong nghiên cứu bao gồm:

  • Độ che phủ tán lá (Canopy Cover - CC): Thông số phản ánh sự phát triển diện tích bề mặt lá và tiến trình khép tán của quần thể cây trồng theo thời gian thực.
  • Hiệu suất sử dụng nước (Water Productivity - WP): Tỷ lệ giữa lượng sinh khối khô tạo thành trên một đơn vị thể tích nước thoát hơi qua cây trồng.
  • Độ dẫn điện của dung dịch đất và nước (Electrical Conductivity - EC): Đơn vị đo độ mặn tính bằng dS/m, phản ánh nồng độ muối hòa tan ảnh hưởng đến thế thẩm thấu của rễ cây.
  • Hệ số thấm bão hòa (Ksat): Chỉ số đánh giá tốc độ chuyển động của nước qua phẫu diện đất trong điều kiện bão hòa theo phân loại O'Neal.

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

Nghiên cứu sử dụng nguồn dữ liệu đa nguồn gồm: chuỗi số liệu khí tượng quan trắc ngày tại trạm Sóc Trăng (nhiệt độ tối cao, nhiệt độ tối thấp, lượng mưa, bức xạ mặt trời, độ ẩm trung bình 83%, lượng bốc thoát hơi tham chiếu ETo tính theo phương pháp FAO Penman-Monteith); số liệu lý hóa phẫu diện đất sét mặn tại xã Tân Quy A; và số liệu thực nghiệm đồng ruộng chi tiết vụ Đông Xuân.

Thí nghiệm đồng ruộng được thiết kế theo khối ngẫu nhiên hoàn toàn (RCBD) trên quy mô 2 công thức so sánh đại diện:

  • Công thức canh tác cải tiến (AWD): Mật độ sạ 120 kg/ha, bón phân cân đối 90 N + 39 P2O5 + 30 K2O kết hợp 500 kg CaO/ha trước khi làm đất, áp dụng quy trình tưới khô ướt xen kẽ với chu kỳ 13 ngày/lần (mực nước 4-6 cm, chỉ tưới bổ sung khi nước ruộng rút xuống mức -15 cm qua ống đo chuyên dụng).
  • Công thức canh tác truyền thống: Mật độ sạ 200 kg/ha, bón phân 120 N + 60 P2O5 + 30 K2O, chu kỳ tưới định kỳ 8 ngày/lần (mực nước duy trì 5-10 cm).

Phương pháp phân tích mô hình hóa bằng AquaCrop được lựa chọn nhờ khả năng mô phỏng chi tiết độ mặn đất và nước ngầm. Quá trình kiểm định và phân tích độ nhạy của mô hình được đánh giá nghiêm ngặt qua các chỉ số thống kê gồm hệ số xác định R2 (yêu cầu đạt trên 0,5), căn bậc hai sai số trung bình bình phương (RMSE), sai số chuẩn hóa (NRMSE) và hệ số hiệu quả mô hình Nash-Sutcliffe (EF).

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

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

Quá trình vận hành và kiểm định mô hình AquaCrop cho thấy độ tương thích vượt trội giữa số liệu mô phỏng và số liệu thực nghiệm với hệ số xác định R2 đạt trên 0,85 và chỉ số sai số chuẩn hóa NRMSE duy trì dưới mức 10%. Các phát hiện định lượng chính bao gồm:

Thứ nhất, độ mặn nước tưới gia tăng gây suy giảm nghiêm trọng đến sinh khối và năng suất hạt. Khi độ dẫn điện của nước tưới tăng từ mức đối chứng lên 4 dS/m, 8 dS/m và 12 dS/m, năng suất lúa OM5451 sụt giảm tương ứng khoảng 12,5%, 24,8% và trên 38,2% so với điều kiện tưới nước ngọt.

Thứ hai, mực nước ngầm và độ mặn nước ngầm có mối tương quan nghịch rõ rệt với năng suất lúa. Ở độ sâu mực nước ngầm nông dưới 0,5 m kết hợp độ mặn nước ngầm từ 8 đến 12 dS/m, lượng muối tích lũy tầng mặt do mao dẫn tăng nhanh, khiến sinh khối cây lúa suy giảm từ 18% đến 32%.

Thứ ba, quy trình canh tác cải tiến chứng minh tính ưu việt vượt bậc khi giúp tiết kiệm từ 22% đến 28% tổng lượng nước tưới trong toàn vụ, giảm 40% lượng giống gieo sạ (từ 200 kg/ha xuống 120 kg/ha) và giảm mức độ tổn thất năng suất từ 7% đến 15% so với tập quán canh tác truyền thống dưới cùng một mức áp lực mặn.

Thứ tư, theo các kịch bản biến đổi khí hậu B1, B2 và A2 trong giai đoạn 2020-2050, nhiệt độ mùa khô tại Sóc Trăng dự báo tăng từ 0,3°C đến 1,2°C trong khi lượng mưa mùa khô giảm từ 1,4% đến 7,7%, làm gia tăng nguy cơ giảm năng suất lúa từ 6,5% (kịch bản B1) lên tới trên 16,8% (kịch bản A2).

Thảo luận kết quả

Cơ chế suy giảm năng suất lúa dưới tác động của mặn xuất phát từ áp suất thẩm thấu dung dịch đất tăng cao, cản trở hệ rễ hấp thu nước và khoáng chất, đồng thời sự tích tụ độc chất của các ion Na+ và Cl- phá vỡ cấu trúc lục lạp, ức chế quá trình quang hợp và phân hóa đòng. Dữ liệu thực nghiệm khi được biểu diễn qua biểu đồ phân tán (scatter plot) và ma trận tương quan đa biến giữa độ mặn nước tưới - độ sâu nước ngầm - năng suất thực tế đã chứng minh rõ xu hướng suy thoái tuyến tính của chỉ số thu hoạch HI khi độ mặn vượt ngưỡng 4 dS/m.

Khi so sánh với các nghiên cứu của Viện Nghiên cứu Lúa Quốc tế (IRRI) và các báo cáo ngành nông nghiệp tại Đồng bằng sông Cửu Long, kết quả mô phỏng từ AquaCrop hoàn toàn tương đồng về quy luật suy giảm năng suất (khoảng 10% cho mỗi độ tăng nhiệt độ kết hợp áp lực hạn mặn). Điều này khẳng định việc áp dụng mô hình toán học là giải pháp thay thế hiệu quả cho các khảo nghiệm đồng ruộng tốn kém, cho phép dự báo chính xác diễn biến mùa vụ trong tương lai dài hạn.

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

Trên cơ sở kết quả mô phỏng và thực nghiệm, luận văn đề xuất 4 nhóm giải pháp hành động cụ thể nhằm bảo vệ sản xuất lúa bền vững tại huyện Long Phú:

  1. Chuyển đổi triệt để sang quy trình kỹ thuật tưới khô ướt xen kẽ (AWD): Triển khai lắp đặt ống nhựa đo mực nước ruộng có đục lỗ (đường kính 100 mm, sâu 20 cm) tại 100% các ô ruộng để nông dân theo dõi ngưỡng rút nước -15 cm trước khi bơm tưới lại; mục tiêu tiết kiệm 25% lượng nước tưới và giảm chi phí nhiên liệu bơm tát; hoàn thành chuyển giao kỹ thuật trong giai đoạn 2024-2026 do Trạm Khuyến nông và Phòng Nông nghiệp huyện Long Phú chủ trì.

  2. Cải tạo nền đất và tối ưu hóa chế độ dinh dưỡng: Áp dụng bắt buộc biện pháp bón lót vôi bột CaO với định lượng 500 kg/ha trước khi làm đất để hạ phèn, đẩy ion Na+ ra khỏi phức hệ hấp thu; giảm lượng phân đạm bón thừa xuống mức chuẩn 90 kg N/ha kết hợp các dạng phân bón chậm tan (Urea 46A+, DAP 46P+); thực hiện ngay từ vụ Đông Xuân tới bởi các Hợp tác xã nông nghiệp.

  3. Hiện đại hóa hệ thống quan trắc và vận hành công trình thủy lợi: Lắp đặt mạng lưới cảm biến đo độ mặn tự động tại các cống đầu mối tiếp giáp sông Hậu, đóng cống triệt để khi độ mặn ngoài sông vượt quá 2 dS/m và vận hành lấy nước ngọt bổ sung vào các kỳ triều thấp; hoàn thành nâng cấp hạ tầng vào năm 2027 do Sở Nông nghiệp và Phát triển nông thôn tỉnh Sóc Trăng điều hành.

  4. Tái cơ cấu lịch thời vụ và chuyển đổi mô hình kinh tế: Dịch chuyển lịch gieo sạ vụ Đông Xuân sớm hơn từ 10 đến 15 ngày (bắt đầu từ tháng 10 dương lịch) để thu hoạch trước đỉnh mặn tháng 3; đồng thời quy hoạch chuyển đổi khoảng 15% đến 20% diện tích đất lúa nhiễm mặn nặng thường xuyên sang mô hình luân canh lúa - tôm sinh thái giai đoạn 2025-2030 dưới sự chỉ đạo của Ủy ban nhân dân huyện Long Phú.

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

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng chính:

  • Cán bộ quản lý nông nghiệp và hoạch định chính sách cấp tỉnh, huyện: Sử dụng kết quả dự báo năng suất theo các kịch bản B1, B2, A2 để xây dựng kế hoạch ứng phó biến đổi khí hậu và quy hoạch vùng sản xuất lúa an toàn.
  • Kỹ sư thủy lợi và cán bộ quản lý tài nguyên nước: Ứng dụng các thông số cân bằng nước, lượng bốc thoát hơi tham chiếu ETo và tính toán lưu lượng tiêu thoát để vận hành tối ưu hệ thống cống đập ngăn mặn, giữ ngọt.
  • Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Khoa học Môi trường, Nông học, Thủy văn: Tiếp cận phương pháp thiết lập tham số, quy trình chuẩn hóa và kiểm định mô hình AquaCrop trong mô phỏng tương tác đất - cây trồng - nước ngầm.
  • Ban chủ nhiệm các Hợp tác xã nông nghiệp và nông dân canh tác lúa vùng ven biển: Áp dụng trực tiếp gói kỹ thuật sạ thưa 120 kg/ha, bón lót vôi tôi khử mặn và kỹ thuật tưới tiết kiệm nước qua ống đo mực nước để giảm 20-30% chi phí đầu vào.

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

Mô hình AquaCrop có điểm gì ưu việt hơn so với các mô hình cây trồng khác như DSSAT khi đánh giá đất nhiễm mặn? AquaCrop tích hợp sẵn mô đun mô phỏng vận chuyển muối và đáp ứng của cây trồng đối với độ mặn của cả nước tưới lẫn nước ngầm, điều mà mô hình DSSAT v4 chưa hỗ trợ đầy đủ. Hơn nữa, mô hình sử dụng chỉ số độ che phủ tán lá CC thay vì chỉ số diện tích lá LAI phức tạp, giúp giảm số lượng tham số đầu vào mà vẫn duy trì độ chính xác cao.

Giống lúa OM5451 chịu ảnh hưởng thế nào khi tưới nước có độ dẫn điện 8 dS/m? Khi độ dẫn điện nước tưới đạt mức 8 dS/m, năng suất thực tế của giống lúa OM5451 bị suy giảm khoảng 24,8% so với đối chứng nước ngọt. Cây lúa xuất hiện các triệu chứng cháy chóp lá, giảm số lượng chồi hữu hiệu và tỷ lệ hạt lép tăng trên 30% do stress thẩm thấu và ngộ độc ion.

Kỹ thuật tưới khô ướt xen kẽ AWD giúp tiết kiệm lượng nước tưới cụ thể là bao nhiêu? Thực nghiệm tại xã Tân Quy A chứng minh kỹ thuật tưới khô ướt xen kẽ kéo dài chu kỳ tưới lên 13 ngày/lần thay vì 8 ngày/lần như tập quán cũ, giúp giảm từ 22% đến 28% tổng lượng nước tưới toàn vụ (tương đương tiết kiệm khoảng 350-450 mm nước) mà không làm sụt giảm năng suất hạt.

Độ sâu mực nước ngầm ảnh hưởng như thế nào đến độ mặn tầng đất mặt ruộng lúa? Khi mực nước ngầm dâng cao cách mặt đất dưới 0,5 m trong điều kiện mùa khô nắng nóng, quá trình bốc hơi bề mặt thúc đẩy mao dẫn đưa các ion muối hòa tan từ nước ngầm lên tích tụ dày đặc ở tầng rễ lúa (0-20 cm), làm tăng độ mặn đất lên gấp 1,5 đến 2 lần so với khu vực có mực nước ngầm sâu hơn 1,0 m.

Theo kịch bản biến đổi khí hậu A2, năng suất lúa tại Long Phú sẽ biến động ra sao trong giai đoạn 2020-2030? Dưới kịch bản phát thải cao A2, nhiệt độ trung bình mùa khô tại Long Phú tăng thêm khoảng 0,4°C đến 0,6°C và lượng mưa mùa khô giảm trên 3,0%, dẫn đến xâm nhập mặn lấn sâu hơn vào nội đồng và làm năng suất lúa OM5451 giảm bình quân từ 8,5% đến 14,2% nếu giữ nguyên tập quán canh tác cũ.

Kết luận

  • Luận văn đã thiết lập thành công quy trình ứng dụng mô hình AquaCrop trong mô phỏng tương tác đất - nước - cây trồng cho giống lúa OM5451 với độ chính xác cao (R2 > 0,85, NRMSE < 10%).
  • Nghiên cứu định lượng rõ mức độ suy giảm năng suất lúa từ 12,5% đến trên 38,2% tương ứng với các dải độ mặn nước tưới từ 4 dS/m đến 12 dS/m tại huyện Long Phú.
  • Khẳng định tính ưu việt của biện pháp canh tác cải tiến (sạ 120 kg/ha, tưới AWD, bón vôi CaO) giúp tiết kiệm hơn 22% lượng nước tưới và giảm thiểu đáng kể tỷ lệ thiệt hại mùa vụ.
  • Cung cấp bức tranh dự báo chi tiết về suy giảm năng suất lúa theo các kịch bản biến đổi khí hậu B1, B2, A2 giai đoạn 2020-2050, làm tiền đề cho quy hoạch nông nghiệp bền vững.
  • Đề xuất đồng bộ gói giải pháp công trình và phi công trình thích ứng với xâm nhập mặn, bảo đảm an ninh lương thực và sinh kế lâu dài cho cộng đồng địa phương.

Đóng góp lớn nhất của công trình là việc chuẩn hóa bộ thông số đầu vào cho mô hình AquaCrop trên đất phèn mặn Đồng bằng sông Cửu Long, mở ra hướng ứng dụng công nghệ mô hình hóa trong quản lý tài nguyên nông nghiệp. Trong giai đoạn tiếp theo từ năm 2024 đến năm 2030, các địa phương cần đẩy mạnh mở rộng mô hình thực nghiệm và tích hợp bản đồ GIS để giám sát mặn theo thời gian thực. Hãy liên hệ với các đơn vị nghiên cứu nông nghiệp và cơ quan khuyến nông địa phương ngay hôm nay để nhận tài liệu chuyển giao quy trình canh tác lúa thông minh thích ứng biến đổi khí hậu.