中国农业科学 ›› 2026, Vol. 59 ›› Issue (14): 3105-3120.doi: 10.3864/j.issn.0578-1752.2026.14.009

• 土壤肥料·节水灌溉·农业生态环境 • 上一篇    下一篇

覆膜微量滴灌条件下旱地谷子临界氮稀释曲线模型构建与氮营养诊断

刘家璇1(), 曹红霞1(), 任利宇2, 马丽娜3, 张斌2, 沈佳琪1, 李志军1   

  1. 1 西北农林科技大学旱区农业水土工程教育部重点实验室, 陕西杨凌 712100
    2 陕西崇仁水利工程有限公司, 西安 710000
    3 榆林农业技术服务中心, 陕西榆林 719000
  • 收稿日期:2025-09-02 接受日期:2025-10-21 出版日期:2026-07-16 发布日期:2026-07-21
  • 通信作者:
    曹红霞,E-mail:
  • 联系方式: 刘家璇,E-mail:liujiaxuan202203@163.com。
  • 基金资助:
    陕西省重点产业创新链-农业领域项目(2023-ZDLNY-46); 陕西省科技厅重点研发项目(2024NC-YBXM-025)

Construction of a Critical Nitrogen Dilution Curve Model and Nitrogen Nutrition Diagnosis for Dryland Foxtail Millet Under Film Mulching with Micro-Drip Irrigation

LIU JiaXuan1(), CAO HongXia1(), REN LiYu2, MA LiNa3, ZHANG Bin2, SHEN JiaQi1, LI ZhiJun1   

  1. 1 Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas of Ministry of Education, Northwest A&F University, Yangling 712100, Shaanxi
    2 Shaanxi Chongren Water Conservancy Engineering Co., Ltd., Xi’an 710000
    3 Yulin Agricultural Technology Service Center, Yulin 719000, Shaanxi
  • Received:2025-09-02 Accepted:2025-10-21 Published:2026-07-16 Online:2026-07-21

摘要:

【目的】明确“四位一体”技术下覆膜微量滴灌条件中不同施肥量和基追比对谷子氮素吸收与地上部生物量累积的影响,构建不同基追比条件下的谷子临界氮稀释曲线,评价贝叶斯分层模型在曲线参数估计中的适用性,并结合氮营养指数与产量响应确定适宜施肥量及基追比,为陕北旱作区谷子水肥一体化条件下精准氮肥管理提供理论依据。【方法】于2023—2024年在陕西榆林开展为期两年的大田滴灌施肥试验,设置B82(基肥﹕追肥=8﹕2)、B64(基肥﹕追肥=6﹕4)和B46(基肥﹕追肥=4﹕6)3个基追比和4个施肥水平(N-P2O5-K2O:0-0-0、90-45-45、120-60-60、150-75-75 kg·hm-2)。研究不同施肥量和基追比对陕北旱作补充灌溉谷子地上部生物量和植株氮浓度的影响,并分别采用经典统计方法和贝叶斯分层模型建立临界氮稀释曲线模型,最终确定最佳的施肥管理策略。【结果】处理B64能够有效促进植株氮素的吸收和地上部生物量的累积,植株地上部生物量和氮浓度总体表现为B64>B82>B46处理。研究发现,不同基追比处理下,临界氮稀释曲线相关参数未表现出统计学显著差异,但降水较少的年份会使曲线稀释曲线参数A2偏大,同时增大了B46处理参数的不确定性。在此基础上,基于贝叶斯统计方法提出2023、2024年不同基追比处理的临界氮稀释曲线B82:Nc=3.58W-0.48、Nc=3.14W-0.36;B64:Nc=3.51W-0.40、Nc=3.65W-0.37;B46:Nc=3.79W-0.62、Nc=2.71W-0.36。基于曲线计算的氮营养指数(NNI)随施肥量增加而升高,NNI在120-60-60 kg·hm-2施肥水平下接近1,此时相对产量(RY)达到最大。综合生物量和产量分析,本试验条件下谷子的最佳施肥量为120-60-60 kg·hm-2,最优基追比为6﹕4。【结论】采用贝叶斯与经典统计构建谷子不同基追比临界氮稀释曲线,贝叶斯估算不确定性更低、精度更优。综合多指标确定当地谷子最佳基追比6﹕4,氮磷钾120-60-60 kg·hm⁻2为适宜施肥量。

关键词: 谷子, 临界氮浓度, 基追比, 氮营养指数, 地上部生物量, 临界氮稀释曲线, 贝叶斯统计模型

Abstract:

【Objective】To quantify the effects of different fertilization rates and basal-to-topdressing ratios on nitrogen uptake and aboveground biomass accumulation of foxtail millet under film-mulched micro-drip irrigation within the “four-in-one” system in northern Shaanxi; to establish critical nitrogen dilution curves under different basal-to-topdressing ratios; to evaluate the applicability of a Bayesian hierarchical model for parameter estimation of the curves; and to determine the optimal fertilization rate and basal-to-topdressing ratio based on the nitrogen nutrition index and yield response, thereby providing a theoretical basis for precise nitrogen management under integrated water-fertilizer conditions in dryland regions of northern Shaanxi.【Method】A two-year field experiment (2023-2024) was conducted in Yulin, Shaanxi Province, using drip irrigation and fertigation treatments. Three basal-to-topdressing fertilizer ratios were set: B82 (8:2), B64 (6:4), and B46 (4:6), along with four fertilizer application levels of N-P2O5-K2O (0-0-0, 90-45-45, 120-60-60, and 150-75-75 kg·hm-2). The effects of NPK rates and basal-to-topdressing ratios on the aboveground dry matter and nitrogen concentration of millet under supplemental irrigation were investigated. Classical statistical analysis and Bayesian hierarchical modeling were employed to establish critical nitrogen dilution curves and to determine optimal fertilization strategies.【Result】The B64 treatment effectively promoted nitrogen uptake and aboveground biomass accumulation, with overall trends of B64>B82>B46 treatment for both dry matter and nitrogen concentration. Although the basal-to-topdressing ratio had no statistically significant effect on the parameters of the critical nitrogen dilution curve, years with lower precipitation resulted in higher dilution coefficient (A2) values and increased parameter uncertainty in the B46 treatment. On this basis, the critical nitrogen dilution curves under different basis tracking treatments in 2023, 2024 were proposed: B82: Nc=3.58W-0.48, Nc=3.14W-0.36; B64: Nc=3.51W-0.40, Nc=3.65W-0.37; and B46: Nc=3.79W-0.62, Nc=2.71W-0.36 under different basis tracking treatments. The nitrogen nutrition index (NNI) calculated from these curves increased with fertilization rate, approaching 1.0 at 120-60-60 kg·hm-2, where relative yield (RY) also reached its maximum. Considering both biomass and yield performance, the optimal fertilizer application rate for millet under these conditions was 120-60-60 kg·hm-2, with the best basal-to-topdressing ratio of 6:4.【Conclusion】Critical nitrogen dilution curves for foxtail millet under different basal-to-topdressing ratios in northern Shaanxi were developed using both Bayesian and classical statistical methods. The Bayesian approach yielded lower uncertainty and higher precision. By integrating multiple indicators, the optimal basal-to-topdressing ratio was determined to be 6:4, with an N-P2O5-K2O rate of 120-60-60 kg·hm-2 recommended as the suitable fertilization regime for local foxtail millet production.

Key words: foxtail millet, critical nitrogen concentration, basal-to-topdressing ratio, nitrogen nutrition index (NNI), aboveground biomass, critical nitrogen dilution curve (CNDC), Bayesian statistical model