Tổng quan về luận án

Biến đổi khí hậu toàn cầu và sự gia tăng nồng độ khí nhà kính (KNK) trong khí quyển là một trong những thách thức môi trường cấp bách nhất của thế kỷ XXI. Tại Việt Nam, sản xuất lúa nước đóng vai trò trụ cột đảm bảo an ninh lương thực quốc gia với diện tích gieo trồng trên 7,7 triệu ha và sản lượng xấp xỉ 44 triệu tấn/năm. Tuy nhiên, theo Báo cáo cập nhật hai năm một lần lần thứ II (BUR2) của Bộ Tài nguyên và Môi trường (2017), canh tác lúa là nguồn phát thải KNK lớn nhất trong nông nghiệp, tạo ra 44,7 triệu tấn $\text{CO}_2\text{tđ}$, chiếm 50,3% tổng phát thải KNK ngành nông nghiệp và 17,2% tổng phát thải quốc gia. Trong đó, trên 85% diện tích lúa áp dụng chế độ tưới ngập thường xuyên – môi trường tối ưu cho quá trình sinh khí mê-tan ($\text{CH}_4$).

Khoảng trống nghiên cứu (research gap) lớn của khoa học môi trường đất và khí hậu tại Việt Nam là sự mất cân đối vùng miền và giới hạn về phương pháp lượng hóa. Hầu hết các nghiên cứu định lượng phát thải KNK đồng ruộng trước đây tập trung chủ yếu tại Đồng bằng sông Hồng (Nguyễn Hữu Thành và nnk, 2011; Vũ Thắng và nnk) và Đồng bằng sông Cửu Long (Wassmann và nnk, 2000). Vùng duyên hải Nam Trung Bộ, đặc biệt là lưu vực sông Vu Gia – Thu Bồn (tỉnh Quảng Nam), với chế độ nhiệt ẩm biến động mạnh, địa hình chia cắt phức tạp và đặc tính đất phù sa chua, phù sa glây đặc thù, hầu như chưa có số liệu quan trắc hệ thống và mô hình hóa phát thải. Hơn nữa, việc kiểm kê KNK quốc gia vẫn phụ thuộc vào phương pháp Bậc 1 (Tier 1) của IPCC với hệ số mặc định toàn cầu, dẫn đến sai số rất lớn.

Để giải quyết khoảng trống này, luận án xác định hai câu hỏi nghiên cứu cốt lõi:

  1. Động thái phát thải $\text{CH}_4$ và $\text{N}_2\text{O}$ từ đất lúa lưu vực sông Vu Gia – Thu Bồn biến thiên như thế nào dưới tác động tương hỗ của chế độ quản lý nước (ngập liên tục vs. tưới tiết kiệm), mức bón phân và các thông số hóa lý đất ($\text{Eh}$, $\text{pH}$, nhiệt độ)?
  2. Mô hình sinh địa hóa DNDC (DeNitrification-DeComposition) có khả năng mô phỏng chuẩn xác phát thải KNK tại địa bàn nghiên cứu hay không, và việc tích hợp DNDC-GIS có thể lượng hóa tiềm năng giảm phát thải quy mô vùng theo các kịch bản tưới như thế nào?

Luận án kiểm định hai giả thuyết khoa học:

  • Giả thuyết $H_1$: Kỹ thuật tưới rút nước giữa vụ (nông lộ phơi/AWD) làm tăng thế oxy hóa khử ($\text{Eh}$) của đất, ức chế vi khuẩn sinh mê-tan, từ đó cắt giảm đáng kể phát thải $\text{CH}_4$ mà không làm bùng phát $\text{N}_2\text{O}$ vượt mức bù trừ tiềm năng ấm lên toàn cầu (GWP), đồng thời duy trì ổn định năng suất lúa.
  • Giả thuyết $H_2$: Mô hình DNDC sau khi được hiệu chỉnh và kiểm định bằng dữ liệu vi khí hậu, thổ nhưỡng và tập quán địa phương có thể mô phỏng chính xác phát thải $\text{CH}_4$ và $\text{N}_2\text{O}$, làm nền tảng xây dựng hệ số phát thải riêng (Tier 2) và bản đồ phát thải số hóa (Tier 3).

Nghiên cứu được triển khai trên quy mô 2 trắc nghiệm đồng ruộng đa yếu tố qua 2 vụ canh tác (Đông Xuân và Hè Thu) tại huyện Đại Lộc và Duy Xuyên, kết hợp điều tra 155 nông hộ và mô hình hóa toàn diện trên hơn 41.000 ha đất lúa toàn lưu vực sông Vu Gia – Thu Bồn.

Literature Review và Positioning

Nghiên cứu phát thải KNK từ đất ngập nước trên thế giới bắt đầu từ các công trình kinh điển của Koyama (1963) tại Nhật Bản, Ehhalt & Schmidt (1978), Cicerone & Shetter (1981) tại Hoa Kỳ, và Schutz (1989) tại Ý. Các nghiên cứu quốc tế đã xác lập nền tảng động học về sự hình thành $\text{CH}_4$ (methanogenesis) trong điều kiện kỵ khí tuyệt đối ($\text{Eh} < -150\text{ mV}$) bởi vi khuẩn cổ sinh mê-tan (Archaea methanogens) sử dụng chất nền acetate hoặc $\text{H}_2/\text{CO}_2$ (Conrad, 1996; Neue và nnk, 1994, 1998). Đối với $\text{N}_2\text{O}$, quá trình phát thải được điều khiển bởi mô hình rò rỉ ống dẫn (Hole-in-the-pipe model) qua hai cơ chế sinh hóa: nitrat hóa (nitrification) trong điều kiện hiếu khí và phản nitrat hóa (denitrification) trong điều kiện kỵ khí (Firestone & Davidson, 1989; Bouwman, 1998; Majumdar, 2003).

Trong y văn tồn tại những tranh luận học thuật sâu sắc về mối quan hệ đánh đổi (trade-off) giữa $\text{CH}_4$ và $\text{N}_2\text{O}$. Trường phái thứ nhất (Smith và nnk, 2007; Yan và nnk, 2003) lập luận rằng việc tháo cạn ruộng để giảm $\text{CH}_4$ sẽ kích hoạt quá trình khoáng hóa nitơ và chuyển hóa nitrat, dẫn đến hiện tượng "bùng nổ phát thải $\text{N}_2\text{O}$" (nitrous oxide pulse), làm triệt tiêu hiệu quả giảm GWP. Ngược lại, trường phái thứ hai (Wassmann và nnk, 2000, 2006; Li và nnk, 2004, 2006) khẳng định rằng nếu kiểm soát chặt chẽ thời điểm rút nước trùng với giai đoạn đẻ nhánh rộ và bón phân đạm cân đối, lượng $\text{CH}_4$ giảm được sẽ vượt trội hoàn toàn lượng $\text{N}_2\text{O}$ gia tăng, do tỷ trọng phát thải $\text{N}_2\text{O}$ trong đất lúa nước nhiệt đới thường chiếm dưới 1% tổng nitơ đầu vào.

So sánh với hai nghiên cứu quốc tế điển hình:

  1. Nghiên cứu của Yagi và nnk (1996) tại Nhật Bản chỉ ra việc rút nước giữa vụ giúp giảm 40-45% lượng phát thải $\text{CH}_4$ nhưng chủ yếu khảo sát trên đất lúa ôn đới với chu kỳ phân hủy chất hữu cơ chậm.
  2. Nghiên cứu của Corton và nnk (2000) tại Viện Lúa Quốc tế (IRRI - Philippines) ghi nhận giảm 50% $\text{CH}_4$ khi áp dụng kỹ thuật tưới khô ướt xen kẽ (AWD) nhưng chưa tích hợp mô hình hóa sinh địa hóa quy mô lưu vực đa dạng thổ nhưỡng.

Luận án định vị tiên phong bằng cách kết hợp trắc nghiệm đo đạc buồng kín tĩnh theo thời gian thực với mô hình hóa sinh địa hóa DNDC và hệ thông tin địa lý GIS cho một lưu vực sông miền Trung Việt Nam, bắc cầu nối từ thực nghiệm vi mô đến chính sách vĩ mô.

Đó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à làm sâu sắc thêm Lý thuyết Động học Oxy hóa Khử Đất ngập nước (Redox Chemistry of Submerged Soils) của Ponnamperuma (1972, 1978, 1985) và Patrick & Reddy (1983). Nghiên cứu cung cấp bằng chứng thực nghiệm khẳng định sự phân tầng thế điện thế $\text{Eh}$ và ngưỡng phát thải $\text{CH}_4$ trong đất phù sa Nam Trung Bộ: khi $\text{Eh}$ hạ dưới $-150\text{ mV}$ đến $-239\text{ mV}$, cường độ phát thải $\text{CH}_4$ tăng vọt đạt đỉnh (lên tới 41,88 $\text{mg CH}_4/\text{m}^2/\text{giờ}$); khi tháo nước đưa $\text{Eh}$ lên vùng dương ($+75\text{ mV}$ đến $+170\text{ mV}$), quá trình sinh mê-tan bị triệt tiêu tức thì.

Đồng thời, công trình đóng góp vào Lý thuyết Quá trình Sinh địa hóa Phức hợp (Biogeochemical Process-based Theory) được tích hợp trong mô hình DNDC của Changsheng Li (1992, 2000, 2006). Luận án chứng minh tính thích ứng của mô hình DNDC đối với hệ sinh thái nông nghiệp nhiệt đới gió mùa có biên độ nhiệt và chế độ tưới phân tán, kiểm thực định chế phân rã chất hữu cơ hòa tan (DOC) và động thái amoni/nitrat trong tầng đất canh tác nhiệt đới.

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

Khung phân tích của luận án là sự hợp nhất của ba trụ cột lý thuyết và công nghệ:

  1. Lý thuyết Động học Vi sinh vật Đất (Microbial Kinetics of Methanogenesis & Denitrification).
  2. Lý thuyết Cân bằng Phát thải và Tiềm năng Ấm lên Toàn cầu ($\text{GWP} = 25 \times \text{CH}_4 + 298 \times \text{N}_2\text{O}$ trên chu kỳ 100 năm theo IPCC AR4).
  3. Lý thuyết Phân tích Không gian Địa lý (Spatial Environmental Modeling via GIS).

Điều kiện biên (boundary conditions) của khung phân tích được xác định rõ: áp dụng cho hệ thống canh tác lúa nước 2 vụ/năm vùng nhiệt đới ven biển, với 5 nhóm đất chính (đất cát biển, đất phù sa chua, phù sa ít chua, phù sa glây, và đất phù sa nhiễm mặn/phèn), dung trọng đất từ 1,02 đến 1,38 $\text{g/cm}^3$, hàm lượng cacbon hữu cơ tầng mặt ($\text{SOC}$) dao động từ 0,5% đến 2,5%.

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

Thiết kế nghiên cứu

Luận án theo đuổi thế giới quan thực chứng (positivism) kết hợp trường phái hiện thực phản biện (critical realism). Thiết kế nghiên cứu đa tầng (multi-level research design) tích hợp chặt chẽ giữa quan sát thực địa, thực nghiệm đối chứng ngẫu nhiên và mô phỏng số học:

  • Cấp độ ô thí nghiệm (Plot level): Bố trí thí nghiệm khối ngẫu nhiên hoàn chỉnh (RCBD) với 3 lần lặp lại tại 2 vùng sinh thái đất điển hình: Điểm 1 tại Đại Quang - Đại Lộc (đất phù sa cổ/phù sa ít chua) và Điểm 2 tại Nam Phước - Duy Xuyên (đất phù sa glây giàu hữu cơ).
  • Cấp độ nông hộ (Household level): Khảo sát 155 hộ nông dân về quy trình quản lý phân bón, nước và tàn dư rơm rạ thông qua hệ cơ sở dữ liệu Microsoft Access chuyên biệt.
  • Cấp độ vùng (Regional level): Mô phỏng lưới không gian toàn lưu vực sông Vu Gia – Thu Bồn trên nền bản đồ thổ nhưỡng và khí hậu số hóa.

Quy trình nghiên cứu rigorous

Quy trình thu thập và phân tích mẫu tuân thủ nghiêm ngặt chuẩn quốc tế của IRRI và mạng lưới GRACEnet:

  • Đo KNK: Sử dụng phương pháp buồng kín tĩnh (closed static chamber) bằng vật liệu acrylic, có trang bị quạt đảo khí, nhiệt kế điện tử và ống thông áp. Mẫu khí được hút bằng bơm tiêm chuyên dụng tại các mốc thời gian 0, 10, 20 và 30 phút sau khi đóng buồng, tiến hành từ 8h00 đến 10h30 sáng để đảm bảo tính đại diện nhiệt độ trung bình ngày. Tần suất đo 7-10 ngày/lần, tăng dày trong các đợt rút nước và bón thúc.
  • Phân tích sắc ký khí: Nồng độ $\text{CH}_4$ được định lượng trên máy sắc ký khí (Gas Chromatography - GC) với đầu dò ion hóa ngọn lửa (FID), nồng độ $\text{N}_2\text{O}$ được định lượng bằng đầu dò cộng kết điện tử ($\mu\text{ECD}$). Đường chuẩn khí tiêu chuẩn có hệ số xác định $R^2 > 0,999$.
  • Đo thông số môi trường đất: Điện thế oxy hóa khử ($\text{Eh}$) được đo trực tiếp bằng điện cực bạch kim (Pt electrode) kết hợp điện cực so sánh calomel/$\text{Ag/AgCl}$ cắm cố định tại độ sâu tầng rễ (0-10 cm và 10-20 cm). Đo $\text{pH}$ đất và nhiệt độ tầng đất tức thời tại mỗi thời điểm lấy mẫu khí.

Data và phân tích

Số liệu khí tượng trích xuất từ trạm khí tượng Tam Kỳ (nhiệt độ tối cao, tối thấp, lượng mưa, bức xạ mặt trời hàng ngày). Mô hình DNDC phiên bản 9.5 được hiệu chỉnh (calibration) và kiểm định (validation) thông qua các chỉ số thống kê khắt khe: Hệ số hiệu quả mô hình Nash-Sutcliffe ($EF$), Sai số bình phương trung bình căn nguyên ($RMSE$), và Hệ số xác định ($R^2$). Kiểm tra độ nhạy (sensitivity analysis) được thực hiện đa biến đối với 12 thông số đầu vào gồm dung trọng đất, $\text{SOC}$, $\text{pH}$, lượng đạm bón, độ sâu cày xới và thời gian tháo nước.

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

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

Luận án mang lại 5 phát hiện đột phá có giá trị thực nghiệm cao:

  1. Quy luật động thái phát thải $\text{CH}_4$ theo sinh trưởng: Cường độ phát thải $\text{CH}_4$ rất thấp ở giai đoạn đầu sau sạ (13-18 NSS: 3,05 - 6,41 $\text{mg/m}^2/\text{giờ}$), tăng vọt và đạt đỉnh tối đa ở giai đoạn đẻ nhánh rộ đến phân hóa đòng - đứng cái (28-48 NSS: 20,65 - 41,88 $\text{mg/m}^2/\text{giờ}$ ở công thức ngập thường xuyên TN), sau đó suy giảm nhanh khi lúa vào chắc và chín (dưới 1,0 $\text{mg/m}^2/\text{giờ}$).
  2. Hiệu ứng cắt giảm phát thải vượt trội của chế độ tưới rút nước giữa vụ: So với công thức tưới ngập liên tục (TN), chế độ tưới tiết kiệm rút nước 1 lần giữa vụ (1RN) và tưới khô ướt xen kẽ (TTK/AWD) làm giảm từ 28,6% đến 46,3% tổng lượng phát thải $\text{CH}_4$ tích lũy toàn vụ. Cụ thể tại Điểm 2 vụ Hè Thu, phát thải $\text{CH}_4$ giảm từ đỉnh 41,88 xuống còn 25,24 $\text{mg/m}^2/\text{giờ}$ trong đợt tháo cạn.
  3. Mối tương quan nghịch phi tuyến giữa $\text{CH}_4$ và $\text{Eh}$ đất: Phân tích hồi quy xác lập tương quan nghịch rất chặt chẽ giữa cường độ phát thải $\text{CH}_4$ và điện thế khử $\text{Eh}$ ($R^2$ từ $-0,75$ đến $-0,88$, $p < 0,01$). Khi tháo nước phơi ruộng, $\text{Eh}$ tăng từ $-234\text{ mV}$ lên mức dương $+85\text{ mV}$ đến $+168\text{ mV}$, vô hiệu hóa hoàn toàn điều kiện sinh mê-tan.
  4. Phát thải $\text{N}_2\text{O}$ biến thiên ngược chiều nhưng không làm tăng tổng GWP: Cường độ phát thải $\text{N}_2\text{O}$ ở điều kiện ngập nước duy trì mức nền rất thấp (0,004 - 0,020 $\text{mg/m}^2/\text{giờ}$), nhưng xuất hiện các đỉnh phát thải ngắn khi rút cạn nước và bón phân đạm (đạt 0,035 - 0,051 $\text{mg/m}^2/\text{giờ}$). Tuy nhiên, do lượng phát thải $\text{CH}_4$ chiếm tỷ trọng áp đảo, tổng tiềm năng gây ấm toàn cầu ($\text{GWP}$) ở các công thức tưới tiết kiệm vẫn giảm từ 24,1% đến 38,5% so với tưới ngập liên tục.
  5. Năng suất lúa không bị suy giảm khi tưới tiết kiệm: Năng suất thực thu ở nghiệm thức rút nước 1 lần giữa vụ (1RN) tương đương hoặc tăng nhẹ 3,2 - 5,8% so với tưới ngập truyền thống nhờ bộ rễ ăn sâu, giảm hiện tượng đổ ngã và ngộ độc hữu cơ rễ lúa.

Implications đa chiều

  • Ý nghĩa lý thuyết: Cung cấp bộ tham số thực nghiệm nhiệt đới chính xác về hệ số chuyển hóa C-N trong môi trường đất lúa phù sa, khẳng định tính đúng đắn của việc mô hình hóa sinh địa hóa đối với các vùng khí hậu biến động mạnh.
  • Ý nghĩa phương pháp luận: Chuẩn hóa quy trình tích hợp DNDC-GIS trong việc nâng bậc kiểm kê KNK từ Tier 1 lên Tier 2 và Tier 3, khắc phục hạn chế tốn kém chi phí của các trạm đo đạc vi khí tượng Eddy Covariance.
  • Ý nghĩa thực tiễn & chính sách: Luận án đề xuất hệ số phát thải nền $\text{EF}_i = 1,48\text{ kg CH}_4/\text{ha/ngày}$ cho đất lúa ngập liên tục không bón chất hữu cơ tại Quảng Nam, đồng thời xây dựng hệ số điều chỉnh chế độ tưới $\text{SF}_w$ chuẩn xác (TN = 1,0; 1RN = 0,71; 2RN = 0,58). Đây là căn cứ khoa học trực tiếp để Sở Nông nghiệp và PTNT tỉnh Quảng Nam và Bộ Nông nghiệp & PTNT phê duyệt lộ trình mở rộng kỹ thuật tưới 1RN trên 75% diện tích đất lúa chủ động tưới vào năm 2025.

Limitations và Future Research

Luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu:

  1. Độ phân giải cơ sở dữ liệu không gian: Dữ liệu thổ nhưỡng đầu vào ở cấp lưu vực được tổng hợp từ bản đồ tỷ lệ 1:50.000, do đó có sự khái quát hóa nhất định về các biến thiên vi địa hình và độ phì tầng đất mặt cục bộ.
  2. Khả năng mô phỏng đỉnh $\text{N}_2\text{O}$ của DNDC: Mô hình DNDC mô phỏng cực kỳ chuẩn xác động thái $\text{CH}_4$ ($R^2 = 0,78 - 0,91$), nhưng có xu hướng đánh giá thấp hơn thực tế các xung phát thải $\text{N}_2\text{O}$ ngắn hạn ngay sau khi ruộng được cấp nước trở lại sau đợt phơi khô đột ngột.
  3. Phạm vi khảo sát chất hữu cơ: Nghiên cứu tập trung vào biện pháp bón phân khoáng vô cơ kết hợp tàn dư rơm rạ truyền thống, chưa mở rộng sang các dạng phân bón hữu cơ chế biến vi sinh hoặc than sinh học (biochar).
  4. Chuỗi số liệu thực địa: Số liệu trắc nghiệm đồng ruộng được đo đạc liên tục qua 2 vụ lúa chính; cần tiếp tục quan trắc dài hạn 5-10 năm để đánh giá tích lũy cacbon hữu cơ trong đất ($\text{SOC}$) dưới tác động lâu dài của hạn hán và xâm nhập mặn.

Chương trình nghiên cứu tương lai cần hướng tới:

  • Ứng dụng kỹ thuật đồng vị bền ($^{13}\text{C}$ và $^{15}\text{N}$) để bóc tách chính xác tỷ lệ phát thải KNK từ nguồn phân bón nhân tạo so với nguồn gốc hữu cơ nội sinh của đất.
  • Tích hợp công nghệ cảm biến IoT đo $\text{Eh}$ và mực nước tự động thời gian thực để điều khiển hệ thống tưới tiêu thông minh thích ứng biến đổi khí hậu.

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

Luận án tạo ra tác động lan tỏa mạnh mẽ trên nhiều phương diện:

  • Học thuật: Cung cấp bộ dữ liệu thực nghiệm chuẩn hóa đầu tiên về phát thải KNK đất lúa tại miền Trung Việt Nam, mở rộng chỉ số trích dẫn cho các nghiên cứu mô hình hóa nông nghiệp Đông Nam Á trên các tạp chí quốc tế chuyên ngành (Agriculture, Ecosystems & Environment; Geoderma).
  • Ngành nông nghiệp & thủy lợi: Chuyển đổi tư duy quản lý từ "tưới ngập đệm nước sâu" sang "quản lý nước tiết kiệm thông minh", giúp tiết kiệm từ 15% đến 25% khối lượng nước tưới mặt ruộng, giảm xung đột nguồn nước giữa sản xuất nông nghiệp và thủy điện trong lưu vực sông Vu Gia – Thu Bồn vào mùa cạn.
  • Môi trường & Xã hội: Nếu áp dụng triệt để chế độ tưới rút nước 1 lần giữa vụ (1RN) cho toàn bộ diện tích lúa có tưới của tỉnh Quảng Nam, ước tính sẽ cắt giảm hàng chục nghìn tấn $\text{CO}_2\text{tđ}$ mỗi năm, đóng góp thiết thực vào Cam kết Đóng góp do quốc gia tự quyết định (NDC) của Việt Nam thực hiện Thỏa thuận Paris.

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

  1. Nghiên cứu sinh & Giới học thuật: Tiếp cận phương pháp luận chuẩn xác về vận hành mô hình DNDC, quy trình đo buồng kín khí sắc ký và kỹ thuật phân tích thống kê tương quan sinh địa hóa phức hợp.
  2. Cơ quan Quản lý Nhà nước (Bộ TN&MT, Bộ NN&PTNT, Chi cục Trồng trọt & BVTV): Sở hữu bộ hệ số phát thải đặc thù vùng ($\text{EF}_i$, $\text{SF}_w$) phục vụ báo cáo kiểm kê KNK quốc gia và xây dựng quy hoạch nông nghiệp carbon thấp.
  3. Các Công ty Khai thác Công trình Thủy lợi: Cơ sở dữ liệu khoa học để tái cơ cấu lịch đóng mở cống tưới, điều tiết nguồn nước liên hồ chứa hiệu quả.
  4. Hợp tác xã Nông nghiệp & Nông dân: Nâng cao năng suất lúa, tiết kiệm chi phí nhiên liệu bơm tưới và tạo tiền đề tham gia thị trường tín chỉ carbon nông nghiệp trong tương lai.

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ì?

Luận án đã mở rộng Lý thuyết Động học Oxy hóa Khử Đất lúa nhiệt đới của Ponnamperuma thông qua việc lượng hóa toán học ngưỡng kích hoạt và triệt tiêu phát thải $\text{CH}_4$ trong mối liên kết động với điện thế $\text{Eh}$ và $\text{pH}$ trên 5 nhóm đất phù sa miền Trung. Công trình chứng minh rằng trạng thái chuyển đổi oxy hóa khử nhanh do rút nước ngắn ngày phá vỡ hoàn toàn chuỗi phản ứng sinh mê-tan của vi khuẩn cổ mà không tạo ra hiện tượng thất thoát nitơ dạng $\text{N}_2\text{O}$ nghiêm trọng.

2. Điểm mới về phương pháp luận so với các nghiên cứu trước đây?

So với nghiên cứu thực nghiệm đơn thuần của Nguyễn Hữu Thành (2011) hay mô hình hóa thuần túy của Vũ Mạnh Quyết (2006), luận án đã thực hiện phương pháp tiếp cận đa bậc hoàn chỉnh: kết hợp thực nghiệm buồng kín sắc ký khí chuẩn hóa $\to$ hiệu chỉnh thông số mô hình DNDC theo đặc thù thổ nhưỡng địa phương $\to$ tích hợp GIS mở rộng quy mô toàn lưu vực. Đây là bước đột phá kỹ thuật nâng mức kiểm kê phát thải từ Tier 1 lên Tier 2/Tier 3 cho một lưu vực sông cụ thể.

3. Phát hiện nào gây bất ngờ nhất từ số liệu thực nghiệm?

Kết quả bất ngờ nhất là tại vùng đất phù sa chua huyện Duy Xuyên, việc áp dụng chế độ tưới rút nước giữa vụ làm gia tăng nhẹ năng suất lúa thực tế (tăng 3,2 - 5,8%) thay vì gây rủi ro giảm năng suất do thiếu hụt nước như lo ngại ban đầu của nông dân. Nguyên nhân sinh học là việc phơi ruộng giúp giải phóng khí độc ($\text{H}_2\text{S}$, axit hữu cơ) tích tụ trong môi trường yếm khí, kích thích hệ rễ lúa phát triển phân nhánh ăn sâu.

4. Luận án có cung cấp quy trình nhân bản nghiên cứu (Replication Protocol) không?

Có. Luận án chi tiết hóa toàn bộ giao thức thực địa: bản vẽ thiết kế buồng kín tĩnh, mã nguồn bảng nhập liệu Microsoft Access (FHHS_DataEntry), tập tin tham số hóa đất đai và khí tượng đầu vào mô hình DNDC (Tamky2011.txt, các file .dnd cho từng nhóm đất), cùng quy trình phân tích sắc ký khí kiểm soát chuẩn.

5. Định hướng chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Xây dựng mạng lưới tháp quan trắc vi khí tượng tự động (Eddy Covariance towers) trên các lưu vực sông trọng điểm miền Trung; phát triển thuật toán vi thám vệ tinh kết hợp AI để giám sát chế độ ngập ruộng lúa thời gian thực; và đánh giá toàn diện chu trình vòng đời (LCA) nông nghiệp tuần hoàn kết hợp biochar và năng lượng sinh học.

Kết luận

  1. Xác lập thành công bộ cơ sở dữ liệu khoa học toàn diện về phát thải KNK ($\text{CH}_4$, $\text{N}_2\text{O}$) và đặc tính hóa lý môi trường đất lúa lưu vực sông Vu Gia – Thu Bồn qua các mùa vụ và chế độ canh tác đối chứng.
  2. Chứng minh thực nghiệm chế độ tưới tiết kiệm rút nước giữa vụ giúp giảm từ 28,6% đến 46,3% lượng phát thải $\text{CH}_4$, cắt giảm 24,1% - 38,5% tổng tiềm năng ấm lên toàn cầu ($\text{GWP}$) mà vẫn duy trì ổn định và gia tăng năng suất lúa.
  3. Hiệu chỉnh và kiểm định thành công mô hình sinh địa hóa DNDC với độ tin cậy thống kê cao ($R^2 > 0,80$), khẳng định tính khả thi của mô hình trong điều kiện nhiệt đới miền Trung Việt Nam.
  4. Xây dựng bản đồ phân bố không gian phát thải $\text{CH}_4$, $\text{N}_2\text{O}$ và tổng $\text{GWP}$ quy mô vùng bằng công nghệ GIS, trực quan hóa tiềm năng giảm thiểu phát thải của từng huyện trong lưu vực.
  5. Đề xuất hệ số phát thải nền $\text{EF}_i = 1,48\text{ kg CH}_4/\text{ha/ngày}$ và bộ hệ số tỷ lệ chế độ tưới $\text{SF}_w$ đặc thù, đóng góp trực tiếp vào việc hoàn thiện hệ thống kiểm kê KNK Bậc 2 và Bậc 3 cấp quốc gia.
  6. Mở ra 3 hướng nghiên cứu mới về công nghệ giám sát nông nghiệp thông minh, thị trường tín chỉ carbon lúa nước và tối ưu hóa chu trình sinh địa hóa C-N thích ứng với biến đổi khí hậu toàn cầu.