Tổng quan luận án
Luận án tiến sĩ với đề tài "Đánh giá nguy cơ môi trường do sử dụng phân đạm trên đất trồng cây hàng năm tại huyện Giao Thủy, tỉnh Nam Định" do nghiên cứu sinh Lý Thị Thu Hà thực hiện tại Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội (năm 2023), thuộc chuyên ngành Môi trường đất và nước (mã số: 9440301.02), dưới sự hướng dẫn khoa học của PGS. Nguyễn Thị Hà và PGS. Ngô Thế Ân. Công trình được hỗ trợ cơ sở vật chất từ Bộ môn Công nghệ Môi trường (Trường Đại học Khoa học Tự nhiên và Học viện Nông nghiệp Việt Nam) cùng nguồn tài trợ từ Dự án DeltAdapt.
Tính cấp thiết và khoảng trống nghiên cứu
Nitơ (N) là nguyên tố dinh dưỡng thiết yếu quyết định sinh trưởng và năng suất cây trồng. Trong các hệ thống nông nghiệp thâm canh, tỷ lệ bón phân đạm trung bình toàn cầu đạt khoảng 220 kg N/ha/vụ trong sản xuất rau. Tuy nhiên, cây trồng thường chỉ hấp thu khoảng 47–50% lượng phân đạm được cung cấp; lượng dư thừa còn lại phân tán vào các thành phần môi trường. Lượng đạm dư thừa thất thoát thông qua các con đường: rửa trôi nitrat ($\text{NO}_3^-$) xuống tầng nước dưới đất, bay hơi amoniac ($\text{NH}_3$), và chuyển hóa thành các dạng khí như dinitơ ($\text{N}_2$), oxit nitric ($\text{NO}$), và nitơ oxit ($\text{N}_2\text{O}$). Trong đó, $\text{NO}_3^-$ là tác nhân gây ô nhiễm nước ngầm và phú dưỡng nguồn nước mặt, còn $\text{N}_2\text{O}$ là khí nhà kính có tiềm năng làm nóng toàn cầu trong 100 năm gấp 265 đến 298 lần so với $\text{CO}_2$, đồng thời phá hủy tầng ôzôn bình lưu. Đất trồng cây hàng năm đóng góp xấp xỉ 60% tổng lượng phát thải $\text{N}_2\text{O}$ toàn cầu từ hoạt động nhân sinh.
Tại huyện Giao Thủy, tỉnh Nam Định – vùng đất ven biển cửa sông Hồng chịu tác động mạnh của biến đổi khí hậu và xâm nhập mặn – quá trình chuyển đổi cơ cấu cây trồng từ đất chuyên lúa sang mô hình luân canh rau – lúa và chuyên rau diễn ra mạnh mẽ trên nền đất pha cát cồn cao. Quá trình thâm canh này kéo theo sự gia tăng lượng phân đạm bón vào đất. Tuy nhiên, các đánh giá định lượng đồng bộ về rủi ro rửa trôi $\text{NO}_3^-$ tới nước dưới đất tầng nông và mức phát thải $\text{N}_2\text{O}$ tại các vùng ven biển nhiệt đới ở Việt Nam còn rất hạn chế. Việc kết hợp các mô hình không gian và phương pháp đo đạc thực địa để xác định rủi ro môi trường do phân đạm là yêu cầu cấp thiết.
Mục tiêu nghiên cứu
- Mục tiêu chung: Đánh giá nguy cơ, rủi ro môi trường từ nitơ (dạng $\text{NO}_3^-$) rửa trôi tới nước dưới đất tầng nông và phát thải khí nhà kính ($\text{N}_2\text{O}$) liên quan đến sử dụng phân đạm trên đất trồng cây hàng năm tại huyện Giao Thủy, tỉnh Nam Định.
- Mục tiêu cụ thể:
- Đánh giá hiện trạng sử dụng phân đạm trên đất trồng cây hàng năm tại huyện Giao Thủy, tỉnh Nam Định.
- Đánh giá nguy cơ, rủi ro rửa trôi $\text{NO}_3^-$ tới nước dưới đất tầng nông và phát thải khí nhà kính $\text{N}_2\text{O}$ liên quan đến sử dụng phân đạm trên đất trồng cây hàng năm tại huyện Giao Thủy, tỉnh Nam Định.
- Đề xuất các giải pháp quản lý rủi ro môi trường thông qua đánh giá cân bằng nitơ trên đất trồng cây hàng năm nhằm giảm thiểu rửa trôi $\text{NO}_3^-$ tới nước dưới đất tầng nông và phát thải khí nhà kính $\text{N}_2\text{O}$ tại địa bàn nghiên cứu.
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu: Các loại hình sử dụng đất trồng cây hàng năm (ĐTCHN); lượng phân đạm vô cơ và hữu cơ sử dụng; nước dưới đất tầng nông (tầng chứa nước lỗ hổng Holocen - $\text{qh}_2$); nồng độ $\text{NO}_3^-$ rửa trôi; dòng phát thải khí nhà kính $\text{N}_2\text{O}$.
- Phạm vi không gian: Địa bàn huyện Giao Thủy, tỉnh Nam Định với 22 xã và thị trấn đại diện cho các vùng sinh thái nông nghiệp ven biển.
- Phạm vi thời gian: Số liệu điều tra nông học, quan trắc môi trường đất, nước và đo đạc khí nhà kính thực địa được thực hiện trong giai đoạn 2015–2017.
Tổng quan tài liệu và vị trí của luận án
Các hướng nghiên cứu chính đã được tổng quan
- Nhu cầu và hiệu quả sử dụng phân đạm (NUE): Báo cáo của FAO (2017, 2018) cho thấy nhu cầu đạm toàn cầu và châu Á gia tăng liên tục (châu Á chiếm hơn 60% tổng lượng tiêu thụ toàn cầu). Lassaletta và cộng sự (2014, 2016) cùng Yan và cộng sự (2020) chỉ ra xu hướng suy giảm NUE toàn cầu từ 68% (năm 1960) xuống còn khoảng 47%. Trong canh tác lúa gạo, các công trình của Peng và cộng sự (2009), Chen và cộng sự (2014), Liang Wu và cộng sự (2015), Ding và cộng sự (2020), Chivenge và cộng sự (2021), Alam (2023) thống kê lượng bón đạm dao động mạnh từ 33–350 kg N/ha/vụ tùy quốc gia. Trong canh tác rau, mức bón đạm trên thế giới rất cao: New Zealand từ 100–230 kg N/ha (Reid & Morton, 2019); Georgia, Hoa Kỳ từ 191–315 kg N/ha (da Silva, 2020); Hàn Quốc từ 487–500,5 kg N/ha (Koh, 2021); Trung Quốc từ 234–870 kg N/ha (Min & Shi, 2018; Jia et al., 2012). Tại Việt Nam, tiêu thụ phân bón đạt 415,3 kg/ha năm 2018 (FAO, 2018); các nghiên cứu của Nguyen (2017), Phan (2018), Trần Thị Xuân Phương (2021) ghi nhận lượng đạm bón thực tế cho lúa tại Đồng bằng sông Cửu Long và Thừa Thiên Huế thường vượt khuyến cáo từ 20–30%.
- Hệ quả môi trường do phân đạm:
- Tích lũy nitrat trong rau: Colla và cộng sự (2018), Kyriacou và cộng sự (2019) cảnh báo nguy cơ hình thành hợp chất N-nitroso gây ung thư và hội chứng methemoglobin khi bón vượt 200 kg N/ha.
- Axit hóa đất: Guo và cộng sự (2010), Chien và cộng sự (2008), Lucas và cộng sự (2011) chứng minh phân bón gốc amoni làm giảm pH đất từ 0,30–0,80 đơn vị sau nhiều năm thâm canh.
- Bay hơi amoniac và rửa trôi nitrat: Zhou và cộng sự (2021) ước tính $\text{NH}_3$ bay hơi chiếm 10–14% lượng đạm thất thoát; Jabloun và cộng sự (2015), Tamagno và cộng sự (2022), Hess và cộng sự (2020) xác định rửa trôi $\text{NO}_3^-$ chiếm khoảng 19% tổng lượng N bón toàn cầu, chịu tác động mạnh bởi lượng mưa và kết cấu đất cát.
- Phát thải $\text{N}_2\text{O}$: IPCC (2006, 2015, 2019a), Snyder và cộng sự (2009), Wang và cộng sự (2011), Chalk và cộng sự (2019), Xing và cộng sự (2021), Zaman và cộng sự (2021) phân tích các yếu tố độ ẩm (WFPS 70–80%), nhiệt độ đất (15–27°C) và hệ số phát thải trực tiếp (EF khoảng 1,0% đối với phân hóa học).
- Mô hình hóa đánh giá tổn thương và rủi ro nước ngầm: Aller và cộng sự (1987) xây dựng mô hình DRASTIC; Foster (1987) phát triển mô hình GOD; Padovani và cộng sự (2002) thiết lập chỉ số nguy cơ nông nghiệp IPNOA. Các nghiên cứu ứng dụng DRASTIC và IPNOA trên thế giới bao gồm Fritch (2000), Pacheco (2015), Barzegar (2019), Alamne (2022), Chakraborty (2022), Yu (2022). Tại Việt Nam, mô hình DRASTIC được áp dụng bởi Bùi Trần Vượng (2004), Nguyễn Đình Tiến (2007), Phạm Quý Nhân (2014).
Khoảng trống nghiên cứu luận án giải quyết
Các nghiên cứu sử dụng DRASTIC tại Việt Nam trước đây chỉ dừng lại ở việc đánh giá tính tổn thương nội tại chung cho toàn khu vực địa chất thủy văn, chưa tích hợp nguồn nguy cơ nông nghiệp cụ thể (IPNOA) từ từng loại hình sử dụng đất trồng cây hàng năm. Bên cạnh đó, chưa có công trình nào kiểm định mức độ tin cậy của mô hình bằng việc đối chiếu chỉ số rủi ro tính toán với nồng độ $\text{NO}_3^-$ thực tế trong tầng chứa nước nông. Ngoài ra, số liệu định lượng trực tiếp đồng thời giữa rửa trôi $\text{NO}_3^-$ và phát thải $\text{N}_2\text{O}$ trên đất trồng cây hàng năm vùng ven biển nhiệt đới còn thiếu hụt.
Cơ sở lý thuyết và phương pháp nghiên cứu
Cơ sở lý thuyết và khung phân tích
Luận án tiếp cận dựa trên lý thuyết chu trình sinh địa hóa nitơ trong hệ thống nông nghiệp, mô hình "lỗ trong đường ống" (hole-in-the-pipe), và nguyên lý đánh giá rủi ro môi trường (Environmental Risk Assessment - ERA). Rủi ro ô nhiễm nước ngầm được xác định theo công thức kết hợp giữa tính tổn thương nội tại của tầng chứa nước và mức độ nguy cơ từ hoạt động bón phân trên bề mặt:
$$\text{Chỉ số Rủi ro (RI)} = \text{Chỉ số Tổn thương (DI)} \times \text{Chỉ số Nguy cơ (HI)}$$
| Thuật ngữ / Khái niệm |
Định nghĩa áp dụng trong luận án |
| Đất trồng cây hàng năm (ĐTCHN) |
Đất gieo trồng, thu hoạch và kết thúc chu kỳ sản xuất không quá 1 năm (theo Thông tư 27/2018/TT-BTNMT), gồm đất lúa và đất trồng cây hàng năm khác. |
| Đất nhiễm mặn |
Đất có nồng độ muối hòa tan cao, độ dẫn điện $\text{EC} > 4\text{ dS/m}$ ở 25°C (tương đương độ mặn khoảng 2,56‰). |
| Nước dưới đất tầng nông |
Tầng chứa nước lỗ hổng trầm tích Holocen ($\text{qh}_2$), là tầng nước ngầm thứ nhất kể từ mặt đất. |
| Hiệu quả sử dụng nitơ (NUE) |
Tỷ lệ phần trăm giữa lượng nitơ tích lũy trong sinh khối thu hoạch so với tổng lượng nitơ cung cấp từ đầu vào. |
Phương pháp nghiên cứu
- Phương pháp điều tra xã hội học nông nghiệp: Thực hiện khảo sát 130 hộ nông dân tại các xã thuộc huyện Giao Thủy thông qua bộ phiếu câu hỏi cấu trúc nhằm thu thập số liệu về diện tích canh tác, cơ cấu giống cây trồng, chế độ luân canh, chủng loại và liều lượng phân đạm, phương thức bón, nguồn nước tưới và phương pháp xử lý phụ phẩm đồng ruộng.
- Phương pháp lấy mẫu và phân tích trong phòng thí nghiệm:
- Mẫu đất: Lấy mẫu định kỳ theo mùa vụ tại 22 xã đại diện cho 4 loại hình sử dụng đất (LUTs): Đất chuyên rau (CV), Đất luân canh rau – lúa (RV), Đất lúa không nhiễm mặn (FRR), Đất lúa nhiễm mặn (SRR). Phân tích các chỉ tiêu: thành phần cơ giới, pH, độ dẫn điện EC, tổng nitơ ($\text{N}_{\text{ts}}$), amoni ($\text{NH}_4^+$-N), nitrat ($\text{NO}_3^-$-N).
- Mẫu nước dưới đất: Lấy mẫu từ các giếng khoan tầng nông ($\text{qh}_2$) tại 22 xã/thị trấn theo các đợt trong năm để phân tích nồng độ $\text{NO}_3^-$.
- Mẫu thực vật: Lấy mẫu hạt giống, cây giống, sản phẩm thu hoạch và phụ phẩm (rơm rạ, thân lá) để phân tích tỷ lệ vật chất khô và hàm lượng nitơ tổng số.
- Phương pháp đánh giá tổn thương nước dưới đất (DRASTIC): Xây dựng 7 lớp bản đồ chuyên đề theo hệ thống Aller và cộng sự (1987) với công thức:
$$\text{DI} = D_r D_w + R_r R_w + A_r A_w + S_r S_w + T_r T_w + I_r I_w + C_r C_w$$
Trong đó: $D$ là độ sâu mực nước ngầm (trọng số $w=5$), $R$ là lượng nước tái nạp thực ($w=4$), $A$ là môi trường tầng chứa nước ($w=3$), $S$ là thành phần cơ giới đất ($w=2$), $T$ là độ dốc địa hình ($w=1$), $I$ là ảnh hưởng của đới thông khí ($w=5$), $C$ là độ thấm thủy lực ($w=3$). Các lớp được số hóa và chồng xếp trên phần mềm ArcGIS.
4. Phương pháp đánh giá nguy cơ ô nhiễm nitrat nông nghiệp (IPNOA): Tính toán chỉ số nguy cơ theo Padovani và cộng sự (2002):
$$\text{HI} = \text{HF} \times \text{CF}_s \times \text{CF}i \times \text{CF}{pa} \times \text{CF}_e$$
Trong đó: $\text{HF}$ là yếu tố nguy cơ từ tải lượng đạm bón (phân khoáng, phân hữu cơ, bùn thải); $\text{CF}_s$ là yếu tố kiểm soát hàm lượng nitơ trong đất; $\text{CF}i$ là kỹ thuật tưới; $\text{CF}{pa}$ là kỹ thuật làm đất và phương pháp bón; $\text{CF}_e$ là điều kiện khí hậu (lượng mưa và nhiệt độ).
5. Phương pháp đo đạc thực địa phát thải $\text{N}_2\text{O}$ và rửa trôi $\text{NO}_3^-$:
- Phát thải $\text{N}_2\text{O}$: Sử dụng phương pháp buồng tĩnh đóng kín (static closed chamber) lắp đặt trực tiếp trên các ô thí nghiệm đại diện cho 4 LUTs với 4 công thức bón đạm: $\text{N}0$ (đối chứng, không bón đạm), $\text{N}{-20}$ (giảm 20% đạm), $\text{N}$ (mức bón thông thường của nông dân), $\text{N}_{+20}$ (tăng 20% đạm). Mẫu khí được hút bằng ống tiêm chuyên dụng và phân tích trên máy sắc ký khí (GC).
- Rửa trôi $\text{NO}_3^-$: Đo lường bằng phương pháp ống bẫy tích lũy tổng hợp SIA (Self-Integrating Accumulators) chôn dưới tầng rễ cây trồng tại thực địa.
- Phương pháp tính toán cân bằng nitơ và xử lý số liệu:
- Lập phương trình cân bằng nitơ tổng thể cho từng loại hình canh tác: $\text{N}{\text{vào}} = \text{N}{\text{ra}} \pm \Delta \text{N}_{\text{tồn lưu}}$.
- Xử lý thống kê, hồi quy đa biến và lựa chọn mô hình trung bình Bayes (bicreg) bằng ngôn ngữ lập trình R; xây dựng bản đồ chuyên đề trên hệ thống thông tin địa lý (ArcGIS).
Nội dung chính theo từng chương
Chương 1. Tổng quan các vấn đề nghiên cứu
Chương 1 hệ thống hóa cơ sở lý luận về chu trình nitơ trong đất nông nghiệp, làm rõ các quá trình chuyển hóa sinh hóa: amin hóa (phân giải protein thành axit amin), amon hóa (chuyển nhóm amin thành $\text{NH}_4^+$), nitrat hóa (oxy hóa $\text{NH}_4^+$ thành $\text{NO}_2^-$ và $\text{NO}_3^-$) và khử nitrat hóa (khử $\text{NO}_3^-$ thành $\text{N}_2\text{O}$ và $\text{N}_2$). Chương này tổng hợp dữ liệu thứ cấp về thực trạng sử dụng phân đạm: mức tiêu thụ đạm cho lúa tại các nước châu Á dao động từ 48–300 kg N/ha/vụ; lượng bón cho rau dao động từ 120–870 kg N/ha.
Phần tổng quan phân tích các hệ quả môi trường gồm: tích lũy nitrat trong rau vượt ngưỡng an toàn; axit hóa đất canh tác (làm tụt giảm 0,3–0,8 pH sau 20 năm thâm canh); thất thoát amoniac qua bay hơi (chiếm 10–14% tổng N mất mát); nguy cơ rửa trôi nitrat gây ô nhiễm nước ngầm tầng nông (trên 19% N bón bị rửa trôi); và phát thải khí nhà kính $\text{N}_2\text{O}$ có khả năng làm ấm toàn cầu gấp 265–298 lần $\text{CO}_2$. Chương 1 cũng trình bày tổng quan điều kiện tự nhiên, đặc điểm địa chất thủy văn tầng chứa nước Holocen ($\text{qh}_2$) tại huyện Giao Thủy, tỉnh Nam Định, từ đó xác định khoảng trống khoa học cần giải quyết.
Chương 2. Đối tượng, phạm vi và phương pháp nghiên cứu
Chương 2 trình bày chi tiết thiết kế nghiên cứu, khung logic và các quy trình kỹ thuật. Phạm vi không gian bao quát 22 xã/thị trấn huyện Giao Thủy; thời gian thu thập dữ liệu và quan trắc trải dài từ năm 2015 đến 2017. Chương này mô tả chi tiết:
- Quy trình điều tra 130 hộ nông dân canh tác trên 4 loại hình sử dụng đất chính: CV, RV, FRR, SRR.
- Quy trình phân tích lý hóa đất, nước ngầm và mẫu sinh khối thực vật theo tiêu chuẩn quốc gia và quốc tế.
- Bảng phân cấp trọng số, điểm số đánh giá tính tổn thương nước ngầm theo mô hình DRASTIC (7 thông số địa chất thủy văn) và chỉ số nguy cơ IPNOA (các yếu tố $\text{HF}, \text{CF}_s, \text{CF}i, \text{CF}{pa}, \text{CF}_e$).
- Quy chuẩn bố trí thí nghiệm buồng tĩnh đo phát thải $\text{N}_2\text{O}$ tại ruộng lúa, ruộng rau, cầu gỗ dẫn ra vị trí lấy khí và phương pháp bẫy tích lũy SIA đo rửa trôi $\text{NO}_3^-$.
- Khung tính toán cân bằng N nông hộ và các phương pháp thống kê trên phần mềm R và ArcGIS.
Chương 3. Kết quả nghiên cứu và thảo luận
Chương 3 trình bày toàn bộ phát hiện thực nghiệm và thảo luận kết quả nghiên cứu qua 4 nhóm nội dung:
1. Hiện trạng sử dụng đất và bón phân đạm tại huyện Giao Thủy
Khảo sát 130 hộ nông dân cho thấy sự phân hóa rõ rệt về lượng phân đạm sử dụng giữa các loại hình canh tác. Lượng đạm vô cơ bón thực tế cho các nhóm cây trồng cao hơn đáng kể so với mức khuyến cáo kỹ thuật:
| Loại hình canh tác / Cây trồng |
Lượng N bón thực tế (kg N/ha/năm) |
Mức khuyến cáo kỹ thuật (kg N/ha/năm) |
Tỷ lệ vượt khuyến cáo |
| Đất chuyên rau (CV) |
532,9 – 1018,1 (TB: ~725) |
250 – 400 |
Vượt 1,5 – 2,5 lần |
| Đất luân canh rau – lúa (RV) |
630,3 – 1085,5 (TB: ~835) |
350 – 500 |
Vượt 1,6 – 2,2 lần |
| Đất lúa không nhiễm mặn (FRR) |
133,3 – 562,9 (TB: ~340) |
200 – 240 (cho 2 vụ) |
Vượt 1,4 – 1,7 lần |
| Đất lúa nhiễm mặn (SRR) |
182,2 – 430,5 (TB: ~335) |
200 – 220 (cho 2 vụ) |
Vượt 1,5 – 1,8 lần |
Hàm lượng $\text{NO}_3^-$-N và $\text{NH}_4^+$-N trong đất biến động mạnh theo thời gian mùa vụ và chế độ bón. Đất chuyên rau và luân canh rau – lúa duy trì tích lũy $\text{NO}_3^-$-N trong tầng đất mặt cao hơn rõ rệt so với đất chuyên lúa nước.
2. Nguy cơ và rủi ro ô nhiễm nitrat trong nước dưới đất tầng nông
- Chỉ số tổn thương nội tại (DI): Tính toán từ mô hình DRASTIC cho thấy các khu vực cồn cát ven biển, nơi có mực nước ngầm nông ($0–1,5\text{ m}$), thành phần cơ giới đất cát pha và độ thấm thủy lực cao, có chỉ số DI thuộc cấp dễ bị tổn thương từ cao đến rất cao.
- Chỉ số nguy cơ nông nghiệp (HI): Mô hình IPNOA ghi nhận các xã có diện tích chuyên rau (Bạch Long, Quất Lâm, Giao Yến) và luân canh rau – lúa (Giao Phong, Giao Long) có chỉ số nguy cơ $\text{HI}$ cao nhất do lượng đạm bón lớn và kỹ thuật tưới rãnh/tưới tràn.
- Bản đồ rủi ro ô nhiễm nitrat ($\text{RI}$): Kết hợp $\text{DI}$ và $\text{HI}$ cho thấy vùng rủi ro cao tập trung tại các xã chuyển đổi mạnh sang rau màu ven biển.
- Kiểm định độ tin cậy mô hình: Quan trắc nồng độ $\text{NO}_3^-$ trong nước dưới đất tầng nông ($\text{qh}_2$) tại 22 xã/thị trấn ghi nhận giá trị từ $0,01\text{ mg/L}$ đến $0,16\text{ mg/L}$ (giá trị trung bình cao nhất tại Giao Phong $0,11\text{ mg/L}$, Bạch Long $0,09\text{ mg/L}$). Kết quả phân tích hồi quy cho thấy mối tương quan thuận rõ rệt giữa chỉ số rủi ro tính toán $\text{RI}$ và nồng độ $\text{NO}_3^-$ thực tế đo được trong nước ngầm tầng nông.
3. Phát thải khí nhà kính nitơ oxit ($\text{N}_2\text{O}$) từ đất trồng cây hàng năm
- Cường độ và tích lũy phát thải: Thí nghiệm đo buồng tĩnh trên 4 LUTs với các công thức bón đạm cho thấy lượng phát thải $\text{N}_2\text{O}$ tích lũy đạt giá trị cao nhất ở đất luân canh rau – lúa (RV) và đất chuyên rau (CV), thấp hơn nhiều ở đất lúa không nhiễm mặn (FRR) và lúa nhiễm mặn (SRR).
- Tác động của mức bón đạm: Khi tăng mức bón từ $\text{N}0 \rightarrow \text{N}{-20} \rightarrow \text{N} \rightarrow \text{N}_{+20}$, lượng $\text{N}2\text{O}$ phát thải tăng theo hàm phi tuyến/tuyến tính thuận. Ở mức bón $\text{N}{+20}$, phát thải $\text{N}_2\text{O}$ gia tăng vượt bậc so với công thức chuẩn $\text{N}$.
- Ảnh hưởng của các yếu tố môi trường đất: Phân tích mô hình hồi quy trung bình Bayes (bicreg trong R) xác định lượng $\text{N}_2\text{O}$ phát thải tương quan thuận chặt chẽ với hàm lượng $\text{NO}_3^-$ trong đất, nhiệt độ đất ($15–26^\circ\text{C}$ kích thích quá trình nitrat/khử nitrat), độ ẩm đất và độ mặn.
4. Đánh giá cân bằng nitơ và đề xuất giải pháp giảm thiểu
- Cân bằng nitơ: Tính toán cân bằng N cho thấy mô hình chuyên rau (CV) và luân canh rau – lúa (RV) có lượng N thặng dư đầu vào rất lớn, dẫn đến thất thoát qua rửa trôi $\text{NO}_3^-$ (lên tới $100–260\text{ kg N/ha/năm}$) và phát thải khí cao, khiến hiệu quả sử dụng nitơ (NUE) bị suy giảm.
- Đề xuất mức bón tối ưu: Luận án xây dựng bảng định mức phân bón khuyến nghị cho từng loại cây trồng (lúa, khoai tây, dưa lê, dưa hấu, bắp cải, củ cải) trên địa bàn huyện Giao Thủy nhằm cắt giảm $20–30%$ lượng đạm dư thừa mà vẫn đảm bảo năng suất cây trồng.
- Giải pháp xử lý tàn dư thực vật: Tái sử dụng rơm rạ và phụ phẩm cây màu thay vì đốt bỏ trên đồng ruộng giúp giảm thiểu phát thải khí nhà kính và bổ sung chất hữu cơ bảo vệ đất.
Kết quả và những đóng góp mới
Đóng góp về mặt lý luận và khoa học
- Thiết lập quy trình tích hợp mô hình tổn thương địa chất thủy văn DRASTIC với chỉ số nguy cơ nông nghiệp IPNOA trên nền GIS để đánh giá rủi ro ô nhiễm nitrat đối với nước dưới đất tầng nông từ nguồn đất trồng cây hàng năm tại vùng ven biển nhiệt đới.
- Kiểm chứng độ tin cậy của mô hình đánh giá rủi ro thông qua việc xác lập mối tương quan thực nghiệm giữa chỉ số rủi ro tính toán ($\text{RI}$) và nồng độ $\text{NO}_3^-$ thực tế trong tầng chứa nước Holocen ($\text{qh}_2$).
- Bổ sung bộ dữ liệu thực nghiệm về hệ số phát thải và động thái phát thải khí nhà kính $\text{N}_2\text{O}$ đo đạc trực tiếp bằng phương pháp buồng tĩnh trên các hệ thống canh tác lúa, luân canh rau – lúa và chuyên rau tại Đồng bằng sông Hồng.
Đóng góp về mặt thực tiễn
- Cung cấp bức tranh hiện trạng định lượng về mức bón phân đạm quá mức trên địa bàn huyện Giao Thủy, xác định các điểm nóng môi trường về rủi ro ô nhiễm nước ngầm và phát thải khí nhà kính.
- Đề xuất bảng định mức bón phân đạm tối ưu cho từng loại cây hàng năm chính, giúp nông dân và cơ quan quản lý nông nghiệp địa phương giảm chi phí đầu vào và giảm thất thoát nitơ ra môi trường.
- Cung cấp dữ liệu cơ sở phục vụ công tác kiểm kê khí nhà kính ngành nông nghiệp theo cam kết Đóng góp do quốc gia tự quyết định (NDC) của Việt Nam.
Hạn chế và hướng nghiên cứu tiếp
- Hạn chế: Các thí nghiệm đo đạc thực địa chuyên sâu và mạng lưới lấy mẫu quan trắc được giới hạn trong phạm vi không gian huyện Giao Thủy, tỉnh Nam Định trong khung thời gian 2015–2017; các thông số mô hình phụ thuộc vào điều kiện thổ nhưỡng phù sa pha cát và địa chất thủy văn ven biển cửa sông Hồng, chưa đánh giá cho các tầng chứa nước sâu hơn.
- Hướng nghiên cứu tiếp: Mở rộng quy mô ứng dụng mô hình tích hợp DRASTIC – IPNOA cho các huyện ven biển khác thuộc vùng Đồng bằng sông Hồng và Đồng bằng sông Cửu Long; thực hiện quan trắc dài hạn dòng phát thải khí nhà kính dưới các kịch bản xâm nhập mặn và biến đổi khí hậu cực đoan; mở rộng đối tượng nghiên cứu sang cây công nghiệp và cây ăn quả ngắn ngày.
Giá trị tham khảo
- Đối với cơ quan quản lý và hoạch định chính sách: Là cơ sở dữ liệu khoa học giúp Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển nông thôn tỉnh Nam Định xây dựng quy hoạch phân vùng nông nghiệp an toàn, quản lý tài nguyên nước ngầm và ban hành hướng dẫn kỹ thuật bón phân hợp lý.
- Đối với nghiên cứu sinh và giảng viên: Là tài liệu tham khảo trong phương pháp tích hợp mô hình không gian GIS (DRASTIC, IPNOA), kỹ thuật đo đạc khí nhà kính bằng buồng tĩnh, phương pháp bẫy tích lũy SIA và kỹ thuật xử lý dữ liệu hồi quy Bayesian (bicreg trong R) trong khoa học môi trường đất và nước.
Câu hỏi thường gặp
1. Bốn loại hình sử dụng đất trồng cây hàng năm (LUTs) chính được luận án nghiên cứu là gì?
Luận án tập trung nghiên cứu 4 loại hình sử dụng đất chính tại huyện Giao Thủy: Đất chuyên rau (CV - Chuyên rau), Đất luân canh rau – lúa (RV - Rau lúa), Đất trồng lúa không nhiễm mặn (FRR - Fresh rice), và Đất trồng lúa bị nhiễm mặn (SRR - Saline rice).
2. Mô hình nào được sử dụng để đánh giá rủi ro ô nhiễm nitrat tới nước dưới đất tầng nông?
Luận án kết hợp mô hình DRASTIC (đánh giá tính dễ bị tổn thương nội tại của tầng chứa nước gồm 7 thông số: Độ sâu mực nước, Lượng nước tái nạp, Tầng chứa nước, Thành phần cơ giới đất, Địa hình, Đới thông khí, Độ thấm thủy lực) với mô hình IPNOA (đánh giá chỉ số nguy cơ ô nhiễm nitrat nông nghiệp từ tải lượng phân bón và các yếu tố kiểm soát canh tác) để tính chỉ số rủi ro $\text{RI} = \text{DI} \times \text{HI}$ trên nền tảng ArcGIS.
3. Thực trạng bón phân đạm trên đất chuyên rau và luân canh rau – lúa tại Giao Thủy có đặc điểm gì?
Lượng phân đạm bón thực tế trên đất chuyên rau (CV) dao động từ 532,9 đến 1018,1 kg N/ha/năm và đất luân canh rau – lúa (RV) dao động từ 630,3 đến 1085,5 kg N/ha/năm, vượt từ 1,5 đến 2,5 lần so với mức khuyến cáo kỹ thuật, dẫn đến lượng thặng dư nitơ trong đất cao và làm gia tăng thất thoát ra môi trường.
4. Phương pháp nào được sử dụng để đo đạc phát thải $\text{N}_2\text{O}$ và rửa trôi $\text{NO}_3^-$ tại thực địa?
Phát thải khí $\text{N}2\text{O}$ được đo đạc trực tiếp bằng phương pháp buồng tĩnh đóng kín (static closed chamber) trên các ô thí nghiệm với các công thức bón phân $\text{N}0, \text{N}{-20}, \text{N}, \text{N}{+20}$ và phân tích qua sắc ký khí. Lượng $\text{NO}_3^-$ rửa trôi được đo bằng phương pháp bẫy tích lũy tổng hợp SIA (Self-Integrating Accumulators) chôn dưới tầng rễ cây trồng.
5. Những phát hiện chính về cân bằng nitơ và hiệu quả sử dụng phân đạm (NUE) trong luận án là gì?
Các loại hình canh tác rau màu (CV và RV) có lượng nitơ đầu vào vượt xa nhu cầu cây trồng, dẫn đến hiệu quả sử dụng nitơ (NUE) thấp, lượng nitrat rửa trôi ra tầng nước ngầm nông và lượng phát thải khí $\text{N}_2\text{O}$ cao hơn nhiều lần so với đất chuyên lúa (FRR và SRR). Việc cắt giảm 20% lượng đạm bón không làm giảm năng suất mà giúp giảm đáng kể phát thải khí nhà kính và nguy cơ ô nhiễm nguồn nước.
Kết luận
Luận án tiến sĩ của tác giả Lý Thị Thu Hà đã hoàn thành các mục tiêu nghiên cứu thông qua việc đánh giá định lượng hiện trạng bón phân đạm, rủi ro rửa trôi $\text{NO}_3^-$ tới nước dưới đất tầng nông ($\text{qh}_2$) và mức phát thải khí nhà kính $\text{N}_2\text{O}$ trên các loại hình đất trồng cây hàng năm tại huyện Giao Thủy, tỉnh Nam Định. Bằng việc kết hợp phương pháp mô hình hóa không gian DRASTIC – IPNOA với quan trắc thực địa buồng tĩnh và bẫy tích lũy SIA, công trình đã xác lập cơ sở khoa học về mối liên hệ giữa thâm canh phân đạm và suy thoái chất lượng môi trường ven biển. Kết quả nghiên cứu cung cấp bộ số liệu thực nghiệm tin cậy và các giải pháp định mức phân bón tối ưu phục vụ công tác quản lý tài nguyên đất, bảo vệ nguồn nước ngầm và giảm thiểu phát thải khí nhà kính trong sản xuất nông nghiệp.