Giới thiệu dự án

Sản xuất lúa cạn (Oryza sativa L.) giữ vai trò chiến lược trong việc đảm bảo an ninh lương thực và sinh kế bền vững cho đồng bào các dân tộc thiểu số tại các vùng trung du và miền núi Việt Nam. Theo thống kê của Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc (FAO), diện tích lúa canh tác nhờ nước trời và lúa rẫy chiếm tỷ trọng đáng kể nhưng năng suất chỉ đạt từ 0,8 - 1,7 tấn/ha (bằng 20 - 40% so với vùng thâm canh lúa nước chủ động tưới tiêu). Tại Việt Nam, vùng Trung du và Miền núi phía Bắc có khoảng 210.000 ha lúa cạn, phân bố chủ yếu trên địa hình đồi dốc phức tạp.

Giống lúa nếp cạn Khẩu Nua Trạng (thu thập và phục tráng tại Hà Giang giai đoạn 2013 - 2015) là nguồn gen bản địa quý hiếm có phẩm chất gạo thơm ngon, dẻo đặc trưng và khả năng chịu hạn cao. Tuy nhiên, tập quán canh tác gieo hốc trên đất dốc khiến ruộng lúa đối mặt với áp lực cạnh tranh gay gắt từ quần thể cỏ dại. Cứ 100 cây cỏ/$m^2$ làm giảm 17% năng suất lúa, mật độ 100 - 200 cây cỏ/$m^2$ khiến năng suất giảm thêm 10%. Lao động thủ công làm cỏ chiếm tới hơn 50% tổng công lao động (khoảng 150 công/ha), làm tăng chi phí sản xuất và hạn chế khả năng mở rộng quy mô.

Mục tiêu nghiên cứu cụ thể

  1. Đánh giá mức độ ảnh hưởng của 5 phương thức trừ cỏ (kết hợp hóa chất tiền - hậu nảy mầm và làm cỏ thủ công) đến các chỉ tiêu nông sinh học và động thái đẻ nhánh của giống lúa nếp cạn Khẩu Nua Trạng.
  2. Xác định tác động của từng phương thức trừ cỏ đến áp lực phát sinh sâu bệnh hại chính (sâu đục thân, sâu cuốn lá, bệnh đạo ôn, bệnh khô vằn).
  3. Đánh giá tác động đến các yếu tố cấu thành năng suất ($P_{1000}$, số bông/$m^2$, số hạt chắc/bông) và năng suất thực thu ($NSTT$), từ đó đề xuất quy trình kỹ thuật tối ưu hóa chi phí và công lao động.

Phương pháp tiếp cận và phạm vi nghiên cứu

  • Tiếp cận: Ứng dụng quy trình quản lý cỏ dại tích hợp (IWM - Integrated Weed Management), kết hợp xử lý hóa chất tiền nảy mầm (Lyphoxim 40SL) triệt tiêu hạt cỏ đa niên/hàng niên trước gieo với thuốc trừ cỏ chọn lọc hậu nảy mầm (Mizin 80WP) giai đoạn 1 - 3 lá.
  • Phạm vi: Thí nghiệm đồng ruộng bố trí tại Trường Đại học Nông Lâm Thái Nguyên trong vụ Mùa 2016 (từ 10/07/2016 đến 21/11/2016) trên nền đất đồi dốc dưới 5%, đất thoát nước tự nhiên hoàn toàn phụ thuộc nước trời.

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

Phân tích hiện trạng quản lý cỏ dại

Tiêu chí so sánh Làm cỏ thủ công truyền thống Sử dụng thuốc trừ cỏ đơn lẻ Quy trình kết hợp đa pha (Đề xuất)
Công lao động Rất cao (120 - 150 công/ha/vụ) Trung bình (15 - 20 công/ha) Thấp nhất (5 - 8 công/ha)
Phổ diệt cỏ Hạn chế, dễ sót rễ cỏ đa niên Phụ thuộc vào loại thuốc đơn Diệt toàn diện hòa bản, lá rộng, cói lác
Ảnh hưởng sinh trưởng lúa Đứt rễ non khi xới xáo Nguy cơ cháy lá nếu sai liều An toàn sinh học, kích thích đẻ nhánh sớm
Hiệu quả kinh tế Chi phí nhân công quá cao Hiệu quả ngắn hạn Năng suất tăng 14,4%, ROI tối ưu

Phân loại yêu cầu kỹ thuật theo mô hình MoSCoW:

  • Must have: Diệt trừ triệt để nhóm cỏ hòa bản và lá rộng giai đoạn 0 - 30 ngày sau gieo; không gây ngộ độc rễ và thân mầm giống lúa Khẩu Nua Trạng.
  • Should have: Giảm thiểu tối thiểu 70% công làm cỏ bằng tay; duy trì tỷ lệ trỗ thoát cổ bông cấp 1.
  • Could have: Hạn chế môi trường cư trú của sâu đục thân (Scirpophaga incertulas) và rầy mềm.
  • Won't have: Sử dụng các hóa chất BVTV tồn lưu lâu trong đất thuộc danh mục cấm hoặc thuốc trừ cỏ không chọn lọc sau khi lúa đã mọc.

Thiết kế hệ thống thí nghiệm

Thí nghiệm được thiết kế theo khối ngẫu nhiên hoàn chỉnh (RCBD - Randomized Complete Block Design) với 5 công thức và 3 lần nhắc lại ($n = 15$ ô thí nghiệm, diện tích mỗi ô là $30\text{ m}^2$ kích thước $5\text{ m} \times 6\text{ m}$, dải cách ly $1\text{ m}$, xung quanh có dải bảo vệ $1\text{ m}$).

[Dải bảo vệ 1m xung quanh]

Thông số kỹ thuật của các công thức:

  1. CT1 (Đối chứng): Làm cỏ thủ công bằng tay 3 lần (25, 45 và 60 ngày sau gieo - NSG).
  2. CT2: Phun Lyphoxim 40SL trước gieo 7 ngày + Làm cỏ tay 1 lần (25 NSG).
  3. CT3: Phun Lyphoxim 40SL trước gieo 7 ngày + Làm cỏ tay 2 lần (25 và 45 NSG).
  4. CT4: Phun Lyphoxim 40SL trước gieo 7 ngày + Phun Mizin 80WP (25 NSG, khi cỏ 1 - 3 lá).
  5. CT5: Phun Lyphoxim 40SL trước gieo 7 ngày + Phun Mizin 80WP (25 NSG) + Làm cỏ tay 1 lần (45 NSG).

Hoạt chất và thông số nông hóa

  • Lyphoxim 40SL: Hoạt chất Glyphosate Isopropylamine salt 41%. Cơ chế tác động: Nội hấp qua thân lá còn xanh, lưu dẫn xuống hệ rễ và thân ngầm làm cỏ chết hoàn toàn sau 4 - 7 ngày. Định mức: $4\text{ lít/ha}$ pha $500\text{ lít nước}$ (hoặc $70 - 85\text{ ml}$ bình 8 lít, phun $6\text{ bình}/1000\text{ m}^2$).
  • Mizin 80WP: Hoạt chất Atrazine 80% + phụ gia 20%. Cơ chế tác động: Nội hấp chọn lọc qua rễ, thân mầm và lá; tác động tiền nảy mầm và hậu nảy mầm sớm, ức chế quang hợp của cỏ hòa bản và lá rộng hàng niên. Định mức: $500\text{ g/ha}$ pha $400\text{ lít nước}$.
  • Quy chuẩn khảo nghiệm: Tuân thủ tiêu chuẩn ngành 10 TCN 558-2011 (Quy phạm khảo nghiệm giá trị canh tác và sử dụng của giống lúa).

Implementation và kết quả

Quy trình kỹ thuật và thuật toán xử lý dữ liệu

Quy trình thâm canh áp dụng thống nhất:

  • Mật độ gieo: Hốc cách hốc $25 - 30\text{ cm}$, hàng cách hàng $20\text{ cm}$, gieo $5 - 7\text{ hạt/hốc}$, lấp đất sâu $3 - 5\text{ cm}$.
  • Bón phân nền cho 1 ha: 1 tấn phân vi sinh + 300 kg vôi bột (bón lót toàn bộ). Bón thúc 1 (15 - 20 NSG): 60% Đạm Urê + 40% Kali Clorua. Bón thúc 2 (50 - 60 NSG): 40% Đạm Urê + 60% Kali Clorua.

Thuật toán phân tích phương sai (ANOVA RCBD) và tính toán các chỉ tiêu cấu thành năng suất được triển khai bằng phần mềm chuyên dụng IRRISTAT 5.0 và mô hình hóa kiểm định qua script phân tích:

import numpy as np
import scipy.stats as stats

def anova_rcbd_analysis(treatments, blocks, yield_data):
    """
    Phân tích ANOVA cho thí nghiệm khối ngẫu nhiên hoàn chỉnh (RCBD)
    treatments: Số công thức (t = 5)
    blocks: Số lần nhắc lại / Khối (r = 3)
    yield_data: Mảng 2 chiều kích thước (t, r)
    """
    t = len(treatments)
    r = blocks
    n_total = t * r
    
    grand_total = np.sum(yield_data)
    cf = (grand_total ** 2) / n_total
    
    ss_total = np.sum(yield_data ** 2) - cf
    ss_treatments = np.sum(np.sum(yield_data, axis=1) ** 2) / r - cf
    ss_blocks = np.sum(np.sum(yield_data, axis=0) ** 2) / t - cf
    ss_error = ss_total - ss_treatments - ss_blocks
    
    df_treatments = t - 1
    df_blocks = r - 1
    df_error = (t - 1) * (r - 1)
    
    ms_treatments = ss_treatments / df_treatments
    ms_error = ss_error / df_error
    
    f_val = ms_treatments / ms_error
    p_val = 1 - stats.f.cdf(f_val, df_treatments, df_error)
    
    # Tính sai khác nhỏ nhất có ý nghĩa (LSD 0.05)
    t_critical = stats.t.ppf(0.975, df_error)
    lsd_05 = t_critical * np.sqrt(2 * ms_error / r)
    cv_percent = (np.sqrt(ms_error) / np.mean(yield_data)) * 100
    
    return {
        "F_value": f_val,
        "P_value": p_val,
        "LSD_0.05": lsd_05,
        "CV_percent": cv_percent
    }

Công thức tính năng suất lý thuyết: $$NSLT \text{ (tạ/ha)} = \frac{\text{Số bông}/m^2 \times \text{Số hạt chắc/bông} \times P_{1000} \text{ (g)}}{10.000}$$

Kết quả nông sinh học và động thái đẻ nhánh

Các công thức xử lý không làm biến đổi chu kỳ sinh trưởng của giống lúa: Thời gian mọc (7 NSG), bắt đầu đẻ nhánh (25 NSG), đẻ nhánh tối đa (55 NSG), trỗ $10%$ (101 NSG), trỗ hoàn toàn $80%$ (106 NSG), và chín hoàn toàn (134 NSG). Chiều cao cây trung bình dao động từ 130,17 - 131,38 cm ($p > 0,05$), độ thoát cổ bông đạt cấp 1 (thoát hoàn toàn) và độ rụng hạt đạt cấp 1 ($< 10%$ hạt rụng khi vuốt).

Công thức Dảnh ban đầu (dảnh/khóm) Nhánh tối đa (nhánh/khóm) Nhánh hữu hiệu (nhánh/khóm) Tỷ lệ nhánh hữu hiệu (%)
CT1 (ĐC) 3 10,17 5,20 50,98
CT2 3 11,10 5,43 47,22
CT3 3 11,10 5,27 47,78
CT4 3 12,13 6,20 50,94
CT5 3 11,03 5,50 46,89
$P$-value - $< 0,05$ ($p = 0,010$) $< 0,05$ ($p = 0,031$) -
$LSD_{0,05}$ - 0,8 0,6 -
CV (%) - 4,1 5,8 -

CT4 (Lyphoxim + Mizin) đạt số nhánh tối đa cao nhất (12,13 nhánh/khóm) và số nhánh hữu hiệu vượt trội (6,20 nhánh/khóm, cao hơn có ý nghĩa thống kê so với đối chứng CT1 ở mức tin cậy 95%).

Mức độ phát sinh sâu bệnh hại

Thang điểm sâu hại (Đẻ nhánh / Chín):
CT1 (Làm cỏ tay)    : [Sâu đục thân: Điểm 5] [Sâu cuốn lá: Điểm 1] [Đạo ôn: Điểm 1]
CT2 (Lyphoxim + 1x) : [Sâu đục thân: Điểm 5] [Sâu cuốn lá: Điểm 1] [Đạo ôn: Điểm 1]
CT3 (Lyphoxim + 2x) : [Sâu đục thân: Điểm 3] [Sâu cuốn lá: Điểm 1] [Đạo ôn: Điểm 1]
CT4 (Lyphoxim+Mizin): [Sâu đục thân: Điểm 3] [Sâu cuốn lá: Điểm 1] [Đạo ôn: Điểm 1]
CT5 (Lyph+Miz+1x)   : [Sâu đục thân: Điểm 3] [Sâu cuốn lá: Điểm 1] [Đạo ôn: Điểm 1]
  • Sâu đục thân: Gây hại giai đoạn đẻ nhánh ở CT1 và CT2 với tỷ lệ 21 - 30% dảnh héo/bông bạc (điểm 5). Ở CT4 và CT5, ruộng thông thoáng hơn, mật độ cỏ giảm làm giảm nơi trú ngụ của sâu, mức độ hại giảm xuống điểm 3 (11 - 20% dảnh héo).
  • Bệnh hại: Bệnh đạo ôn hại lá (Pyricularia oryzae) và bệnh khô vằn (Rhizoctonia solani) đều ở mức nhẹ (điểm 1) trên toàn bộ 5 công thức thí nghiệm.

Năng suất và các yếu tố cấu thành năng suất

Công thức Số bông/$m^2$ Tổng số hạt/bông Số hạt chắc/bông Khối lượng $P_{1000}$ (g) $NSLT$ (tạ/ha) $NSTT$ (tạ/ha) Chênh lệch so với ĐC (%)
CT1 (ĐC) 147 146,74 32,27 29,07 22,77 27,33 -
CT2 158 158,32 35,33 31,80 21,53 26,80 -1,94%
CT3 160 160,18 51,20 29,77 22,41 26,77 -2,05%
CT4 167 167,33 52,40 33,17 34,59 31,27 +14,42%
CT5 180 179,53 49,13 31,50 32,60 30,07 +10,03%
$P$-value $> 0,05$ $> 0,05$ $< 0,05$ $< 0,05$ ($p = 0,044$) $< 0,05$ ($p = 0,017$) $< 0,05$ ($p = 0,032$) -
$LSD_{0,05}$ - - 3,3 8,92 - 3,1 -
CV (%) 8,5 23,3 5,7 18,0 - 5,9 -

Phương thức trừ cỏ CT4 cho năng suất thực thu cao nhất ($31,27\text{ tạ/ha} = 3,127\text{ tấn/ha}$), vượt đối chứng làm cỏ thủ công $3,94\text{ tạ/ha}$ (+14,42%, với $p = 0,032 < 0,05$ và $LSD_{0,05} = 3,1\text{ tạ/ha}$). Khối lượng 1000 hạt của CT4 cũng đạt giá trị cực đại ($33,17\text{ g}$).


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

  1. Quy trình diệt cỏ kép "Tiền nảy mầm – Hậu nảy mầm sớm": Khắc phục triệt để nhược điểm của biện pháp phun thuốc đơn pha. Lyphoxim 40SL làm sạch hoàn toàn cỏ dại lưu niên (cỏ tranh, cỏ ống) trước khi tra hạt; Mizin 80WP tạo màng chắn hóa học ngăn chặn lứa hạt cỏ thứ cấp nảy mầm trong 30 ngày đầu mà không ức chế sinh trưởng của mầm lúa nếp cạn.
  2. Cắt đứt chuỗi ký chủ trung gian của sâu hại: Giảm thiểu 40% tỷ lệ dảnh héo do sâu đục thân so với tập quán làm cỏ bằng tay muộn (hạ cấp gây hại từ điểm 5 xuống điểm 3).
  3. Tối ưu hóa năng suất và hiệu suất lao động: Gia tăng 14,42% năng suất thực thu trên diện tích đất dốc nhờ bảo toàn dinh dưỡng và mật độ bông hữu hiệu ($6,20\text{ cành/khóm}$ so với $5,20\text{ cành/khóm}$ ở đối chứng).
  4. Đóng góp vào cơ sở dữ liệu canh tác lúa cạn: Bổ sung dẫn liệu khoa học chính xác cho bộ giống lúa đặc sản phục tráng Khẩu Nua Trạng tại vùng Trung du và Miền núi phía Bắc.

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

Kịch bản ứng dụng thực tế

Quy trình phù hợp cho các vùng canh tác lúa nương rẫy định canh, ruộng bậc thang và đất đồi dốc dưới 15% tại Hà Giang, Thái Nguyên, Cao Bằng, Bắc Kạn, Yên Bái và Sơn La.

                    LỘ TRÌNH TRIỂN KHAI VỤ MÙA

Phân tích chi phí - lợi ích (Tính trên 1 ha canh tác)

Khoản mục Làm cỏ thủ công (CT1) Quy trình CT4 (Lyphoxim + Mizin) Hiệu quả chênh lệch
Chi phí thuốc trừ cỏ 0 VNĐ 450.000 VNĐ +450.000 VNĐ
Công phun thuốc 0 công 4 công (1.000.000 VNĐ) +1.000.000 VNĐ
Công làm cỏ thủ công 120 công (30.000.000 VNĐ) 0 công -30.000.000 VNĐ
Tổng chi phí trừ cỏ 30.000.000 VNĐ 1.450.000 VNĐ Tiết kiệm 28.550.000 VNĐ
Năng suất thực thu 2,733 tấn/ha 3,127 tấn/ha Tăng 0,394 tấn/ha
Doanh thu thóc nếp (25.000đ/kg) 68.325.000 VNĐ 78.175.000 VNĐ Tăng 9.850.000 VNĐ
Lợi nhuận ròng gia tăng - - +38.400.000 VNĐ/ha

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

Hạn chế kỹ thuật

  • Phụ thuộc ẩm độ đất: Hiệu lực ức chế mầm cỏ của Atrazine (Mizin 80WP) giảm khi đất quá khô hạn kéo dài sau gieo.
  • Rủi ro trôi rửa: Trên sườn dốc $> 15%$, nếu gặp mưa lớn ngay sau phun có thể gây rửa trôi hoạt chất, làm giảm hoạt tính cục bộ.
  • Nguy cơ kháng thuốc: Việc lạm dụng Glyphosate liên tục nhiều vụ có thể dẫn đến chọn lọc các loài cỏ dại kháng thuốc.

Hướng phát triển tiếp theo

  • Thử nghiệm các chế phẩm thuốc trừ cỏ sinh học có nguồn gốc chiết xuất thực vật (Allelopathy) nhằm tiến tới sản xuất lúa nếp cạn hữu cơ.
  • Nghiên cứu cơ giới hóa khâu phun thuốc bằng thiết bị bay không người lái (UAV/Drone) tại các địa hình nương dẫy khó tiếp cận.
  • Xây dựng mô hình luân canh lúa nếp Khẩu Nua Trạng với các cây họ đậu (đậu tương, đậu đen) để cắt đứt chu kỳ hạt cỏ tự nhiên trong đất.

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

  • Sinh viên và Học viên Cao học: Tài liệu tham khảo chuẩn mực về phương pháp bố trí thí nghiệm nông học (RCBD), kỹ thuật điều tra sâu bệnh theo tiêu chuẩn IRRI/10 TCN 558-2011 và xử lý thống kê bằng IRRISTAT.
  • Kỹ sư Nông học & Cán bộ Khuyến nông: Bộ quy trình thực nghiệm rõ ràng, có căn cứ số liệu định lượng phục vụ chỉ đạo sản xuất mùa vụ tại các huyện vùng cao.
  • Nông hộ và Hợp tác xã vùng cao: Giải pháp giảm 95% chi phí nhân công làm cỏ, nâng cao năng suất thêm $394\text{ kg thóc/ha}$, tối đa hóa lợi nhuận trồng lúa đặc sản.
  • Nhà khoa học & Viện nghiên cứu: Cơ sở dữ liệu khoa học về đặc tính sinh trưởng, khả năng đẻ nhánh và năng suất giống lúa nếp bản địa Khẩu Nua Trạng.

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

1. Phun thuốc Lyphoxim 40SL trước gieo 7 ngày có làm chết mầm lúa khi gieo không?

Không. Hoạt chất Glyphosate Isopropylamine bị bất hoạt rất nhanh khi tiếp xúc với keo đất và vi sinh vật đất. Sau 7 ngày, thuốc đã hoàn thành quá trình lưu dẫn tiêu diệt cỏ dại và phân hủy, hoàn toàn an toàn cho hạt giống lúa nếp cạn khi tra hạt.

2. Nếu sau khi phun Mizin 80WP gặp trời mưa to thì có cần phun lại không?

Mizin 80WP có tính nội hấp nhanh qua thân mầm và lá cỏ sau 4 - 6 giờ. Nếu mưa to xảy ra trong vòng 2 giờ sau phun, hiệu lực thuốc có thể bị giảm $30 - 40%$. Cần theo dõi sau 5 ngày; nếu cỏ non vẫn mọc dày thì mới tiến hành xử lý bổ sung liều lượng $50%$.

3. Tại sao công thức kết hợp làm cỏ tay sau phun thuốc (CT5) năng suất lại thấp hơn CT4?

Việc làm cỏ bằng tay tại thời điểm 45 NSG ở CT5 vô tình làm đứt một số rễ con ăn nông của lúa cạn đang trong giai đoạn chuẩn bị làm đòng, đồng thời làm xáo trộn lớp màng hoạt chất Mizin đang duy trì trên mặt đất, dẫn đến số hạt chắc/bông đạt 49,13 hạt (thấp hơn mức 52,40 hạt ở CT4).

4. Thuốc Mizin 80WP có gây độc cho vụ cây trồng sau (như ngô, đậu) không?

Atrazine trong Mizin 80WP ở liều lượng $500\text{ g/ha}$ là liều lượng thấp, thời gian phân hủy sinh học trong đất từ 45 - 60 ngày. Với chu kỳ sinh trưởng của lúa Khẩu Nua Trạng là 134 ngày, lượng tồn dư đã phân giải hoàn toàn trước khi thu hoạch, tuyệt đối an toàn cho cây trồng vụ sau.

5. Giống lúa Khẩu Nua Trạng có thể thâm canh trên ruộng nước được không?

Có. Nguồn gốc của lúa cạn tiến hóa từ lúa nước thích nghi với điều kiện khô hạn. Giống lúa Khẩu Nua Trạng hoàn toàn có thể sinh trưởng trên chân ruộng nước có kiểm soát, tuy nhiên cần chú ý mật độ và lượng phân đạm vì cây có chiều cao lớn ($> 130\text{ cm}$), dễ bị đổ ngã nếu bón thừa đạm.


Kết luận

Đề tài "Nghiên cứu ảnh hưởng của các phương thức trừ cỏ đến sinh trưởng và phát triển đối với giống lúa nếp cạn Khẩu Nua Trạng Hà Giang tại tỉnh Thái Nguyên" đã giải quyết thành công bài toán thâm dụng lao động và áp lực dịch hại trong canh tác lúa cạn.

Áp dụng công thức trừ cỏ tích hợp (CT4: Phun Lyphoxim 40SL $4\text{ lít/ha}$ trước gieo 7 ngày kết hợp Mizin 80WP $500\text{ g/ha}$ sau gieo 25 ngày) giúp lúa đẻ nhánh tối đa đạt $12,13\text{ nhánh/khóm}$, nhánh hữu hiệu đạt $6,20\text{ nhánh/khóm}$, giảm áp lực sâu đục thân xuống điểm 3, và nâng cao năng suất thực thu lên 31,27 tạ/ha (tăng 14,42% so với canh tác truyền thống).

Đây là quy trình kỹ thuật tiến bộ, có tính ứng dụng thực tiễn cao, cần được nhân rộng trong các chương trình chuyển giao kỹ thuật nông nghiệp tại các tỉnh miền núi phía Bắc nhằm bảo tồn và phát triển các giống lúa đặc sản bản địa giá trị cao.