Giới thiệu dự án

Sản xuất lúa gạo (Oryza sativa L.) giữ vị trí chiến lược trong việc bảo đảm an ninh lương thực toàn cầu và ổn định kinh tế - xã hội tại Việt Nam. Theo thống kê từ Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc (FAO), lúa gạo là nguồn lương thực chủ lực nuôi sống hơn 40% dân số thế giới. Tại Việt Nam, ngành lúa gạo duy trì sản lượng bình quân hàng năm trên 42,69 triệu tấn (Tổng cục Thống kê, 2020), đóng góp kim ngạch xuất khẩu 3,07 tỷ USD với khối lượng 6,15 triệu tấn gạo chất lượng cao. Tuy nhiên, hiệu quả sử dụng phân bón hóa học trong canh tác lúa tại Đồng bằng sông Hồng và Đồng bằng sông Cửu Long vẫn còn ở mức thấp; ước tính có tới 1,2 – 1,3 triệu tấn đạm (trên tổng số 2 triệu tấn phân bón bón vào đất mỗi năm) bị thất thoát qua các con đường bốc hơi amoniac ($NH_3$), rửa trôi nitrat ($NO_3^-$) và thẩm thấu sâu, gây lãng phí kinh tế và gia tăng phát thải khí nhà kính ($N_2O, CH_4$).

Vấn đề cốt lõi đặt ra trong thực tiễn sản xuất là tập quán bón thừa phân đạm vô cơ của nông dân nhằm tối đa hóa sinh trưởng thân lá, dẫn đến hiện tượng lúa lốp đổ, giảm tỷ lệ tích lũy chất khô vào hạt, kéo dài thời gian sinh trưởng sinh dưỡng và làm bùng phát các đối tượng dịch hại nguy hiểm như bệnh đạo ôn (Magnaporthe oryzae), bệnh bạc lá (Xanthomonas oryzae) và sâu cuốn lá nhỏ (Cnaphalocrocis medinalis). Đối với giống lúa thuần cảm ôn mới chọn tạo TBR125 của Công ty Cổ phần Giống Cây trồng Thái Bình (ThaiBinh Seed), việc thiếu hụt một quy trình định lượng phân đạm chuẩn xác theo từng giai đoạn phát triển là rào cản lớn khiến giống chưa thể phát huy tối đa tiềm năng năng suất (vốn có thể đạt 70 – 90 tạ/ha).

Dự án khóa luận tốt nghiệp tập trung giải quyết các mục tiêu cụ thể sau:

  1. Xác định động thái tăng trưởng chiều cao cây, tốc độ ra lá và khả năng đẻ nhánh của giống lúa tẻ TBR125 qua 4 mức đạm bón trong vụ Xuân 2021.
  2. Định lượng các chỉ tiêu sinh lý then chốt bao gồm chỉ số diện tích lá (Leaf Area Index - LAI), hàm lượng diệp lục qua chỉ số SPAD và động thái tích lũy khối lượng chất khô (Dry Matter - DM) tại 3 thời kỳ xung yếu: đẻ nhánh rộ, trỗ 80% và chín sáp.
  3. Đánh giá mức độ phát sinh và gây hại của các loại sâu bệnh chính tương ứng với từng mức độ cung cấp dinh dưỡng đạm.
  4. Xác lập mức bón đạm tối ưu nhằm tối đa hóa năng suất lý thuyết (NSLT), năng suất thực thu (NSTT) và nâng cao hiệu quả kinh tế trên đơn vị diện tích canh tác.

Giải pháp được tiếp cận thông qua phương pháp bố trí thí nghiệm đồng ruộng chuẩn hóa theo khối ngẫu nhiên đầy đủ (Randomized Complete Block Design - RCB), kết hợp quan trắc không phá hủy bằng thiết bị quang phổ diệp lục SPAD-502Plus và phương pháp cân sinh khối chuẩn. Kết quả kỳ vọng là xác định chính xác ngưỡng đạm tối ưu (120 kg N/ha) giúp giống lúa TBR125 đạt năng suất thực thu vượt 71 tạ/ha, đảm bảo cấu trúc khóm thông thoáng, hạn chế đổ ngã và giảm thiểu chi phí đầu vào.

Phạm vi nghiên cứu được thực nghiệm trong vụ Xuân 2021 tại vùng sinh thái nông nghiệp huyện Đông Hưng, tỉnh Thái Bình; kết quả làm cơ sở khoa học để mở rộng khuyến nông cho các chân đất vàn, vàn cao vùng Đồng bằng Bắc Bộ.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Thực trạng bón phân đạm cho lúa tại miền Bắc Việt Nam hiện nay tồn tại khoảng cách lớn giữa khuyến cáo kỹ thuật chung và đặc tính sinh học của từng giống lúa mới. Bảng so sánh dưới đây phân tích ưu và nhược điểm của các phương thức quản lý dinh dưỡng đạm hiện hành:

Phương thức quản lý dinh dưỡng Ưu điểm Nhược điểm Rủi ro kỹ thuật
Tập quán nông dân (Bón thừa N > 150 kg/ha) Cây phát triển nhanh giai đoạn đầu, lá xanh đậm Tốn kém chi phí, lúa mềm yếu, đẻ nhiều nhánh vô hiệu Lúa bị lốp đổ khi trỗ, bùng phát rầy nâu và sâu bệnh hại nặng
Quy trình chuẩn chung (Khuyến cáo tĩnh 90 - 100 kg N/ha) Tương đối an toàn, chi phí đầu tư ở mức trung bình Không phát huy hết tiềm năng năng suất của giống chịu thâm canh Năng suất lúa chỉ đạt mức trung bình (55 - 60 tạ/ha)
Quy trình thâm canh định lượng theo giống (Nghiên cứu TBR125) Cân đối dinh dưỡng, tối ưu hóa $AE_N$, năng suất vượt trội $\ge 71$ tạ/ha Đòi hỏi kỹ thuật bón đúng thời điểm và chia đợt nghiêm ngặt Cần tập huấn và chuyển giao công nghệ cho nông dân

Độ ưu tiên yêu cầu kỹ thuật của đề tài được phân loại theo mô hình MoSCoW:

  • Must have (Bắt buộc phải có): Bố trí thí nghiệm RCB 3 lần nhắc lại; xác định rõ 4 mức đạm (60, 90, 120, 150 kg N/ha); bón cân đối nền $60\text{ kg } P_2O_5 + 60\text{ kg } K_2O/\text{ha}$; theo dõi chi tiết 4 yếu tố cấu thành năng suất.
  • Should have (Nên có): Đo quang phổ diệp lục bằng máy SPAD-502Plus; xác định chỉ số diện tích lá LAI theo phương pháp trọng lượng mẫu chuẩn 1 $dm^2$; phân tích động thái tích lũy chất khô qua sấy đẳng trọng ở 80°C.
  • Could have (Có thể mở rộng): Tính toán hiệu quả kinh tế và biên lợi nhuận ròng giữa các công thức đạm.
  • Won't have (Chưa thực hiện đợt này): Đánh giá dư lượng nitrat trong nước ngầm và đo phát thải khí nhà kính tại ruộng.
+-----------------------------------------------------------------------------------+
|                        SƠ ĐỒ BỐ TRÍ THÍ NGHIỆM ĐỒNG RUỘNG                         |
|                             (Kiểu RCB - 3 Lần nhắc lại)                           |
+-----------------------------------------------------------------------------------+
|  Dải bảo vệ xung quanh (Border rows)                                              |
|  +-------------------+  +-------------------+  +-------------------+  +---------+ |
|  | Khối I (Rep 1):   |  | Khối II (Rep 2):  |  | Khối III (Rep 3): |  | Bờ bao  | |
|  | [N1] [N2] [N3] [N4|  | [N2] [N3] [N1] [N4|  | [N3] [N4] [N2] [N1|  | ngăn    | |
|  |                   |  |                   |  |                   |  | nước    | |
|  | Ô = 15m² (3m x 5m)|  | Ô = 15m² (3m x 5m)|  | Ô = 15m² (3m x 5m)|  | 10 cm   | |
|  +-------------------+  +-------------------+  +-------------------+  +---------+ |
|  Kênh cấp / thoát nước chủ động                                                    |
+-----------------------------------------------------------------------------------+

Thiết kế hệ thống

Hệ thống nghiên cứu thực nghiệm đồng ruộng và phân tích số liệu sinh học được thiết lập theo các thông số kỹ thuật chuẩn:

  • Vật liệu cây trồng: Giống lúa thuần TBR125 cảm ôn, ngắn ngày, thời gian sinh trưởng vụ Xuân 125 – 135 ngày, mật độ cấy 35 khóm/$m^2$, cấy 2 dảnh/khóm.
  • Vật liệu phân bón: Phân đạm Urê tinh khiết (46.67% N nguyên chất), Super lân ($16% P_2O_5$), Kali Clorua ($60% K_2O$).
  • Thiết bị đo đạc: Thiết bị đo độ xanh lá lúa Konica Minolta SPAD-502Plus (độ chính xác $\pm 1.0$ SPAD unit), cân phân tích điện tử độ nhạy 0.01g, tủ sấy mẫu thực vật Memmert UN110 (Đức).
  • Phần mềm xử lý thống kê sinh học: Microsoft Excel 365, IRISTAT version 4.0 và ngôn ngữ lập trình R version 4.1.2 phục vụ phân tích phương sai (ANOVA) và kiểm định sai khác trung bình (LSD test ở mức ý nghĩa $\alpha = 0.05$).

Cấu trúc lưu trữ dữ liệu thực nghiệm được chuẩn hóa dưới dạng bảng thuộc tính nhằm đảm bảo tính toàn vẹn dữ liệu:

Table: FieldExperimentRecords
--------------------------------------------------------------------------------
Field_Name          Data_Type       Constraint      Description
--------------------------------------------------------------------------------
Plot_ID             VARCHAR(10)     PRIMARY KEY     Mã định danh ô thí nghiệm
Block_Rep           INT             NOT NULL        Lần nhắc lại (1, 2, 3)
Treatment_N         INT             NOT NULL        Mức đạm bón (60, 90, 120, 150)
Date_Sampled        DATE            NOT NULL        Ngày thu thập số liệu
Weeks_After_Trans   INT             NOT NULL        Số tuần sau cấy (2 đến 9 TSC)
Plant_Height_cm     DECIMAL(5,2)    CHECK(>0)       Chiều cao thân cây (cm)
Tillers_Per_Hill    DECIMAL(4,2)    CHECK(>=1)      Số nhánh trên một khóm
Leaf_Count_Main     DECIMAL(4,2)    CHECK(>=1)      Số lá trên thân chính
SPAD_Value          DECIMAL(4,2)    NULL            Chỉ số SPAD diệp lục lá đòng
LAI_Value           DECIMAL(4,2)    NULL            Chỉ số diện tích lá (m²/m²)
Dry_Matter_g_m2     DECIMAL(6,2)    NULL            Chất khô tích lũy (g/m²)
Panicles_Per_m2     DECIMAL(5,2)    NULL            Số bông trên 1 m² (A)
Spikelets_Per_Pan   DECIMAL(5,2)    NULL            Tổng số hạt trên bông (B)
Filled_Grain_Ratio  DECIMAL(4,2)    NULL            Tỷ lệ hạt chắc (%) (C)
Weight_1000_Grain   DECIMAL(4,2)    NULL            Khối lượng 1000 hạt (g) (D)
Actual_Yield_taha   DECIMAL(4,2)    NULL            Năng suất thực thu (tạ/ha)
--------------------------------------------------------------------------------

Methodology

Phương pháp nghiên cứu tuân thủ nghiêm ngặt quy trình chuẩn của Viện Nghiên cứu Lúa Quốc tế (IRRI):

  • Bố trí công thức: Thí nghiệm đơn nhân tố gồm 4 mức đạm:
    • $N_1$: 60 kg N/ha (Công thức kiểm chứng mức thấp)
    • $N_2$: 90 kg N/ha (Mức trung bình thấp)
    • $N_3$: 120 kg N/ha (Mức thâm canh đề xuất)
    • $N_4$: 150 kg N/ha (Mức thâm canh cao)
  • Kỹ thuật bón phân: Áp dụng phương pháp bón đạm phân kỳ nhằm tránh thất thoát dinh dưỡng:
    • Bón lót trước khi cấy: $100% P_2O_5 + 50% N$.
    • Bón thúc lần 1 (10 ngày sau cấy - hồi xanh đẻ nhánh): $30% N + 50% K_2O$.
    • Bón thúc lần 2 (20 ngày sau cấy - đón đòng): $20% N + 50% K_2O$.
  • Lấy mẫu quan trắc: Mỗi ô cố định 5 khóm đại diện để theo dõi động thái sinh trưởng hàng tuần (7 ngày/lần). Đối với chỉ tiêu phá hủy (LAI và chất khô), lấy ngẫu nhiên 3 khóm/ô tại 3 mốc: đẻ nhánh rộ, trỗ 80% và chín sáp.

Implementation và kết quả

Development process

Quá trình triển khai thực nghiệm đồng ruộng và phân tích nông học được phân chia thành các giai đoạn rõ ràng:

  1. Giai đoạn chuẩn bị mạ và làm đất: Gieo mạ ngày 29/01/2021, chăm sóc mạ xuân đủ 19 ngày tuổi, đạt 4.0 lá, cây đanh dảnh, sạch bệnh trước khi cấy vào ngày 17/02/2021.
  2. Giai đoạn điều tiết dinh dưỡng và nước: Phân bổ đạm theo đúng tiến trình phân kỳ, giữ mực nước chân ruộng 3 – 5 cm giai đoạn đẻ nhánh, phơi ruộng rút nước nhẹ giai đoạn cuối đẻ nhánh để hạn chế nhánh vô hiệu, sau đó cấp nước nông khi lúa làm đòng và trỗ bông.
  3. Giai đoạn thu thập và tính toán thuật toán nông học:

Công thức tính chỉ số diện tích lá (LAI): $$LAI = \frac{P_2 \times M}{P_1 \times 100}$$ Trong đó: $P_1$ là khối lượng tươi của diện tích lá mẫu $1\text{ dm}^2$ ($0.01\text{ m}^2$), $P_2$ là tổng khối lượng lá tươi của khóm lúa ($g$), $M$ là mật độ cấy ($35\text{ khóm/m}^2$).

Công thức tính năng suất lý thuyết (NSLT): $$NSLT (\text{tạ/ha}) = \frac{A \times B \times C \times D}{10.000}$$ Trong đó: $A$ là số bông/$m^2$; $B$ là số hạt/bông; $C$ là tỷ lệ hạt chắc ($0 < C \le 1$); $D$ là khối lượng 1000 hạt ($g$).

Đoạn mã phân tích thống kê R chuẩn hóa việc kiểm định phương sai ANOVA và so sánh giá trị trung bình bằng thuật toán LSD:

# R Script: Phân tích phương sai ANOVA và kiểm định LSD cho thí nghiệm bón đạm
library(agricolae)

# 1. Nhập dữ liệu năng suất thực thu của 4 công thức đạm TBR125
data_yield <- data.frame(
  Block = factor(rep(1:3, each = 4)),
  Treatment = factor(rep(c("N1_60N", "N2_90N", "N3_120N", "N4_150N"), times = 3)),
  Yield_taha = c(60.8, 66.5, 71.0, 69.2,   # Block 1
                 61.5, 67.0, 71.3, 69.8,   # Block 2
                 61.3, 66.9, 71.0, 69.2)   # Block 3
)

# 2. Xây dựng mô hình tuyến tính phân tích phương sai khối ngẫu nhiên (RCB)
model <- lm(Yield_taha ~ Block + Treatment, data = data_yield)
anova_table <- anova(model)
print(anova_table)

# 3. Phân tích so sánh bội Least Significant Difference (LSD test at alpha = 0.05)
lsd_result <- LSD.test(model, "Treatment", p.adj = "none", alpha = 0.05, group = TRUE)
print(lsd_result$groups)

Testing và validation

Số liệu thí nghiệm được kiểm định qua hệ số biến thiên ($CV%$) và khoảng sai khác có ý nghĩa nhỏ nhất ($LSD_{0.05}$). Các phép thử đều cho thấy độ tin cậy thực nghiệm cao:

  • Hệ số biến thiên chiều cao cây ($CV% = 1.6%$), số lá ($CV% = 0.0%$), số nhánh đẻ ($CV% = 1.7% - 2.4%$), phản ánh điều kiện đồng ruộng đồng đều và sai số ngoại cảnh được kiểm soát nghiêm ngặt.
  • Hệ số biến thiên chất khô tích lũy ($CV% = 5.2% - 7.9%$) và LAI ($CV% = 2.3% - 9.2%$) nằm trong ngưỡng cho phép đối với các thí nghiệm sinh học đồng ruộng ($CV < 10%$).

Kết quả đạt được

Kết quả tổng hợp ảnh hưởng của lượng đạm bón đến các chỉ tiêu sinh trưởng, sinh lý và năng suất của giống lúa TBR125 được thể hiện chi tiết tại bảng số liệu thực nghiệm dưới đây:

Chỉ tiêu theo dõi Mức $N_1$ (60 kg N) Mức $N_2$ (90 kg N) Mức $N_3$ (120 kg N) Mức $N_4$ (150 kg N) $CV (%)$ $LSD_{0.05}$
Chiều cao cuối cùng (cm) $113.3^c$ $114.7^{bc}$ $117.3^b$ $121.6^a$ 1.6 3.83
Số lá trên thân chính (lá) $14.0^a$ $14.0^a$ $14.0^a$ $14.0^a$ 0.0 0.00
Số nhánh hữu hiệu (nhánh/khóm) $5.9^b$ $6.0^{ab}$ $6.2^{ab}$ $6.3^a$ 2.4 0.30
LAI trỗ 80% ($m^2\text{ lá}/m^2\text{ đất}$) $4.3^b$ $5.7^a$ $5.9^a$ $6.1^a$ 9.2 1.02
Chất khô chín sáp ($g/m^2$) $826.9^c$ $951.8^b$ $1023.1^a$ $1024.4^a$ 6.6 25.78
Số bông hữu hiệu/m² ($A$) 206.5 210.0 217.0 220.5 - -
Tổng số hạt/bông ($B$) 158.2 165.4 175.5 171.2 - -
Tỷ lệ hạt chắc ($%$) ($C$) $86.1%$ $87.4%$ $88.3%$ $85.0%$ - -
Khối lượng 1000 hạt ($g$) ($D$) 23.2 23.4 23.5 23.1 - -
Năng suất thực thu (tạ/ha) $61.2^c$ $66.8^b$ $71.1^a$ $69.4^a$ 2.1 2.45

Ghi chú: Các chữ cái khác nhau trong cùng một hàng thể hiện sự sai khác có ý nghĩa thống kê ở mức tin cậy 95% ($p < 0.05$).

Biểu đồ so sánh Năng suất thực thu và Tỷ lệ hạt chắc theo 4 mức đạm
75 +-------------------------------------------------------+ 90%
   |                                    [71.1]             |
70 |                     [66.8]          (88.3)    [69.4]  | 88%
   |      [61.2]          (87.4)                   (85.0)  |
65 |       (86.1)                                          | 86%
   |                                                       |
60 +-------------------------------------------------------+ 84%
       N1 (60kg N)     N2 (90kg N)     N3 (120kg N)   N4 (150kg N)
       [--] NSTT (tạ/ha)               (--) Tỷ lệ hạt chắc (%)

Phân tích chuyên sâu cho thấy:

  • Động thái tăng trưởng chiều cao và đẻ nhánh: Chiều cao cây tăng nhanh nhất ở tuần thứ 4 sau cấy (đạt mức tăng 24.7 – 30.0 cm/tuần). Mức bón $N_4$ cho chiều cao lớn nhất (121.6 cm), làm tăng nguy cơ ngã đổ cuối vụ. Số nhánh đạt cực đại vào tuần thứ 6 (10.2 – 13.6 nhánh/khóm) sau đó giảm dần do sự tiêu biến của nhánh vô hiệu.
  • Động thái sinh lý: Chỉ số diện tích lá (LAI) đạt đỉnh tại thời kỳ trỗ 80%, trong đó mức bón $N_3$ ($5.9\text{ m}^2/\text{m}^2$) và $N_4$ ($6.1\text{ m}^2/\text{m}^2$) sai khác không có ý nghĩa thống kê nhưng vượt trội so với $N_1$ ($4.3\text{ m}^2/\text{m}^2$). Chất khô tích lũy tăng dần và đạt tối đa ở giai đoạn chín sáp; mức $N_3$ đạt $1023.1\text{ g/m}^2$, tương đương mức $N_4$ ($1024.4\text{ g/m}^2$).
  • Năng suất và cấu thành năng suất: Mức bón $N_3$ (120 kg N/ha) đạt năng suất thực thu cao nhất (71.1 tạ/ha), cao hơn $N_1$ là $16.18%$ và cao hơn $N_2$ là $6.44%$. Mặc dù mức bón $N_4$ (150 kg N/ha) có số bông/$m^2$ cao hơn (220.5 bông) nhưng tỷ lệ hạt chắc giảm mạnh xuống $85.0%$ (so với $88.3%$ ở $N_3$), dẫn đến năng suất thực thu bị suy giảm ($69.4$ tạ/ha) mà không mang lại hiệu quả gia tăng.

Đổi mới và đóng góp

Nghiên cứu mang lại những đóng góp khoa học và kỹ thuật quan trọng trong canh tác lúa thuần năng suất cao:

  1. Xác lập chuẩn định lượng phân đạm cho giống TBR125: Chứng minh rõ ràng quy luật bão hòa đạm nông học; nâng mức đạm từ 120 kg N lên 150 kg N/ha không làm tăng năng suất có ý nghĩa mà làm tăng chi phí và rủi ro sâu bệnh hại.
  2. So sánh đối chứng với các nghiên cứu tiền nhiệm:
    • So với nghiên cứu của Nguyễn Quốc Khương và cs. (2017) trên đất phù sa Vĩnh Long (hiệu quả nông học tối ưu đạt tại mức 120 kg N/ha cho năng suất 7.05 tấn/ha), kết quả trên giống TBR125 tại Thái Bình đạt 7.11 tấn/ha, khẳng định mức bón 120 kg N/ha có tính ổn định sinh thái cao.
    • So với nghiên cứu của Jiuxin Guo và cs. (2017) tại Giang Tô - Trung Quốc (cắt giảm đạm từ mức nông dân 336.7 kg N/ha xuống 240 kg N/ha giúp tăng 5.9% năng suất), nghiên cứu này chứng minh việc tiết giảm 30 kg N/ha (từ $N_4$ xuống $N_3$) giúp giữ vững năng suất đỉnh cao mà tiết kiệm 20% lượng phân bón hóa học.
    • Kết quả cũng tương thích với công bố của Gewaily và cs. (2018) tại Ai Cập khi chỉ ra sự gia tăng liều đạm quá mức ngưỡng làm giảm tỷ lệ quang hợp thuần của bộ lá đòng và làm tăng tỷ lệ hạt lép lửng.
  3. Hiệu quả cải tiến: Tiết kiệm chính xác 30 kg đạm nguyên chất/ha (tương đương 64.3 kg phân Urê/ha), hạn chế tỷ lệ bệnh bạc lá và bệnh khô vằn phát sinh trong vụ Xuân ẩm ướt.

Ứng dụng thực tế và triển khai

Kịch bản triển khai sản xuất đại trà

Quy trình bón phân tối ưu được chuyển giao cho các hộ nông dân và Hợp tác xã Nông nghiệp tại Thái Bình theo lộ trình cụ thể:

+-------------------------------------------------------------------------------+
|                       LỘ TRÌNH TRIỂN KHAI VÀ BÓN PHÂN                         |
+-------------------------------------------------------------------------------+
| Giai đoạn 1: Làm đất & Bón lót (0 ngày sau cấy - 17/02/2021)                  |
| - Bón 100% Lân Super (375 kg/ha) + 50% Đạm Urê (128.5 kg/ha)                  |
|                                                                               |
| Giai đoạn 2: Bón thúc lần 1 - Đẻ nhánh (10 ngày sau cấy - 27/02/2021)         |
| - Bón 30% Đạm Urê (77.1 kg/ha) + 50% Kali Clorua (50 kg/ha)                   |
| - Giữ nước 3-4 cm, làm sạch cỏ dại hoàn toàn                                 |
|                                                                               |
| Giai đoạn 3: Bón thúc lần 2 - Đón đòng (20 ngày sau cấy - 09/03/2021)         |
| - Bón 20% Đạm Urê (51.4 kg/ha) + 50% Kali Clorua (50 kg/ha)                   |
| - Điều tiết nước khô ướt xen kẽ, kiểm tra sâu cuốn lá                         |
+-------------------------------------------------------------------------------+

Phân tích chi phí - lợi ích (Cost-Benefit Analysis)

Tính toán hiệu quả kinh tế trên quy mô 1 hecta (giá thóc thương phẩm tạm tính 8.000 VNĐ/kg, giá Urê 15.000 VNĐ/kg):

  • Công thức đề xuất $N_3$ (120 kg N/ha):
    • Chi phí phân Urê: 257.1 kg $\times$ 15.000 = 3.856.500 VNĐ.
    • Tổng thu từ thóc: 7.110 kg $\times$ 8.000 = 56.880.000 VNĐ.
    • Doanh thu trừ chi phí đạm: 53.023.500 VNĐ/ha.
  • Công thức bón thừa $N_4$ (150 kg N/ha):
    • Chi phí phân Urê: 321.4 kg $\times$ 15.000 = 4.821.000 VNĐ.
    • Tổng thu từ thóc: 6.940 kg $\times$ 8.000 = 55.520.000 VNĐ.
    • Doanh thu trừ chi phí đạm: 50.699.000 VNĐ/ha.
  • Hiệu quả chênh lệch: Áp dụng mức bón $N_3$ giúp tăng thêm lợi nhuận ròng 2.324.500 VNĐ/ha so với mức bón $N_4$ nhờ tiết kiệm phân bón và đạt năng suất hạt chắc cao hơn.

Hạn chế và hướng phát triển

Dự án ghi nhận một số giới hạn kỹ thuật và định hướng nghiên cứu tiếp theo:

  1. Hạn chế thực nghiệm: Thí nghiệm mới được tiến hành trong 1 vụ Xuân 2021 trên một loại đất phù sa tại Đông Hưng, Thái Bình; chưa khảo nghiệm so sánh đồng thời trong vụ Mùa (vốn có điều kiện bức xạ nhiệt và áp lực sâu bệnh khác biệt).
  2. Ảnh hưởng ngoại cảnh: Điều kiện thời tiết vụ Xuân 2021 có độ ẩm cao, mưa phùn kéo dài giai đoạn đầu làm chậm tốc độ thoát ẩm của đất khi rút nước phơi ruộng.
  3. Hướng phát triển:
    • Nghiên cứu kết hợp bón đạm với các dạng phân bón tan chậm (Slow-release fertilizers) hoặc phân bón hữu cơ vi sinh nhằm nâng cao hơn nữa hiệu suất hấp thu đạm của cây lúa vượt ngưỡng 50%.
    • Ứng dụng bảng so màu lá lúa (Leaf Color Chart - LCC) số hóa hoặc cảm biến viễn thám NDVI không người lái (UAV) để điều chỉnh lượng đạm bón thời gian thực (Real-time Nitrogen Management).

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

  • Sinh viên và Học viên Nông học: Tài liệu mẫu chuẩn mực về phương pháp nghiên cứu cây trồng, bố trí thí nghiệm đồng ruộng RCB và quy trình đo đạc sinh lý (LAI, SPAD, chất khô).
  • Kỹ sư nông nghiệp và Cán bộ Khuyến nông: Cung cấp thông số định lượng chính xác để ban hành sổ tay hướng dẫn kỹ thuật canh tác giống TBR125 cho các vùng thâm canh lúa miền Bắc.
  • Hợp tác xã Nông nghiệp và Nông dân: Nắm bắt quy tắc bón phân đúng liều lượng và thời điểm, giúp tăng thu nhập trực tiếp từ 2,3 – 3,5 triệu VNĐ/ha/vụ và giảm thiểu rủi ro đổ ngã lúa.
  • Nhà khoa học và Chọn tạo giống: Bổ sung dữ liệu sinh lý và dinh dưỡng của giống TBR125 vào ngân hàng dữ liệu các giống lúa thuần chất lượng cao của Việt Nam.

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

  1. Vì sao mức đạm $N_4$ (150 kg N/ha) có sinh khối chất khô cao nhưng năng suất thực thu lại thấp hơn mức $N_3$ (120 kg N/ha)? Trả lời: Do bón thừa đạm ở mức $N_4$ kích thích thân lá phát triển quá mạnh, tạo nhiều nhánh vô hiệu, tán lá rậm rạp che khuất quang hợp tầng dưới. Điều này làm cản trở quá trình vận chuyển hydratcarbon không cấu trúc từ thân lá về bông giai đoạn làm chín, khiến tỷ lệ hạt chắc giảm từ $88.3%$ xuống còn $85.0%$ và hạt lép tăng lên.

  2. Chỉ số SPAD đo bằng máy SPAD-502Plus có ý nghĩa gì trong việc điều tiết phân đạm? Trả lời: Chỉ số SPAD phản ánh mật độ diệp lục tố trên bề mặt lá, tương quan thuận trực tiếp với hàm lượng đạm trong cây. Ở giai đoạn trỗ 80%, giá trị SPAD của mức bón $N_3$ đạt trạng thái lý tưởng giúp duy trì khả năng quang hợp bền bỉ của lá đòng mà không làm lá bị mềm yếu hay nhiễm bệnh bạc lá.

  3. Có thể áp dụng nguyên vẹn mức bón 120 kg N/ha này cho vụ Mùa được không? Trả lời: Không nên áp dụng nguyên vẹn. Vụ Mùa có nền nhiệt độ và ánh sáng cao, kèm theo mưa nhiều, cây lúa sinh trưởng nhanh hơn và đất giải phóng đạm khoáng hóa nhiều hơn. Lượng đạm khuyến cáo cho vụ Mùa nên giảm khoảng $15 - 20%$ (xuống mức 95 – 100 kg N/ha) để chống hiện tượng lốp đổ non.

  4. Tỷ lệ phân chia đạm bón lót và bón thúc trong nghiên cứu có gì đặc biệt? Trả lời: Tỷ lệ $50%$ bón lót + $30%$ bón thúc đẻ nhánh + $20%$ bón thúc đón đòng giúp cung cấp đạm đúng theo nhu cầu sinh lý của lúa: đạm bón lót giúp bén rễ nhanh, thúc 1 giúp tăng số nhánh hữu hiệu sớm, và thúc 2 giúp tăng số hạt trên bông và tỷ lệ chắc hạt, giảm thiểu việc rửa trôi đạm trong môi trường nước.

  5. Giống lúa TBR125 có khả năng chống chịu sâu bệnh như thế nào khi bón đạm ở mức tối ưu? Trả lời: Ở mức bón 120 kg N/ha, giống TBR125 thể hiện tính cứng cây, lá đòng đứng thẳng, mức độ nhiễm các bệnh đạo ôn, bạc lá, khô vằn và sâu cuốn lá chỉ ở cấp độ nhẹ (Cấp 1 - Cấp 3 theo thang phân cấp IRRI), hoàn toàn nằm trong tầm kiểm soát sinh học mà không cần phun xịt thuốc bảo vệ thực vật quá mức.

Kết luận

Khóa luận tốt nghiệp đã hoàn thành toàn diện các mục tiêu nghiên cứu thông qua phương pháp thực nghiệm khoa học chuẩn xác. Kết quả khẳng định mức bón đạm 120 kg N/ha (kết hợp nền $60\text{ kg } P_2O_5 + 60\text{ kg } K_2O/\text{ha}$) là định mức dinh dưỡng kỹ thuật tối ưu nhất cho giống lúa tẻ TBR125 trong vụ Xuân tại Thái Bình. Tại mức bón này, cây lúa đạt chiều cao hợp lý (117.3 cm), tỷ lệ nhánh hữu hiệu cao (6.2 nhánh/khóm), chỉ số diện tích lá cân đối ($5.9\text{ m}^2/\text{m}^2$), chất khô tích lũy lớn ($1023.1\text{ g/m}^2$), đưa năng suất thực thu đạt đỉnh 71.1 tạ/ha và tối đa hóa lợi nhuận ròng cho người trồng lúa.

Kết quả nghiên cứu cung cấp cơ sở dữ liệu quan trọng để các cơ quan khuyến nông, viện nghiên cứu và doanh nghiệp giống hoàn thiện quy trình canh tác chuẩn, hướng tới nền nông nghiệp thâm canh lúa gạo bền vững, hiệu quả và thân thiện với môi trường sinh thái.