Yield and water use responses of winter wheat to irrigation and nitrogen application in the North China Plain

Journal of Integrative Agriculture 2018, 17(5): 1194–1206 Available online at www.sciencedirect.com

ScienceDirect

RESEARCH ARTICLE

Yield and water use responses of winter wheat to irrigation and nitrogen application in the North China Plain ZHANG Ming-ming1, 2, DONG Bao-di1, QIAO Yun-zhou1, SHI Chang-hai3, YANG Hong1, 2, WANG Ya-kai1, 2, LIU Meng-yu1 1

Key Laboratory of Agricultural Water Resources, Hebei Laboratory of Agricultural Water-Saving, Center for Agricultural Resources Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Shijiazhuang 050021, P.R.China 2 University of Chinese Academy of Sciences, Beijing 100049, P.R.China 3 Qingdao Agricultural University, Qingdao 266109, P.R.China

Abstract With increasing water shortage resources and extravagant nitrogen application, there is an urgent need to optimize irrigation regimes and nitrogen management for winter wheat (Triticum aestivum L.) in the North China Plain (NCP). A 4-year field experiment was conducted to evaluate the effect of three irrigation levels (W1, irrigation once at jointing stage; W2, irrigation once at jointing and once at heading stage; W3, irrigation once at jointing, once at heading, and once at filling stage; 60 mm each irrigation) and four N fertilizer rates (N0, 0; N1, 100 kg N ha−1; N2, 200 kg N ha−1; N3, 300 kg N ha−1) on wheat yield, water use efficiency, fertilizer agronomic efficiency, and economic benefits. The results showed that wheat yield under W2 condition was similar to that under W3, and greater than that under W1 at the same nitrogen level. Yield with the N1 treatment was higher than that with the N0 treatment, but not significantly different from that obtained with the N2 and N3 treatments. The W2N1 treatment resulted in the highest water use and fertilizer agronomic efficiencies. Compared with local traditional practice (W3N3), the net income and output-input ratio of W2N1 were greater by 12.3 and 19.5%, respectively. These findings suggest that two irrigation events of 60 mm each coupled with application of 100 kg N ha–1 is sufficient to provide a high wheat yield during drought growing seasons in the NCP. Keywords: winter wheat, irrigation regime, nitrogen application, grain yield, water use efficiency

1. Introduction

Received 14 September, 2017 Accepted 21 December, 2017 ZHANG Ming-ming, E-mail: [email protected]. ac.cn; Correspondence DONG Bao-di, Tel: +86-311-85825949, Fax: +86-311-85815093, E-mail: [email protected] © 2018 CAAS. Publishing services by Elsevier B.V. All rights reserved. doi: 10.1016/S2095-3119(17)61883-5

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

ZHANG Ming-ming et al. Journal of Integrative Agriculture 2018, 17(5): 1194–1206

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

1195

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). 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 benefits.

2. Materials and methods 2.1. Experimental site Field experiments were carried out during four growing seasons from 2012 to 2016 at the Luancheng AgroEcosystem 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). 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 profile. 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 In this study, the winter wheat cultivar Shixin 828 was used in a maize-wheat rotation system over four growing

1196

Precipitation, temperature, and sunshine hours in October were the monthly mean from sowing to October 31. Precipitation, temperature, and sunshine hours in June were the monthly mean from June 1 to harvested day. ET0 is the mean monthly reference evapotranspiration.

5 Jun.

1)

6 8

9 7

9 6

2 5

5 5

8 6

6

24

3

20

26 24 23

20 11

1 2

39

4

May

14 15

Apr.

2

21

20

5

6

7 8

7 6

6 8

5 3

4 4

3 3

4 4

3 9

16

9 9 7

12 9

0 3

25

1

29

Mar.

2

15

14

1 2 0 –1 9 4 5 10 Feb.

15

4

7 5 2 3 2 2

4

1

1

5

5 4

4 2

3 1

1 1

2 1

1 1 –1

1 0

–4

–1 –1

–2 –5

–3 0

1 0

0

5 Jan.

0

4 Dec.

0

6

1 4

3 3

6 5

7 3

1 2

2 2

2 2 4 6

3 13 14

6

12 14

4 29

10 0

4 4

7 1

23 Nov.

2012–2013 2013–2014 2014–2015 2015–2016 2012–2013 2013–2014 2014–2015 2015–2016 2012–2013 2013–2014 2014–2015 2015–2016 2012–2013 2013–2014 2014–2015 2015–2016

Oct.

Month

Precipitation (mm)

Temperature (°C)

Meteorological factor1)

Table 1 Monthly meterological data during the winter wheat growing seasons in 2012–2016

ET0 (mm d–1)

Sunshine hours (h)

ZHANG Ming-ming et al. Journal of Integrative Agriculture 2018, 17(5): 1194–1206

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.

2.3. Measurements 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. 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. Total water use or evapotranspiration (ET) during the winter wheat growing season was calculated using the following water

ZHANG Ming-ming et al. Journal of Integrative Agriculture 2018, 17(5): 1194–1206

1197

Table 2 Experimental design Treatment1) W1N0 W1N1 W1N2 W1N3 W2N0 W2N1 W2N2 W2N3 W3N0 W3N1 W3N2 W3N3 1)

Irrigation time

Irrigation stage

1 1 1 1 2 2 2 2 3 3 3 3

Jointing Jointing Jointing Jointing Jointing, heading Jointing, heading Jointing, heading Jointing, heading Jointing, heading, filling Jointing, heading, filling Jointing, heading, filling Jointing, heading, filling

Irrigation amount (mm) 60 60 60 60 120 120 120 120 180 180 180 180

Total fertilization amount (kg ha−1) N P 0 0 100 172 200 172 300 172 0 0 100 172 200 172 300 172 0 0 100 172 200 172 300 172

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.

balance formula modified from that described by Zhang et al. (2012): ET=I+P+ΔSWC (1) 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. Water use efficiency at the grain yield level (WUEy) was defined 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): WUEy=Yield/ET (2) WUEbm=BM/ET (3) 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): AE (kg kg–1)=(YF–Y0)/AF (4) 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 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). 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 also

1198

ZHANG Ming-ming et al. Journal of Integrative Agriculture 2018, 17(5): 1194–1206

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) Source of variation Year (Y) Nitrogen (N) Irrigation (I) N×I Y×I Y×N Y×N×I

Height 101.7** 69.1** 153.1* 3.0* 12.4** 2.8** NS

BM 59.5** 190.6** 156.5** NS 7.2** 2.6** NS

Yield 105.7** 224.4** 225.5** 2.9* 12.9** 5.5** NS

SN 93.6** 29.6** 8.1** NS NS 3.8** NS

GW 266.9** 70.6** 108.0** 9.4** 9.2** 8.0** NS

GN 67.8** NS 41.2** NS 3.8** NS NS

ET 36.8** 17.1** 235.8** NS NS 10.5** NS

WUEy 26.1** 40.2** 20.8** NS 4.7** 3.9** NS

WUEbm 67.2** 69.6** 20.8** NS 10.7** 7.3** NS

AE 52.6** 49.5** 26.8** NS 2.6* NS NS

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.

W1

80

a

c

75

c

c

a

c

d

70

ab

a bc

a

a

65 60

80 75

e

d

cd cd

d

a

65

1N W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

1 2N 2 W 2N 3 W 3N 0 W 3N 1 W 3N 2 W 3N 3

0

W

W

W

2N

3 W

2N

2 W

1N

1 W

1N

0

1N

1N

Treatment

80 75

f

70

g

ef

de

cd

abc

bcd cd

a

a

ab

ef

65 60

Plant height (cm)

D 85 80 75

e

d

cd bcd d

abc bcd cd abc

ab

a

a

70 65

Treatment

1N W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

W

1N

0

1N 1 W 1N 2 W 1N 3 W 2N 0 W 2N 1 W 2N 2 W 2N 3 W 3N 0 W 3N 1 W 3N 2 W 3N 3

60

W

Plant height (cm)

a

70

Treatment C 85

W

abc abc abc bcd ab

60

W

W

W3

B 85 Plant height (cm)

Plant height (cm)

A 85

W2

Treatment

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.

significant (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). 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

1199

ZHANG Ming-ming et al. Journal of Integrative Agriculture 2018, 17(5): 1194–1206

a

a

bc

a

a

a

W3

B 22 000 20 000 18 000 16 000 14 000 12 000 10 000 8 000 6 000 4 000 2 000 0

f

cd cd c

e

Treatment

b

b

a

de

ab ab a

a

a

b

a

D 22 000 20 000 18 000 16 000 14 000 12 000 10 000 8 000 6 000 4 000 2 000 0

f

e

f

a

abc ab

abc bcd abc

de cde

f

W 1N W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

b

a

a

Biomass (kg ha–1)

c

b

ab ab

Treatment

W 1N W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

Biomass (kg ha–1)

C 22 000 20 000 18 000 16 000 14 000 12 000 10 000 8 000 6 000 4 000 2 000 0

b

W 1N W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

a

c

W 1N W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

Biomass (kg ha–1)

d

bc bc b

W2 Biomass (kg ha–1)

W1 A 22 000 20 000 18 000 16 000 14 000 12 000 10 000 8 000 6 000 4 000 2 000 0

Treatment

Treatment

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.

W1

6 000 c

b

b

b

b

a

a

a

b

a

a

a

4 000 2 000

Yield (kg ha–1)

Yield (kg ha–1)

10 000 8 000

W2

W3

B 12 000

A 12 000

0

10 000 8 000 6 000

d

b

b

b

a c

a

cd

2 000

1N W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

W

1N W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

W

Treatment

b

a

a

a

a

a

b

4 000 2 000

W

1N W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

0

Treatment

10 000 8 000 d 6 000

c

c

a

bc

a

ab a

a

a

d

d

4 000 2 000 0

1N W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

c

bc bc

a

W

10 000

D 12 000 Yield (kg ha–1)

Yield (kg ha–1)

a

4 000

Treatment

6 000 d

a

a

0

C 12 000 8 000

a

Treatment

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.

1200

ZHANG Ming-ming et al. Journal of Integrative Agriculture 2018, 17(5): 1194–1206

Table 4 Yield component factors of winter wheat as affected by irrigation and nitrogen during four growing seasons Treatment1) W1N0 W1N1 W1N2 W1N3 W2N0 W2N1 W2N2 W2N3 W3N0 W3N1 W3N2 W3N3

2012–20132) SN GW GN (no. (×104 ha−1) (g) spike–1) 547.7 d 28.9 f 31.8 d 621.5 bc 29.0 f 34.4 abc 643.4 ab 30.2 def 32.1 cd 638.3 ab 29.8 ef 32.1 cd 575.8 cd 31.9 cde 34.9 ab 628.6 abc 34.0 bc 35.4 a 685.1 a 33.7 bc 35.1 a 647.9 ab 33.3 bc 35.2 a 555.1 d 32.4 cd 34.9 ab 615.7 bc 36.9 a 34.3 abc 604.1 bcd 35.4 ab 35.1 a 642.6 ab 35.5 ab 32.6 bcd

2013–20142) 2014–20152) 2015–20162) SN GW GN (no. SN GW GN (no. SN GW GN (no. (×104 ha−1) (g) spike–1) (×104 ha−1) (g) spike–1) (×104 ha−1) (g) spike–1) 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 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 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 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 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 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 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 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 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 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 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 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

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.

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).

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. 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. 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 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.

1201

ZHANG Ming-ming et al. Journal of Integrative Agriculture 2018, 17(5): 1194–1206

ef

abc ab bc bc cd cd

W3

a

f

ef

f

f

bcd bcdecde de ab ab abc a

W

W

1N W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

ef

de

W2 B 400 350 300 250 200 150 100 50 0

ET (mm)

f

e

W1

1N W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

ET (mm)

A 400 350 300 250 200 150 100 50 0

Treatment

cd

1N W

Treatment

e

d

d

d

d

c

bc

a

a

ab

d

W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

cd de

bc

1N

ef

bc

D 400 350 300 250 200 150 100 50 0

W

ef

c

ab

ET (mm)

f

cd

a

W 0 1N W 1 1N W 2 1N W 3 2N W 0 2N W 1 2N W 2 2N W 3 3N W 0 3N W 1 3N W 2 3N 3

ET (mm)

C 400 350 300 250 ef 200 150 100 50 0

Treatment

Treatment

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 Treatment1) W1N0 W1N1 W1N2 W1N3 W2N0 W2N1 W2N2 W2N3 W3N0 W3N1 W3N2 W3N3

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

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.

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 The economic benefit of different treatments was analyzed using the average grain yield of the four years (Table 7).

1202

ZHANG Ming-ming et al. Journal of Integrative Agriculture 2018, 17(5): 1194–1206

W1

y=0.363x+6.106 R2=0.627 45

WUEy (kg ha–1 mm–1)

C 38 36 34 32 30 28 26 24 22 20 40

B

50 55 60 WUEbm (kg ha–1 mm–1)

65

70

y=0.562x–1.898 R2=0.926 45

50 55 60 WUEbm (kg ha–1 mm–1)

65

WUEy (kg ha–1 mm–1)

38 36 34 32 30 28 26 24 22 20 40

W2

70

W3 38 36 34 32 30 28 26 24 22 20 40

D 38 36 34 32 30 28 26 24 22 20 40

y=0.689x–7.678 R2=0.931 45

50 55 60 WUEbm (kg ha–1 mm–1)

WUEy (kg ha–1 mm–1)

WUEy (kg ha–1 mm–1)

A

65

70

y=0.497x+1.215 R2=0.913 45

50 55 60 WUEbm (kg ha–1 mm–1)

65

70

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 Treatment1) W1N1 W2N1 W3N1 W1N2 W2N2 W3N2 W1N3 W2N3 W3N3

AE (kg ha–1) 2012–2013 2013–2014 2014–2015 2015–2016 4.0 abc 5.8 cde 4.0 abc 6.7 cd 4.2 ab 9.1 ab 5.3 ab 9.9 a 5.6 a 9.3 a 5.6 ab 9.6 ab 3.2 bcd 4.3 ef 3.6 abc 5.0 de 4.1 abc 6.6 cd 3.8 abc 6.5 cd 3.3 bcd 7.3 bc 5.8 a 7.9 bc 2.2 d 3.6 f 3.4 bc 4.0 e 2.3 d 5.0 def 3.1 c 5.2 de 2.5 cd 5.7 cde 4.2 abc 6.1 cd

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.

Results indicated that within the different N treatments, W2 groups showed the greatest benefits. 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, fulfills 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

1203

ZHANG Ming-ming et al. Journal of Integrative Agriculture 2018, 17(5): 1194–1206

Table 7 The output-input ratio of different irrigation and nitrogen treatments Treatment1) W1N0 W1N1 W1N2 W1N3 W2N0 W2N1 W2N2 W2N3 W3N0 W3N1 W3N2 W3N3

Grain yield Output (kg ha–1)2) (CNY ha–1) 5 649.3 7 043.5 7 148.0 7 216.2 6 667.2 8 601.7 8 622.2 8 503.5 6 590.5 8 632.6 8 852.0 8 770.9

12 428.5 15 495.8 15 725.6 15 875.6 14 667.7 18 923.8 18 968.7 18 707.8 14 499.0 18 991.8 19 474.5 19 295.9

Input (CNY ha–1) Electricity for Seed Fertilizer Pesticide Machine Labor irrigation 600 0.0 450 300 3 000 900 600 1 541.4 450 300 3 000 900 600 2 019.3 450 300 3 000 900 600 2 497.8 450 300 3 000 900 600 0.0 450 600 3 000 1 275 600 1 541.4 450 600 3 000 1 275 600 2 019.3 450 600 3 000 1 275 600 2 497.8 450 600 3 000 1 275 600 0.0 450 900 3 000 1 650 600 1 541.4 450 900 3 000 1 650 600 2 019.3 450 900 3 000 1 650 600 2 497.8 450 900 3 000 1 650

Total input (CNY ha–1)

Net income (CNY ha–1)

5 250.0 6 791.4 7 269.3 7 747.8 5 925.0 7 466.4 7 944.3 8 422.8 6 600.0 8 141.4 8 619.3 9 097.8

7 178.5 8 704.4 8 456.3 8 127.8 8 742.7 11 457.4 11 024.4 10 285.0 7 899.0 10 850.4 10 855.2 10 198.1

Outputinput ratio 2.4 2.3 2.2 2.0 2.5 2.5 2.4 2.2 2.2 2.3 2.3 2.1

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.

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 benefits in terms of economizing the use of irrigation water and nitrogen fertilizer. 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 soil

1204

ZHANG Ming-ming et al. Journal of Integrative Agriculture 2018, 17(5): 1194–1206

moisture (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. An important goal of agricultural practices is to achieve continuous and efficient development of farmland, which takes into consideration both production and environmental benefits. 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

ZHANG Ming-ming et al. Journal of Integrative Agriculture 2018, 17(5): 1194–1206

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.

1205

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 benefits 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 deficit 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

1206

ZHANG Ming-ming et al. Journal of Integrative Agriculture 2018, 17(5): 1194–1206

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)

Sect ion editor ZHANG Wei-li Managing editor SUN Lu-juan