1. Introduction

Water and fertilizer are the two most critical factors affecting wheat yield (Wang et al. 2014; Sui et al. 2015;Mon et al. 2016). In the North China Plain (NCP), one of the most important agricultural areas in China (Sun et al.2015), the winter wheat growing season extends from approximately October to the following June. During this period, an average of 120 mm of precipitation in the NCP, is considerably less than the required amount for winter wheat,i.e., 400–500 mm (Dong et al. 2011). Thus, to achieve a high grain yield in the NCP, there has been over-extraction of groundwater, which has resulted in a rapid decline of the groundwater table in the region (Moiwo et al. 2010; Sun et al. 2015). Traditional high-yield cultivation practices for wheat in the NCP include full irrigation and a large input of N fertilizer. Three or more irrigations per wheat cycle is still a routine practice in this region. However, in the NCP, sharp rising temperatures in late May and early June shorten the grain- filling period of winter wheat, resulting in a shortfall of potential yield under full irrigation (Shao et al.2011). In addition, farmers rely heavily on high rates of N fertilization in their attempts to obtain high crop yields(Sui et al. 2015). However, excess N fertilization does not contribute to a significant increase in crop grain yields (Ju et al. 2009; Coventry et al. 2011; Sui et al. 2015). Survey data show that the amount of fertilizer (particularly nitrogen fertilizer) used remains at a high level in the NCP, with an average of 500 kg N ha–1 per year, of which fertilization of winter wheat accounts for over 55% (Zhang et al. 2011). In the Beijing region of the NCP, nitrogen application of 290 kg ha−1 is the local conventional level of fertilization (Sui et al.2015). Over-extraction of groundwater and excessive use of nitrogen fertilizer cause serious environmental problems,such as ground subsidence, degradation of soil, nitrate pollution of groundwater, acid rain, ammonia volatilization,and other forms of air pollution (Zhang et al. 2003; Yang et al. 2006; Ju et al. 2009; Moiwo et al. 2010; Miao et al.2011; Duan et al. 2014; Sun et al. 2015).

那堂妹妹还没问清客人的来由，就往外跑，说是去找她们的祖父去，请他等一等。大概她想是凡男客就是来会祖父的。

中国海油节能管理的实践，不仅是完成国家约束性指标下的自我进步与更新，还有效提升了企业管理水平，更是成为集团公司低成本发展的内在发力点和降本增效的持续推动力，为实现建设有中国特色的国际一流能源公司和集团公司的可持续发展贡献了价值。

With an increase in irrigation and fertilizer application,wheat water use efficiency and agronomic efficiency are generally decreased and grain yield is not increased (Qiu et al. 2008; Zhang et al. 2010, 2011; Shao et al. 2011).Two rounds of irrigation (60 mm each time), one at the jointing and the other at the anthesis stage, during the wheat-growing season have been found to be associated with the highest water use efficiency, and resulted in a grain yield that was not significantly different to that obtained with three rounds of irrigation (Chu et al. 2009). The application of 120–240 kg N ha−1, in combination with traditional flood irrigation methods, is sufficient to maintain high yields of winter wheat (Ju et al. 2009; Sui et al. 2015). Under experimental conditions, the wheat yield did not increase significantly at N rates above 200 kg N ha–1 (Fang et al.2006). However, these previous conclusions are mainly based on short-term (1–2 year) studies, which are not always consistent. Weather conditions during most studies were normal, and there have been fewer studies that have investigated cropping in successive drought years. Studies on nitrogen fertilizer management in the NCP have generally been carried out under traditional irrigation with relatively large irrigation amount and high irrigation frequencies. In addition, the optimal amounts of irrigation and nitrogen in past studies may not be the same as those currently used at the local level, particularly under the circumstance of accumulated soil fertility that has resulted from excessive fertilizer use (Fang et al. 2006; Zhang et al. 2011; Wang et al. 2015).

我从脑海里赶走家人的身影，把思绪拉回现实。双脚分开，与肩同宽，小心翼翼地抓住枪柄。枪很重，很难举起来，但我希望它离我的脸越远越好。我扣住扳机，开始很迟疑，后来慢慢用力，畏畏缩缩扣了下去。砰的一声，我的耳朵被震得很疼，双手在后坐力的作用下几乎打到鼻子上。我踉跄着，急忙双手向后按在身后的墙上，否则肯定会摔个仰面朝天。子弹不知飞到哪里去了，但我很肯定它离靶心不近。

WUEy varied among years and treatments (Tables 3 and 5). With regards to irrigation treatments, the WUEy values of the W1 and W2 groups were higher than those of the W3 groups, and some of the differences reached a significant level during the four years. There were no marked differences between WUEy values of the W1 and W2 groups.In terms of N application, with the exception of N0, WUEy did not differ significantly among N fertilization treatments within the same growing season. Lower values of WUEy were obtained in the N0 treatments than in N-fertilized treatments under the same number of irrigations. In all years, the WUEy value in the W2N1 treatment were performed at the highest level, being 25.4, 29.6, 34.4, and 28.5 kg ha–1 mm–1 for the 2012–2013, 2013–2014, 2014–2015, and 2015–2016 seasons, respectively.

2. Materials and methods

2.1. Experimental site

Water use efficiency at the grain yield level (WUEy) was de fined as the grain yield (kg ha–1) divided by ET (mm), and water use efficiency at the BM level (WUEbm) was calculated as BM (kg ha–1) divided by ET (mm):

“初日照杨柳，玉楼含翠阴。啼春独鸟思，望远佳人心。幽怨贮瑶瑟，韶光凝碧林。所思曾不见，芳草意空深。”[5]

Soil at the study site was a loam with a field capacity of 38% (v/v), a permanent wilting point of 13% (v/v), and a bulk density of 1.53 g cm–3 for a 2-m pro file. Soil fertility was measured from the soil surface to a depth of 30 cm prior to the first growing season (2012–2013). Soil organic matter was 16.3 g kg–1, total nitrogen was 1.2 g kg–1,and available phosphorus and potassium were 38.0 and 158.5 mg kg–1, respectively.

2.2. Experimental design

AE was significantly affected by both irrigation and N(Table 3). At the same N level, AE generally increased with an increasing number of irrigation events (Table 6). At all irrigation levels, the rate of fertilizer application negatively affected AE, AE decreased with increasing fertilization rates.

2.3. Measurements

?

Total water use or evapotranspiration (ET) during the winter wheat growing season was calculated using the following waterbalance formula modified from that described by Zhang et al. (2012):

从中国 《宪法》《教育法》《公共文化服务保障法》和《公共图书馆法》看，公民享受图书馆服务已是不争事实，公民利用图书馆是应该享有的基本权利，已经以法律化确认下来的。

Soil moisture and water use efficiency Neutron access tubes were installed in each plot to measure soil water content. Soil volumetric water content was measured using a neutron probe(CNC503B; Beijing Aozuo Ecology Instrumentation Ltd., China)down to a soil depth of 160 cm at 20-cm interval before sowing and after harvest. The water content in the top 20 cm of soil was measured using the oven-drying method.

Plant height, aboveground biomass, grain yield, and yield components At harvest, 20 winter wheat plants were selected at random from each plot to measure the plot’s average plant height.Plants were harvested manually and three 1-m2 samples were selected at random from each plot. All plant samples were air-dried to a grain water content of approximately 11%. The spike number(SN) and aboveground biomass (BM) per 1 m2 were recorded.Plants were then threshed using a stationary thresher to obtain the grain yield. The 1 000-grain weight (GW) was measured using the average of three samples of 1 000 grains. The grain number (GN)per spike was estimated using SN, GW, and grain yield.

由图5知，金矿细磨至-400目占85%所需细磨时间为10min，结合矿物细磨过程中能耗变化曲线，在立磨机实际处理量为500t/d的情况下，通过计算可知选用装机功率为280kW、容积为10.0m3的KLM-280立磨机可满足要求。

Table 2 Experimental design

W1N0 1 Jointing 60 0 0 W1N1 1 Jointing 60 100 172 W1N2 1 Jointing 60 200 172 W1N3 1 Jointing 60 300 172 W2N0 2 Jointing, heading 120 0 0 W2N1 2 Jointing, heading 120 100 172 W2N2 2 Jointing, heading 120 200 172 W2N3 2 Jointing, heading 120 300 172 W3N0 3 Jointing, heading, filling 180 0 0 W3N1 3 Jointing, heading, filling 180 100 172 W3N2 3 Jointing, heading, filling 180 200 172 W3N3 3 Jointing, heading, filling 180 300 172

Treatment1) Irrigation time Irrigation stage Irrigation amount (mm) Total fertilization amount (kg ha−1)N P

Where, ET (mm) is the total water consumption during a given growth period, I (mm) is the amount of irrigation,P (mm) is precipitation, and ΔSWC (mm) is the initial soil water content minus final soil water content.

Field experiments were carried out during four growing seasons from 2012 to 2016 at the Luancheng Agro-Ecosystem Experimental Station of the Chinese Academy of Sciences (37°53´N, 114°41´E; elevation 50.1 m).Luancheng Agro-Ecosystem Experimental Station lies in the northern part of the NCP, which has a continental monsoon climate, and is representative of typical NCP high-production areas. Mean annual rainfall in the region is 482 mm, with more than 70% of precipitation concentrated in the summer season. Mean precipitation during the winter wheat growing season for the period 1990 to 2010 was approximately 120 mm (Wu et al. 2012). During the present study,precipitation in the first, second, third, and fourth growing seasons was 92, 49, 77, and 70 mm, respectively. According to Fu et al. (2014), all four experimental growing seasons would be considered drought years. Meteorological data were recorded at a weather station approximately 500 m from the experiment field (Table 1).

安：我在克利夫兰音乐学院教过不少中国学生，在中国也为很多学生上过大师课。我认为，中国学生大多十分认真、敬业，也十分好学。但绝大多数中国学生的理念存在问题，他们会把自己认为的重点作为训练的终极目标—如速度、力度，以及准确性，等等，往往忽视对音色的探索，也并没有把表达内心作为一项重要的任务来看待。当然这只是对群体的描述，我也遇到过很多极有音乐性的中国学生。但大体而论，中国学生把“弹快、弹准”看得比音乐性更加重要，他们需要理解帕内拉与齐科里尼教授所言的“音乐意识”。

Agronomic efficiency Agronomic efficiency (AE) is an index used to estimate the grain yield increase from fertilizer application (Duan et al. 2014). AE was calculated using the following formula (Mon et al. 2016):

Where, YF (kg ha−1) is the grain yield in a fertilized plot,Y0 (kg ha−1) is the grain yield in an unfertilized plot, and AF(kg ha−1) is the total amount of fertilizer applied to the plot.

Statistical analysis Data were analyzed using an analysis of variance (ANOVA) to determine differences between treatments. ANOVA was performed using IBM SPSS 22.Differences among mean values were calculated using least significance differences (LSD) at the 5% level.

3. Results

3.1. Plant height and aboveground biomass

As the yield of winter wheat is positively correlated with SN, GW, and GN, fertilizer and irrigation can improve yield by regulating these components of yield. Measurements associated with yield components (SN, GW, and GN) for each treatment group are presented in Table 4. Irrigation level had a significant effect on all yield components(Table 3). Although SN and GW were significantly affected by N level, GN did not show a similar response (Table 3).F-values were larger for GW than for SN and GN, and irrigation had a larger effect than N (Table 3). As shown in Table 4, GW increased significantly with an increasing number of irrigations. In the 2012–2013 growing season, the mean GWs for the W1, W2, and W3 groups were 29.5, 33.2,and 35.1 g, respectively. The corresponding values in the 2013–2014, 2014–2015, and 2015–2016 growing season were 39.5, 40.5, and 42.9 g; 32.3, 37.5, and 38.6 g; and 32.0,32.1, and 34.1 g, respectively. Lv et al. (2013) indicated that sunshine duration plays a key role in determining grain weight, and that longer daylight hours at the filling stage are beneficial for an enhancement in grain weight. In the present study, the values of GW in the 2013–2014 season ranked the highest, followed by those in 2014–2015, 2012–2013,and 2015–2016, with 8.7, 8.4, 6.3, and 6.3 h of sunshine in May, respectively (Table 1). The W2N1 treatment resulted in higher SN and GN and moderate GW, and ultimately produced a higher yield. This result was likely mainly attributable to the higher mobilization efficiency and selfadjustment mechanism among the three yield components of the grain (Zhang et al. 2008; Lv et al. 2013).

Investigation and Analysis of Wind Energy Resources on Global Ocean Shipping Lines

Total BM of the wheat varied among the different irrigation regimes, N levels, and years (Table 3; Fig. 2). BM varied with irrigation level under each N level, except N0 in 2015–2016 (Fig. 2). Total BM under W2 and W3 conditions was similar and higher than that in the W1 treatment for the same N treatment. On average, over the four growing seasons,total BM was 13.1, 16.9, 15.6, and 12.8% greater under W2 than that under W1 for N0, N1, N2, and N3, respectively.The BM of the N0 treatment was significantly lower than that of all three N fertilization treatments. N fertilization resulted in an increase in BM by 14.5–55.7% over the N0 treatment.Maximum BM differed across the 4 years with values of 15 509, 16 337, 17 193, and 18 315 kg ha–1 in 2013, 2014,2015, and 2016, respectively (Fig. 2).

3.2. Grain yield and yield components

Grain yields among the different treatment groups showed significantly different responses to irrigation and N levels in all four growing seasons (Fig. 3). F-values were similar among irrigation and N treatments, and the N×Irrigation,Year×Irrigation, and Year×N interactions were alsosignificant (Table 3). The yield of the W1 group was very low, regardless of N level, and was significantly lower than that of the W2 and W3 groups. The yield of W2 was high and similar to that of W3. The yield was greater with the addition of N than with the N0 treatment; however, there were no significant differences in grain yield among the different levels of N fertilization. Overall, grain yield showed a year-after-year rise over the course of the investigation.The pattern of variation in yield was similar for all 4 years.The higher grain yields for the years 2013, 2014, 2015, and 2016 (7 500, 8 800, 8 800, and 9 200 kg ha–1, respectively)were achieved using the W2N1 treatment (Fig. 3).

Table 3 F-values of analysis of variance combined across year, nitrogen and irrigation on winter wheat height, BM, yield, SN, GW,GN, ET, WUEy, WUEbm and AE1)

1) BM, aboveground biomass; SN, spike number; GW, 1 000-grain weight; GN, grain number per spike; ET, evapotranspiration; WUEy,water use efficiency at the grain yield level; WUEbm, water use efficiency at the biomass level; AE, agronomic efficiency. ＊, significant at the 0.05 probability level; ＊＊, significant at the 0.01 probability level; NS, non-significant.

Source of variation Height BM Yield SN GW GN ET WUEy WUEbm AE Year (Y) 101.7＊＊ 59.5＊＊ 105.7＊＊ 93.6＊＊ 266.9＊＊ 67.8＊＊ 36.8＊＊ 26.1＊＊ 67.2＊＊ 52.6＊＊Nitrogen (N) 69.1＊＊ 190.6＊＊ 224.4＊＊ 29.6＊＊ 70.6＊＊ NS 17.1＊＊ 40.2＊＊ 69.6＊＊ 49.5＊＊Irrigation (I) 153.1＊ 156.5＊＊ 225.5＊＊ 8.1＊＊ 108.0＊＊ 41.2＊＊ 235.8＊＊ 20.8＊＊ 20.8＊＊ 26.8＊＊N×I 3.0＊ NS 2.9＊ NS 9.4＊＊ NS NS NS NS NS Y×I 12.4＊＊ 7.2＊＊ 12.9＊＊ NS 9.2＊＊ 3.8＊＊ NS 4.7＊＊ 10.7＊＊ 2.6＊Y×N 2.8＊＊ 2.6＊＊ 5.5＊＊ 3.8＊＊ 8.0＊＊ NS 10.5＊＊ 3.9＊＊ 7.3＊＊ NS Y×N×I NS NS NS NS NS NS NS NS NS NS

Fig. 1 Plant height in 2012–2013 (A), 2013–2014 (B), 2014–2015 (C), and 2015–2016 (D) for winter wheat under different irrigation and nitrogen application conditions. W1 indicates irrigation once at jointing stage; W2 indicates irrigation once at jointing and once at heading stage; W3 indicates irrigation once at jointing, once at heading, and once at filling stage. N0, N1, N2, and N3 represent nitrogen fertilization at 0, 100, 200, and 300 kg N ha−1, respectively. Different letters indicate significant differences at P=0.05.Error bars indicate standard deviation.

Fig. 2 Aboveground biomass in 2012–2013 (A), 2013–2014 (B), 2014–2015 (C), and 2015–2016 (D) for winter wheat under different irrigation and nitrogen application conditions. W1 indicates irrigation once at jointing stage; W2 indicates irrigation once at jointing and once at heading stage; W3 indicates irrigation once at jointing, once at heading, and once at filling stage. N0, N1, N2, and N3 represent nitrogen fertilization at 0, 100, 200, and 300 kg N ha−1, respectively. Different letters indicate significant differences at P=0.05. Error bars indicate standard deviation.

Fig. 3 Grain yield in 2012–2013 (A), 2013–2014 (B), 2014–2015 (C), and 2015–2016 (D) under different irrigation and nitrogen application conditions. W1 indicates irrigation once at jointing stage; W2 indicates irrigation once at jointing and once at heading stage; W3 indicates irrigation once at jointing, once at heading, and once at filling stage. N0, N1, N2, and N3 represent nitrogen fertilization at 0, 100, 200, and 300 kg N ha−1, respectively. Different letters indicate significant differences at P=0.05. Error bars indicate standard deviation.

Table 4 Yield component factors of winter wheat as affected by irrigation and nitrogen during four growing seasons

1) W1 indicates irrigation once at jointing stage; W2 indicates irrigation once at jointing and once at heading stage; W3 indicates irrigation once at jointing, once at heading, and once at filling stage. N0, N1, N2, and N3 represent nitrogen fertilization at 0, 100, 200, and 300 kg N ha−1, respectively. 2) SN represents spike number, GW represents 1 000-grain weight, GN represents grain number per spike. Within each column, different letters indicate significant differences between measurements at P=0.05, based on LSD tests.

GN (no.spike–1)W1N0 547.7 d 28.9 f 31.8 d 443.2 e 36.4 e 35.5 cdef 521.1 b 32.1 c 32.0 c 563.3 b 31.6 e 36.5 d W1N1 621.5 bc 29.0 f 34.4 abc 497.7 bcde 42.0 c 35.0 def 611.6 ab 32.5 c 32.4 bc 715.6 a 31.8 de 36.6 cd W1N2 643.4 ab 30.2 def32.1 cd 527.2 abcd 40.5 cd 34.5 ef 651.4 a 31.9 c 32.2 bc 688.5 a 32.6 d 37.2 bcd W1N3 638.3 ab 29.8 ef 32.1 cd 577.4 a 39.0 d 33.0 f 625.6 a 32.9 c 33.8 abc 708.2 a 32.0 cd 37.1 bcd W2N0 575.8 cd 31.9 cde 34.9 ab 498.6 bcde 34.4 e 37.0 bcde 647.8 a 34.0 c 33.7 abc 583.5 b 30.4 cd 36.9 bcd W2N1 628.6 abc 34.0 bc 35.4 a 550.7 ab 42.4 c 37.8 bcde 660.1 a 39.4 ab 34.0 abc 742.8 a 32.2 cd 38.6 a W2N2 685.1 a 33.7 bc 35.1 a 554.2 ab 42.7 bc 37.1 bc 668.9 a 38.3 b 34.6 abc 726.6 a 32.5 cd 37.9 abc W2N3 647.9 ab 33.3 bc 35.2 a 538.6 abc 42.7 bc 38.2 bc 674.1 a 38.3 b 34.4 abc 714.6 a 33.1 bc 38.0 ab W3N0 555.1 d 32.4 cd 34.9 ab 462.8 de 35.6 e 37.9 bc 694.9 a 31.0 c 34.1 abc 538.1 b 32.2 cd 37.9 abc W3N1 615.7 bc 36.9 a 34.3 abc 515.4 abcd 45.7 a 37.4 bcd 602.0 ab 41.8 a 35.3 a 688.8 a 35.0 a 38.0 ab W3N2 604.1 bcd 35.4 ab 35.1 a 506.6 bcde 45.0 ab 39.3 ab 668.9 a 40.4 ab 35.0 ab 717.6 a 35.0 a 38.0 ab W3N3 642.6 ab 35.5 ab 32.6 bcd 476.3 cde 45.5 a 41.4 a 630.5 a 41.0 ab 35.8 a 712.2 a 34.0 ab 39.0 a Treatment1)2012–20132) 2013–20142) 2014–20152) 2015–20162)SN(×104 ha−1)GW(g)GN (no.spike–1)SN(×104 ha−1)GW(g)GN (no.spike–1)SN(×104 ha−1)GW(g)GN (no.spike–1)SN(×104 ha−1)GW(g)

Plant height responded to both irrigation and N levels, but the effect of N was not as strong as that of irrigation (Table 3;Fig. 1). For all four years, plant height for a given N level tended to increase with increasing amounts of irrigation. The plant height with treatment W3 was 4.0–7.1%, 2.2–4.9%,5.2–12.9%, and 1.9–2.7% higher than that with W1 for the same N treatment during the 2012–2013, 2013–2014,2014–2015, and 2015–2016 seasons, respectively. The plant height with treatment N0 was significantly lower than that with the three N fertilized treatments. However, plant height within each irrigation level did not differ significantly among N1, N2, and N3. The N×Irrigation interaction for plant height was significant (Table 3). The higher N level,the higher was the plant height with an increasing number of irrigations (Fig. 1). The lower plant heights recorded for the 2014–2015 season may probably be attributed to low precipitation before April 2015 (11 mm).

3.3. Evapotranspiration and water use efficiency

ET was more strongly affected by irrigation than by N treatment (Table 3). ET increased significantly with the number of irrigation events per season (Fig. 4). With the exception of the 2014–2015 growing season, N fertilization increased ET to a greater extent than that observed in the N0 treatment. However, there was no significant difference among N-fertilization treatments. Among the different irrigation treatments, the highest ET values were obtained in the W3 groups for each N treatment and growing season.

The aim of this study was to define a suitable irrigation regime and fertilization rates for wheat production in the drought-prone NCP. We performed experiments over four years, focusing on the effects of different irrigation approaches and fertilizer application regimes on grain yield,water use efficiency, fertilizer agronomic efficiency, and economic benefit under field conditions. We hypothesized that amounts of irrigation and nitrogen applied at the local level currently are excessive, and lower amounts would be satisfactory for wheat growth, improving yield, as well as have economic and ecological bene fits.

The pattern of differences in WUEbm among treatment groups was nearly identical to that for WUEy (Tables 3 and 5). With an increasing number of irrigations, there appeared to be a downward trend in WUEbm. With the exception of N0, no significant differences were observed among N treatments. There was a clear linear correlation between WUEy and WUEbm over the four-year experimental period(Fig. 5). These results indicate that the high WUEbm of winter wheat may contribute to high WUEy, and this may ultimately result in a high yield.

3.4. Agronomic efficiency

In this study, the winter wheat cultivar Shixin 828 was used in a maize-wheat rotation system over four growing seasons (2012–2016). It is an efficient and high-quality winter wheat cultivar widespread extended in northern China (Xu et al.2013). The experiment was carried out during the winter wheat growing seasons for the years 2012–2016 in the same field using a completely randomized block design. The size of each experimental plot was 4 m×6 m. The experiment included 12 treatments that consisted of each combination of four levels of fertilization and three levels of irrigation. Each treatment was replicated three times. The application of 300 kg N ha−1 and three rounds of irrigations is the local conventional practice. The fertilization treatments used in the present study were 0, 100, 200, and 300 kg N ha−1 for N0, N1, N2,and N3, respectively. All the treatments incorporated fertilization with 172 kg P ha−1. Nitrogen fertilizer was added to the soil of each plot in two stages: 70% was added at sowing and 30% at the jointing stage. P was applied as 100% basal. Detailed protocols are shown in Table 2. The three irrigation treatments were one(W1), two (W2), and three (W3) irrigations based on the wheat development stage (Table 2). For each irrigation, 60 mm of water was applied to the soil by border irrigation using a 7-inch hose pipe connected to a pumping well. Water meters were used to record the irrigation applied to each plot. In each year, pre-sowing water was applied to ensure the same soil water content before sowing.Seeds were sown on October 9 using a 2BJM Bevel-type Precision Seeder (Laiwu Hualong Machinery Factory, China) at a seeding rate of 195 kg ha–1 and row spacing of 15 cm. Winter wheat was harvested on June 8th, 7th, 11th, and 10th in 2013, 2014, 2015,and 2016, respectively. Other field management practices were conducted using the same standards as those used in local fields to ensure good crop growth throughout the growing season.

Fig. 4 Evapotranspiration (ET) in 2012–2013 (A), 2013–2014 (B), 2014–2015 (C), and 2015–2016 (D) for winter wheat under different irrigation and nitrogen application conditions. W1 indicates irrigation once at jointing stage; W2 indicates irrigation once at jointing and once at heading stage; W3 indicates irrigation once at jointing, once at heading, and once at filling stage. N0, N1,N2, and N3 represent nitrogen fertilization at 0, 100, 200, and 300 kg N ha−1, respectively. Different letters indicate significant differences at P=0.05. Error bars indicate standard deviation.

Table 5 Water use efficiency at the grain yield level (WUEy) and at the aboveground biomass level (WUEbm) of winter wheat under different irrigation and nitrogen treatments during four growing seasons

1) W1 indicates irrigation once at jointing stage; W2 indicates irrigation once at jointing and once at heading stage; and W3 indicates irrigation once at jointing, once at heading, and once at filling stage. N0, N1, N2, and N3 represent nitrogen fertilization at 0, 100, 200,and 300 kg N ha−1, respectively. Within each column, different letters indicate that differences between measurements are significant at the level of P=0.05, based on the LSD test.

WUEy(kg ha–1 mm–1)W1N0 22.3 de 47.1 cde 23.5 d 46.5 cd 22.7 f 45.5 d 28.8 abc 57.0 ab W1N1 24.3 abc 51.6 abc 27.1 bc 51.4 abc 29.2 cde 56.8 bc 29.9 ab 56.6 abc W1N2 24.6 abc 52.0 ab 29.2 ab 54.7 a 28.3 de 55.8 bc 30.2 a 58.1 ab W1N3 24.6 abc 54.0 a 29.4 ab 55.6 a 29.5 bcde 58.4 ab 30.6 a 59.9 a W2N0 24.0 bcd 46.5 de 20.7 e 41.3 de 27.4 e 49.8 cd 25.1 de 50.0 de W2N1 25.4 abc 50.7 abcd 29.6 ab 51.5 abc 34.4 a 64.4 a 28.5 abc 54.1 abcd W2N2 25.8 ab 50.6 abcd 29.9 a 53.8 ab 32.3 abc 60.4 ab 28.6 abc 54.6 abcd W2N3 26.2 a 54.1 a 29.8 a 55.0 a 32.9 ab 60.4 ab 27.6 abcd 53.1 bcd W3N0 21.9 e 45.3 e 19.1 e 40.1 e 23.5 f 44.4 d 23.8 e 45.9 e W3N1 24.9 abc 48.8 bcde 27.0 bc 48.6 bc 31.3 abcd 59.1 ab 25.7 cde 48.9 de W3N2 23.8 cde 48.5 bcde 27.8 abc 50.2 abc 31.0 bcd 56.4 bc 26.5 cde 50.8 cde W3N3 22.1 e 45.6 de 26.6 c 48.6 bc 32.3 abc 59.0 ab 27.1 bcd 52.5 bcd Treatment1)2012–2013 2013–2014 2014–2015 2015–2016 WUEy(kg ha–1 mm–1)WUEy(kg ha–1 mm–1)WUEy(kg ha–1 mm–1)WUEy(kg ha–1 mm–1)WUEy(kg ha–1 mm–1)WUEy(kg ha–1 mm–1)WUEy(kg ha–1 mm–1)

The AE values of N1 groups were higher than those of N2 and N3 groups for all irrigation treatments. Averaged over the four growing seasons, AE values under W2 irrigation conditions were 32.6 and 82.1% greater under N1 than under N2 and N3, respectively.

3.5. Economic efficiency and ecological sustainability

Fig. 5 Linear regression between water use efficiency (WUE) at the grain yield level (WUEy) and at the biomass level (WUEbm) of winter wheat in 2012–2013 (A), 2013–2014 (B), 2014–2015 (C), and 2015–2016 (D).

Table 6 Agronomic efficiency (AE) of winter wheat as affected by irrigation and nitrogen application during four growing seasons

1) W1 indicates irrigation once at jointing stage; W2 indicates irrigation once at jointing and once at heading stage; W3 indicates irrigation once at jointing, once at heading, and once at filling stage. N0, N1, N2, and N3 represent nitrogen fertilization at 0, 100, 200, and 300 kg N ha−1, respectively. Within each column, different letters indicate that differences between measurements are significant at the level of P=0.05,based on the LSD test.

?

The economic bene fit of different treatments was analyzed using the average grain yield of the four years (Table 7).Results indicated that within the different N treatments, W2 groups showed the greatest bene fits. With the exception of treatment N0, net income decreased with increasing N fertilization treatments under the same irrigation treatment.Compared with the traditional cultivation practice in local fields (W3N3), the net income of W2N1 was higher by 12.3%, and the output-input ratio was 19.5% higher.

From the standpoint of environmental sustainability, the W2N1 treatment economizes on irrigation and fertilizer usage. Compared with the traditional practices of irrigation and fertilization, using the W2N1 regime could save 600 m3 ha–1 of irrigation water and 200 kg N ha–1 without impairing farmers’ benefits. Given that irrigation water is becoming less available in many areas and that the effects of global climate change (e.g., global warming, haze) are becoming increasingly evident (Jiru and Van Ranst 2010), the W2N1 treatment employed in this experiment meets the acute water challenges facing those who farm in the NCP, ful fills the requirements for crop production, and is environmentally friendly.

4. Discussion

In the present study, we determined that the optimal fertilization rate and irrigation regime for wheat production in the drought-prone NCP was 100 kg N ha−1 coupled with two irrigations of 60 mm each at the jointing and heading stages,respectively. This W2N1 treatment consistently produced the highest yield in all four years of the study. Similarly, Dong et al. (2011) reported that two irrigation events of 75 mm in a drought year would result in optimal WUEy and yield in the NCP. Wang et al. (2015) reported that fertilization with 210–240 kg N ha−1 under supplemental irrigation at the jointing and anthesis stages produced the highest grain yield. In contrast, the results of our study indicate that there were no marked increases in yield with increasing nitrogen fertilizer application in excess of 100 kg N ha–1 under each of the irrigation treatments. Lv et al. (2013)showed that ET increases gradually when irrigation water supply is increased and that ET values with four irrigations were the highest, being 25–113 mm higher than with two irrigations. Consistently, we showed that the ET of the W2N1 treatment was approximately 30 mm lower than that obtained with local traditional practice (W3N3). Excessive application of irrigation water and nitrogen cannot increase crop yields further, and thus represent a considerable waste of resources (Shirazi et al. 2014). Under the conditions in the present study, the optimal W2N1 treatment yielded the maximum bene fits in terms of economizing the use of irrigation water and nitrogen fertilizer.

3个月面授结束之后，学员可回原单位继续临床实践工作，但仍然和指导教师保持联系，由指导教师制定个案书写及开题报告书写要求，并及时回答学员的在这些方面的困惑。督促学员积极查阅文献、书写个案及科研报告，及时交付并适当修正，以备毕业答辩。答辩前2周进行预答辩。

Across all years in the present study, grain yields under W2N1 were similar to those under W3N3. Many studies have shown that biomass and the harvest index (HI) are the two main factors determining grain yield (Sayre et al.1997; Shearman et al. 2005; White and Wilson 2006).The W2N1 treatment produced the same biomass and HI as the W3N3 treatment (W2N1, HI=0.53; W3N3, HI=0.52;P＞0.05), with lower water and nitrogen inputs, which resulted in maintenance of a higher grain yield, WUE, and AE. In the present study, W2N1 resulted in a higher WUEbm and WUEy than W3N3 did. It has been shown that it is possible to simultaneously achieve both satisfactory yield and WUE by regulating the number of irrigations (Lv et al. 2013).Consistent with this, we found positive relationships between WUEy and WUEbm. Hsiao (1993) reported that WUEbm is generally closely related to the photosynthetic WUE of single leaves (short-term gas exchange on a photosynthesis basis), but that WUEy is not closely related to WUEbm. In contrast, Qiu et al. (2008) showed that WUEbm is positively related to WUEy, with increases in WUEbm having positive effects on WUEy, thereby enhancing ultimate yields (Qiu et al. 2008). The results of the present study are consistent with those of Qiu et al. (2008) as we found that WUEy was positively associated with WUEbm. A positive relationship between WUEy and WUEbm could be used to provide valuable information to guide agricultural production and farming practices. By monitoring WUEbm throughout the growing season, and adjusting it by reasonable irrigation,farmers may simultaneously achieve relatively high yield and optimal WUEy. Our results showed that in all the four years we investigated, AE decreased with increasing fertilizer application. Moreover, there was an upward trend in AE with increasing irrigation. These results are consistent with those of Mon et al. (2016), who found that AE was negative at low irrigation levels and high N rates.

Xue et al. (2006) reported that in the Texas High Plain,irrigation increases ET, and that biomass and grain yield also increases as ET increased. However, although our results showed that ET increased as irrigation increased in the NCP,biomass and grain yield did not show similar increases.Although the yield of wheat with two irrigations was not significantly different from that with three irrigations, the ET of two-irrigation wheat was lower than that of three-irrigation wheat. Similar results were obtained by Sun et al. (2011)who indicated that ET tended to increase with an increasing level of irrigation. Increases in irrigation frequency and amount of water might increase soil evaporation. The major factor affecting soil evaporation fluctuation was surface soilmoisture (Gao et al. 2005). Compared with three irrigations,the two-irrigation treatment was associated with less surface soil moisture and lower vapor flux, which decreases soil evaporation. Furthermore, larger amounts of irrigated water may cause unfavorable delayed senescence in winter wheat and thus raise the ET.

Table 7 The output-input ratio of different irrigation and nitrogen treatments

1) W1 indicates irrigation once at jointing stage; W2 indicates irrigation once at jointing and once at heading stage; W3 indicates irrigation once at jointing, once at heading, and once at filling stage. N0, N1, N2, and N3 represent nitrogen fertilization at 0, 100, 200, and 300 kg N ha−1, respectively. 2) Grain yield was the average of grain yield for four years.

Treatment1)Grain yield(kg ha–1)2)Output(CNY ha–1)Input (CNY ha–1) Total input(CNY ha–1)Net income(CNY ha–1)Outputinput ratio Seed Fertilizer Pesticide Electricity for irrigation Machine Labor W1N0 5 649.3 12 428.5 600 0.0 450 300 3 000 900 5 250.0 7 178.5 2.4 W1N1 7 043.5 15 495.8 600 1 541.4 450 300 3 000 900 6 791.4 8 704.4 2.3 W1N2 7 148.0 15 725.6 600 2 019.3 450 300 3 000 900 7 269.3 8 456.3 2.2 W1N3 7 216.2 15 875.6 600 2 497.8 450 300 3 000 900 7 747.8 8 127.8 2.0 W2N0 6 667.2 14 667.7 600 0.0 450 600 3 000 1 275 5 925.0 8 742.7 2.5 W2N1 8 601.7 18 923.8 600 1 541.4 450 600 3 000 1 275 7 466.4 11 457.4 2.5 W2N2 8 622.2 18 968.7 600 2 019.3 450 600 3 000 1 275 7 944.3 11 024.4 2.4 W2N3 8 503.5 18 707.8 600 2 497.8 450 600 3 000 1 275 8 422.8 10 285.0 2.2 W3N0 6 590.5 14 499.0 600 0.0 450 900 3 000 1 650 6 600.0 7 899.0 2.2 W3N1 8 632.6 18 991.8 600 1 541.4 450 900 3 000 1 650 8 141.4 10 850.4 2.3 W3N2 8 852.0 19 474.5 600 2 019.3 450 900 3 000 1 650 8 619.3 10 855.2 2.3 W3N3 8 770.9 19 295.9 600 2 497.8 450 900 3 000 1 650 9 097.8 10 198.1 2.1

An important goal of agricultural practices is to achieve continuous and efficient development of farmland, which takes into consideration both production and environmental bene fits. We found that among the different treatments, the highest net income (11 457 CNY ha–1) was obtained using the W2N1 treatment. The corresponding output-input ratio for this treatment (2.5) was also the highest among those treatments examined.

There remains considerable potential for further savings in the NCP with regards to water and fertilizer use. In the late 1970s and early 1980s, there were as many as 8–10 irrigation events during the winter wheat growing season,which maintained a moist soil surface during the entire growing period. Subsequent research on water-saving irrigation systems for winter wheat, has indicated that the number of irrigation events per growing season could be reduced to five (Zhang 2011). In the 21st century,Zhang et al. (2003, 2008) devised irrigation regimes for high-yielding and water-saving winter wheat cultivation that entailed one irrigation in a wet year, two irrigations in a normal year, and three irrigations in a drought year.In the present study, two irrigation treatments sufficed for wheat production in the drought-prone NCP. This is still potentially higher than that found in previous studies. We found that the difference in grain yield between N0 and N1 application was larger under W2 irrigation than that under W1 irrigation. It has previously been shown that a water deficit can restrict the effectiveness of N fertilizer (Gu et al.2015). Under two irrigations, 100 kg N ha–1 can satisfy nutritional requirements for the period of crop growth.Compared with current practices, 60 mm of irrigated water and 200 kg N ha–1 could be saved without yield reductions in a drought year. In a normal or wet year, the water-saving potential could be improved even further, when rainfall coincides with the key wheat growing period. In addition,the irrigation-saving effect of W2N1 may have important implications for relieving groundwater over-exploitation and improving the ecological environment. In the NCP,the area planted with winter wheat covers over 5 million ha, and on the basis of our experimental results, there is tremendous potential to save 3 billion m3 of irrigation water in this region. Further, the low fertilizer input of the W2N1 treatment may reduce the risk of environmental degradation produced by excess nitrogen including nitrate pollution of groundwater, ammonia volatilization, and other forms of air pollution.

5. Conclusion

In the present study, we found that different irrigation and N fertilization treatments significantly affected grain yield,water use efficiency, and fertilizer agronomic efficiency.This four-year experiment showed that in drought years in the North China Plain: (1) a high yield of 7 500–9 200 kg ha–1 can be achieved with two irrigation events (one at the jointing stage and one at the heading stage) coupled with 100 kg N ha–1; (2) two rounds of irrigation per growth season can produce a higher WUEy compared with the local conventional three rounds of irrigation; and (3) the W2N1 treatment achieved the maximum economic returns and has important implications for the sustainable development of agriculture and the ecological environment.

Acknowledgements

This study was supported by the National Key Research and Development Program of China (2016YFD0300808), the National Key Technologies R&D Program of China during the 12th Five-Year Plan period (2013BAD05B02), the National Natural Science Foundation of China (31571612 and 31100191), the Science and Technology Service Network Initiative of Chinese Academy of Sciences (KFJ-STSZDTP-001), and the Hebei Key Research and Development Program, China (15226407D and 17227006D).

References

Chu P F, Wang D, Zhang Y L, Wang X Y, Wang X Z, Yu Z W. 2009. Effects of irrigation stage and amount on water consumption characteristics, grain yield and content of protein components of wheat. Scientia Agricultura Sinica,42, 1306–1315. (in Chinese)

Coventry D R, Yadav A, Poswal R S, Sharma R K, Gupta R K, Chhokar R S, Gill S C, Kumar V, Kumar A, Mehta A,Kleemann S G L, Cummins J A. 2011. Irrigation and nitrogen scheduling as a requirement for optimising wheat yield and quality in Haryana, India. Field Crops Research, 123, 80–88.

Dong B, Shi L, Shi C, Qiao Y, Liu M, Zhang Z. 2011. Grain yield and water use efficiency of two types of winter wheat cultivars under different water regimes. Agricultural Water Management, 99, 103–110.

Duan Y H, Xu M G, Gao S D, Yang X Y, Huang S M, Liu H B, Wang B. 2014. Nitrogen use efficiency in a wheat-corn cropping system from 15 years of manure and fertilizer applications. Field Crops Research, 157, 47–56.

Fang Q X, Yu Q, Wang E L, Chen Y H, Zhang G L, Wang J,Li L H. 2006. Soil nitrate accumulation, leaching and crop nitrogen use as influenced by fertilization and irrigation in an intensive wheat-maize double cropping system in the North China Plain. Plant and Soil, 284, 335–350.

Fu Q, Wang Q, Shen X, Fan J. 2014. Optimizing water and nitrogen inputs for winter wheat cropping system on the Loess Plateau. Journal of Arid Land, 6, 230–242.

Gao L, Hu C S, Chen S Y, Zhang L F. 2005. Experimental study on soil evaporation in winter wheat field under sprinkling irrigation condition. Transactions of the Chinese Society of Agricultural Engineering, 21, 183–185. (in Chinese)

Gu L M, Liu T N, Zhao J, Dong S T, Liu P, Zhang J W, Zhao B.2015. Nitrate leaching of winter wheat grown in lysimeters as affected by fertilizers and irrigation on the North China Plain. Journal of Integrative Agriculture, 14, 374–388.

Hsiao T C. 1993. Effects of drought and elevated CO2 on plant water use efficiency and productivity. In: Jackson M B, Black C R, eds., Interacting Stresses on Plants in a Changing Climate. Springer, Berlin, Heidelberg. pp. 435–465.

Jiru M, Ranst E V. 2010. Increasing water productivity on Vertisols: Implications for environmental sustainability.Journal of the Science of Food and Agriculture, 90,2276–2281.

Ju X T, Xing G X, Chen X P, Zhang S L, Zhang L J, Liu X J, Cui Z L, Yin B, Christie P, Zhu Z L, Zhang F. 2009. Reducing environmental risk by improving N management in intensive Chinese agricultural systems. Proceedings of the National Academy of Sciences of the United States of America,106, 3041–3046.

Lv L, Yao Y, Zhang L, Dong Z, Jia X, Liang S, Ji J. 2013. Winter wheat grain yield and its components in the North China Plain: Irrigation management, cultivation, and climate.Chilean Journal of Agricultural Research, 73, 233–242.

Miao Y X, Stewart B A, Zhang F S. 2011. Long-term experiments for sustainable nutrient management in China. A review.Agronomy for Sustainable Development, 31, 397–414.

Moiwo J P, Yang Y H, Li H L, Han S M, Yang Y M. 2010. Impact of water resource exploitation on the hydrology and water storage in Baiyangdian Lake. Hydrological Processes, 24,3026–3039.

Mon J, Bronson K F, Hunsaker D J, Thorp K R, White J W,French A N. 2016. Interactive effects of nitrogen fertilization and irrigation on grain yield, canopy temperature, and nitrogen use efficiency in overhead sprinkler-irrigated durum wheat. Field Crops Research, 191, 54–65.

Qiu G Y, Wang L M, He X H, Zhang X Y, Chen S Y, Chen J, Yang Y H. 2008. Water use efficiency and evapotranspiration of winter wheat and its response to irrigation regime in the North China Plain. Agricultural and Forest Meteorology,148, 1848–1859.

Sayre K D, Rajaram S, Fischer R A. 1997. Yield potential progress in short bread wheats in northwest Mexico. Crop Science, 37, 36–42.

Shao L W, Zhang X Y, Sun H Y, Chen S Y, Wang Y M. 2011.Yield and water use response of winter wheat to winter irrigation in the North China Plain. Journal of Soil and Water Conservation, 66, 104–113.

Shearman V J, Sylvester-Bradley R, Scott R K, Foulkes M J.2005. Physiological progress associated with wheat yield progress in the UK. Crop Science, 45, 175–185.

Shirazi S M, Zardari N H, Yusop Z, Ismail Z, Othman F. 2014.Performance of wheat crop under different irrigation regimes and nitrogen levels. A field experiment. Journal of Environmental Protection & Ecology, 15, 973–982.

Sui J, Wang J D, Gong S H, Xu D, Zhang Y Q. 2015. Effect of nitrogen and irrigation application on water movement and nitrogen transport for a wheat crop under drip irrigation in the North China Plain. Water, 7, 6651–6672.

Sun H, Zhang X, Wang E, Chen S, Shao L. 2015. Quantifying the impact of irrigation on groundwater reserve and crop production - A case study in the North China Plain.European Journal of Agronomy, 70, 48–56.

Sun H Y, Zhang X Y, Chen S Y, Shao L W, Wang Y Z, Dong B F.2011. Analysis of field water consumption, its pattern, impact and driving factors. Chinese Journal of Eco-Agriculture, 19,1032–1038. (in Chinese)

Wang S J, Tian X H, Liu T, Lu X C, You D H, Li S. 2014.Irrigation, straw, and nitrogen management bene fits wheat yield and soil properties in a dryland agro-ecosystem.Agronomy Journal, 106, 2193–2201.

Wang X, Shi Y, Guo Z, Y Zhang, Yu Z. 2015. Water use and soil nitrate nitrogen changes under supplemental irrigation with nitrogen application rate in wheat field. Field Crops Research, 183, 117–125.

White E M, Wilson F E A. 2006. Responses of grain yield,biomass and harvest index and their rates of genetic progress to nitrogen availability in ten winter wheat varieties.Irish Journal of Agricultural and Food Research, 45, 85–101.

Wu Q, Wang G, Lin W, Zhang F. 2012. Estimating groundwater recharge of Taihang Mountain Piedmont in Luancheng County, Hebei Province, China. Geological Science &Technology Information, 31, 99–105. (in Chinese)

Xue Q, Zhu Z, Musick J T, Stewart B A, Dusek D A. 2006.Physiological mechanisms contributing to the increased water-use efficiency in winter wheat under de ficit irrigation.Journal of Plant Physiology, 163, 154–164.

Xu J, Li Y, Sun J, Du L, Zhang Y, Yu Q, Liu X. 2013. Comparative physiological and proteomic response to abrupt low temperature stress between two winter wheat cultivars differing in low temperature tolerance. Plant Biology, 15,292–303.

Yang Y H, Watanabe M, Zhang X Y, Zhang J Q, Wang Q X,Hayashi S. 2006. Optimizing irrigation management for wheat to reduce groundwater depletion in the piedmont region of the Taihang Mountains in the North China Plain.Agricultural Water Management, 82, 25–44.

Zhang X, Chen S, Sun H, Pei D, Wang Y. 2008. Dry matter,harvest index, grain yield and water use efficiency as affected by water supply in winter wheat. Irrigation Science,27, 1–10.

Zhang X, Chen S, Sun H, Wang Y, Shao L. 2010. Water use efficiency and associated traits in winter wheat cultivars in the North China Plain. Agricultural Water Management,97, 1117–1125.

Zhang X Y. 2011. A review of agricultural water-saving research at Luancheng Agro-Ecosystem Experimental Station of Chinese Academy of Sciences over the last 30 years. Chinese Journal of Eco-Agriculture, 19, 987–996.(in Chinese)

Zhang X, Wang Y, Sun H, Chen S, Shao L. 2012. Optimizing the yield of winter wheat by regulating water consumption during vegetative and reproductive stages under limited water supply. Irrigation Science, 31, 1103–1112.

Zhang X Y, Pei D, Hu C S. 2003. Conserving groundwater for irrigation in the North China Plain. Irrigation Science, 21,159–166.

Zhang Y M, Zhang J B, Hu C S, Zhao B Z, Zhu A N, Zhang L Z, Song L N. 2011. Effect of fertilization and irrigation on wheat-maize yield and soil nitrate nitrogen leaching in high agricultural yield region in North China Plain. Chinese Journal of Eco-Agriculture, 19, 532–539. (in Chinese)