Nitrogen Fertilizer and Sowing Density Affect Flag Leaf Photosynthetic Characteristics, Grain Yield, and Yield Components of Oat in a Semiarid Region of Northwest China

Abstract

Oat has been gaining renewed interest due to its role in a healthy human diet, in animal feed, and as a source of high value compounds with industrial applications. Nitrogen fertilization and planting density are two of the most important crop management practices that affect the formation of yield components and final yield of oat. A 2 year 3 × 5 factorial field experiment was conducted to investigate the effects of nitrogen (N) fertilizer and planting density on the flag leaf photosynthetic characteristics, grain yield, and yield components under rainfed conditions. The experiment consisted of three sowing densities (60, 180, and 300 kg·ha⁻¹) and five nitrogen fertilizer rates (0, 45, 90, 135, and 180 kg·ha⁻¹). Results showed that the grain yield was significantly (p < 0.05) correlated with the leaf net photosynthesis rate (Pn), stomatal conductance (Gs), intercellular CO₂ concentration (Ci), transpiration rate (Tr), water-use efficiency (WUE), stomatal limitation value (Ls), chlorophyll content (SPAD value), leaf area index (LAI), panicle length, number of spikelets per panicle (NSP), number of grains per panicle (NGP), weight of grains per panicle (WGP), and 1000-kernel weight. Among the yield components, grain yield was driven by number of spikelets per panicle (NSP) and number of grains per panicle (NGP) at low (or high) planting density with low N supply, whereas, at high N supply, 1000-kernel weight was also an important factor for yield. Nitrogen fertilizer and sowing density had significant (p < 0.05) effects on the flag leaf photosynthetic characteristics, grain yield, and yield components of oat. The yield components increased and then decreased with the increase in nitrogen fertilizer, while they decreased with the increase in planting density. The maximum values (p < 0.05) of grain yield were observed in the nitrogen fertilization of 90 kg·ha⁻¹ and sowing density of 180 kg·ha⁻¹ treatment in both growing seasons, mainly contributing to the improved leaf photosynthesis traits (Pn, Gs, Tr, Ls, SPAD, and LAI). The combination of nitrogen fertilization of 90 kg·ha⁻¹ and sowing density of 180 kg·ha⁻¹ is suitable for oat production on a cool semiarid plateau or other agroecozones with similar environmental conditions.

1. Introduction

Oat (Avena sativa) is distributed in more than 40 countries around the world, mainly concentrated in the northern regions of 40° N latitude including Asia, Europe, and North America, with a total output of more than 430 million tons [1]. In China, oat is extensively planted in northern and western regions. Traditionally the primary regions of cultivation include Inner Mongolia, Ganshu, Hebei, Shanxi, and Qinghai provinces, which account for 85% of the total oat production area [2]. Oat is a low-input cereal used for human food consumption, as well as for animal feed and as a source of high value compounds with industrial applications [3]. In recent decades, discussion on oat grain dietetic value and suitability for the production of functional food has increased in the public and scientific literature [4,5]. With increased attention on consuming healthy grains and a growing understanding of oat, the planting area and range of oat in China are gradually expanding, especially in Qinghai province.

Qinghai is located in the west of China with a plateau continental climate and is one of the important provinces on the Qinghai–Tibet Plateau. Oat, a gramineous herbaceous plant, is suitable for growing in cold areas with long days of sunshine, short frost-free periods, and low temperatures. There is, however, a lack of knowledge about sustainable production strategies suitable for local production conditions.

The use of nitrogen fertilizer is a major factor in the profitable production of most crops in temperate environments, because it affects dry matter production by influencing leaf area development and maintenance, as well as photosynthetic efficiency [6,7]. Amounts of fertilizer nitrogen applied to most cereals have increased in recent decades, and progress has been made in improving the prediction of optimum nitrogen requirement [6]. However, the benefits of nitrogen fertilizer for oat are less clear-cut than for other cereal species. The relatively taller stature of oat compared to other cereals makes it susceptible to lodging. The sensitivity of oat cultivars to lodging is an important factor which limits the application and effectiveness of nitrogen fertilization in maximizing grain yield [1]. Nitrogen fertilizer can affect both vegetative and reproductive development of crops, while nitrogen deficiency delays both vegetative and reproductive phonological development, as well as reduces leaf emergence rate, yield, and yield components [7]. Peltonen-Sainio and Järvinen [8] found that additional nitrogen applied before anthesis increased floret set and survival and resulted in more grains per panicle compared with lower nitrogen treatments. Moreover, a closely positive correlation was found between the photosynthetic capacity of leaves and nitrogen fertilizer; the net photosynthesis rate (Pn), stomatal conductance (Gs), transpiration rate (Tr), stomatal limitation value (Ls), chlorophyll content (SPAD value), and leaf area index (LAI) increased and then decreased with the increase in nitrogen fertilizer [7]. The flag leaf stays longest on the plant and makes a major contribution to the grain yield in cereals [9]. Under favorable conditions, approximately 70–90% of the total grain yield is derived from the photosynthates accumulated during grain filling [10].

Planting density is an important crop management that affects the grain yield by regulating growth, photosynthesis, and yield components, which are the target traits closely related to crops [11]. Increases in plant density for crop production generally have a large positive impact on the canopy LAI and, as a consequence, on biomass and dry matter accumulation [12,13]. However, excess LAI often causes leaf shading and lower leaf nutrient concentration, reducing canopy photosynthesis [14], accelerating leaf senescence, and reducing grain yield [15], which means that high planting density results in strong competition. In potato and maize, the light-saturated rate of leaf photosynthesis can be described as a logistic function of the nitrogen concentration per unit leaf area [14]. An understanding of how yield accumulation is influenced by plant density and nitrogen application is necessary to guide farmers’ practice toward achieving high grain yield of oat.

Most previous studies with nitrogen fertilizer and planting density have focused on the N uptake, grain yield, yield components, grain quality, forage yield, and quality of soil for intercropped oat and various legumes [6,9,16,17,18,19,20,21]. However, the influence of nitrogen fertilizer and planting density on the flag leaf physiological characteristics of oat has not been fully studied. We hypothesized that climate plays a decisive role in the grain yield of oat in semiarid environments, while nitrogen fertilizer and planting density have an important regulating effect on the grain yield of oat. Thus, the objectives of the present study were to assess the influence of nitrogen fertilizer and planting density on the leaf photosynthetic characteristics, grain yield, and yield components of oat.

2. Materials and Methods

2.1. Study Site

A field experiment was conducted during the 2016 and 2017 growing seasons under rainfed conditions at the experimental field station of Qinghai Academy of Animal and Veterinary Science (101°37′49″ E, 36°30′14″ N, Huangzhong County, Xining City, Qinghai, China). The station is located 30 km southwest of Qinghai Academy of Animal and Veterinary Science at an elevation of 2670 m above sea level. The region is a typical semiarid area with 540 mm annual precipitation that mainly occurs during May to September. It has a temperate, semiarid, and continental climate, with an annual mean temperature of 3.7 °C and up to 2774 °C of effective accumulated temperature. The experiment was established in a chestnut soil with pH 8.32, organic matter 18.8 g·kg⁻¹ organic matter, 1.54 g·kg⁻¹ total N, 136 mg·kg⁻¹ available N, 0.73 g·kg⁻¹ total P, 35 mg·kg⁻¹ available P, 24.9 g·kg⁻¹ total K, and 127 mg·kg⁻¹ available K (0 to 30 cm depth). The previous crop was common buckwheat (Fagopyrum esculentum M.), which was harvested in mid-September of 2015. The same field was used in both years. The long-term average monthly precipitation and temperature (1986–2015) from April to September and the deviations during the two growing seasons of the experimentation are given in

Table 1. Long-term average monthly precipitation and temperature (1986–2015) and deviations during the 2016 and 2017 growing seasons of experimentation at Huangzhong, Qinghai Province.

Month	Precipitation (mm)			Temperature (°C)
		30 Year Average	Deviations		30 Year Average	Deviations
		2016	2017		2016	2017
April	32.1	−2.5	5.8	5.5	1.9	0.2
may	70.4	−14.9	41.2	9.9	−0.1	0.0
June	86.7	−32.5	−41.1	13.1	1.4	0.3
July	106.0	3.5	−57.0	15.2	1.4	2.8
August	100.2	−59.2	38.8	14.2	3.4	0.4
September	80.7	0.7	2.1	9.9	0.7	1.4

. Monthly precipitation was higher in July 2016 than 2017, but close to the long-term norm in May and August 2017, comparable to values in 2016. Overall, the experimental year 2016 was comparatively dry.

2.2. Experimental Design and Field Management

A 3 × 5 factorial experiment was arranged in a randomized complete block design with three replications. There were three different sowing densities: 60 kg·ha⁻¹ (D1), 180 kg·ha⁻¹ (D2), and 300 kg·ha⁻¹ (D3). For each sowing density, five different nitrogen fertilizer rates application were designated as 0 kg·ha⁻¹ (N0), 45 kg·ha⁻¹ (N1), 90 kg·ha⁻¹ (N2), 135 kg·ha⁻¹ (N3), and 180 kg·ha⁻¹(N4).

Qingyan No.1 is a new variety for grain and feed production that has been bred by the Qinghai Academy of Animal and Veterinary Science through hybrid breeding technology. The variety, with high grain yield and biomass, as well as good lodging resistance, was used in this study. Seedbed preparation included moldboard ploughing, disc harrowing, and cultivation. Phosphorus, as calcium superphosphate, was applied at 45 kg·ha⁻¹ P₂O₅ during land preparation. The nitrogen fertilizer at specific rates was applied to the target plots as basal fertilizer. Three seeding rates were used, with target plant densities of 150, 450, and 750 plants··m⁻². In this study, the intermediate density (450 plants·m⁻²) represents the normal density commonly used in farm fields. Each plot had an area of 20 m² (4 m × 5 m), consisted of 20 rows 5 m long, with the interrow distances of 0.25 m, while border rows were not included for any sampling. The neighboring plots were separated by a 1 m buffer zone. Similarly, the blocks were separated by a 2 m buffer zone. Sowing was conducted by hand and seeds were placed in every row at a depth of 2–3 cm. Sowing was performed on 23 April 2016 and 29 April 2017.

2.3. Leaf Photosynthetic Characteristics

At the beginning of the boot stage (2 July 2016 and 6 July 2017), five plants in each plot were selected randomly, and the flag leaf was labeled with white thread to investigate the gas exchange and chlorophyll content [22]. The gas exchange was measured using an LI-6400 photosynthesis system (Li-Cor Inc., Lincoln, NE, United States of America) between 9:00 and 11:00 a.m. Light intensity, temperature, CO₂ concentration, flow rate, and relative humidity were maintained at 1200 μmol·m⁻²·s⁻¹, 25 °C, 400 ± 5.0 μmol·mol⁻¹, 0.5 L·min⁻¹, and 30% ± 1.0% respectively. The Pn, Gs, Tr, ambient CO₂ concentration (Ca), and intercellular CO₂ concentration (Ci) were automatically recorded. Water use efficiency (WUE) was calculated as Pn/Tr [23]. Ls was calculated as 1 − Ci/Ca [7]. At the same time, five labeled leaves in each plot were chosen to measure the chlorophyll content (with a SPAD-502 Plus chlorophyll meter, Konica Minolta, Japan), according to the method of Abdelhamidg et al. (2003). For each leaf, measurements of SPAD values were made on five points from the leaf bottom to the tip and averaged. The total leaf area of each labeled plant was measured with a YMJ-A leaf area meter (Zhejiang Top Yunnong Technology Co., Ltd., Hangzhou, China). At that time, the number of plants from a randomly selected 1 m² area of each plot was counted for the calculation of LAI. The LAI was calculated as the total leaf area of one plant (m²·plant⁻¹) × plant density (plants·m⁻²) [24].

2.4. Yield and Yield Components Measurement

Grain yield at harvest was evaluated. Whole plot plants were hand-harvested on 25 August 2016 and 5 September 2017. Twenty plants in each plot were selected randomly to investigate the plant height, panicle length, number of spikelets per panicle (NSP), number of grains per panicle (NGP), and weight of grains per panicle (WGP) at harvest [25]. The plant height refers to the distance from the stem base to the top of spikes. The 1000-kernel weight was determined by measuring the weight of 200 kernels from each plot and multiplying by 5. Grain yield was determined by harvesting plants from 1 m² area in each plot. For grain yield and 1000-grain weight, seeds were air-dried for 2 weeks before measurement.

2.5. Data Analyses

The mean grain yield data from each treatment were plotted against leaf photosynthetic characteristics and yield components. Nitrogen fertilizer rates and planting densities were plotted against leaf photosynthetic characteristics, grain yield, and yield components. Linear (y = a + bx), quadratic (y = a + bx − cx²), hyperbolic (x²a⁻² − y²b⁻² = 1 or y²a⁻² − x²b⁻² = 1), and logarithmic (y = log_ax) equations were tested for their suitability to describe the relationship for grain yield response with leaf photosynthetic characteristics and yield components, as well as for nitrogen fertilizer rates and planting densities with leaf photosynthetic characteristics, grain yield, and yield components. The equation with the highest coefficient of determination (R²) value was judged to be the most appropriate. In these regression equations, grain yield, nitrogen fertilizer rates, or planting densities was the independent variable (x), leaf photosynthetic characteristics and yield components or leaf photosynthetic characteristics, grain yield, and yield components were the dependent variable (y). The R² comparisons among the models tested showed that the quadratic regression equations (y = a + bx − cx²) had the best fit for grain yield of oat over flag leaf photosynthetic characteristics and yield components, as well as nitrogen fertilizer or planting density over flag leaf photosynthetic characteristics, grain yield, and yield components. SPSS 17.0 programs (SPSS Institute Inc., Chicago, IL, USA) were used to conduct the analyses of variance (ANOVA) and the regression analysis. Treatment mean differences were separated by the least significant difference (LSD) test at the 0.05 probability level.

3. Results

3.1. Flag Leaf Photosynthetic Characteristics

3.1.1. Flag Leaf Gas Exchange

The Pn was significantly (p < 0.01) affected by nitrogen fertilizer, planting density, and their interaction during the two growing seasons (

Table 2. Net photosynthetic rate (Pn), stomatal conductance (Gs), intercellular CO₂ concentration (Ci), transpiration rate (Tr), water-use efficiency (WUE), stomatal limitation value (Ls), and leaf area index (LAI) of oat in different nitrogen fertilizer and planting density treatments in two cropping seasons. Data are expressed as the mean of three replications (n = 3). Data within a column in the same year, sharing the same letter, are not significantly different at p < 0.05; ns, not significant; ** significant at p < 0.01. N, D, and N × D represent nitrogen fertilizer, planting density, and the interaction between nitrogen fertilizer and planting density, respectively. N0, N1, N2, N3, and N4 represent nitrogen fertilizer rates at five levels of 0, 45, 90, 135, and 180 kg·ha⁻¹, respectively. D1, D2, and D3 represent planting density at three sowing densities of 60, 180, and 300 kg·ha⁻¹, respectively.

Years	Treatments	Pn (μmol·m−2·s−1)	Gs (mol·m−2·s−1)	Ci (μmol·mol−1)	Tr (mmol·m−2·s−1)	WUE (%)	LS	SPAD	LAI
2016	N0D1	9.95 ± 0.66 fg	0.23 ± 0.02 ghi	343.83 ± 4.57 bc	2.83 ± 0.24 fg	0.35 ± 0.01 cd	0.33 ± 0.01 g	40.17 ± 1.94 f	1.75 ± 0.11 i
	N1D1	11.36 ± 0.56 e	0.25 ± 0.03 fg	336.27 ± 5.22 cde	3.23 ± 0.24 def	0.35 ± 0.01 cd	0.36 ± 0.01 f	41.20 ± 1.51 f	2.01 ± 0.17 h
	N2D1	12.97 ± 1.08 d	0.31 ± 0.03 e	320.83 ± 3.52 ghi	4.27 ± 0.31 c	0.30 ± 0.01 cde	0.41 ± 0.01 e	46.03 ± 1.47 e	3.21 ± 0.17 fg
	N3D1	9.55 ± 0.72 fg	0.22 ± 0.02 ghi	344.46 ± 2.81 bc	2.29 ± 0.35 ghi	0.42 ± 0.04 b	0.31 ± 0.01 gh	35.97 ± 1.55 g	3.30 ± 0.13 f
	N4D1	8.94 ± 0.29 g	0.19 ± 0.01 hi	354.71 ± 3.43 a	1.84 ± 0.33 i	0.50 ± 0.09 a	0.25 ± 0.01 i	33.10 ± 0.44 g	3.54 ± 0.08 e
	N0D2	11.41 ± 0.56 e	0.34 ± 0.03 de	328.52 ± 3.17 efg	3.45 ± 0.24 de	0.33 ± 0.02 cd	0.42 ± 0.01 e	45.20 ± 2.67 e	3.02 ± 0.13 g
	N1D2	14.32 ± 0.21 bc	0.49 ± 0.03 b	321.55 ± 7.24 ghi	4.75 ± 0.09 bc	0.30 ± 0.01 cde	0.50 ± 0.02 c	57.63 ± 0.64 c	3.36 ± 0.11 ef
	N2D2	16.76 ± 0.48 a	0.63 ± 0.04 a	295.99 ± 6.51 j	6.61 ± 0.52 a	0.25 ± 0.01 e	0.62 ± 0.01 a	70.90 ± 2.01 a	3.57 ± 0.15 e
	N3D2	14.65 ± 0.76 b	0.51 ± 0.02 b	317.88 ± 5.15 hi	4.92 ± 0.56 b	0.30 ± 0.02 de	0.55 ± 0.02 b	64.80 ± 2.36 b	3.99 ± 0.16 d
	N4D2	13.58 ± 0.22 cd	0.38 ± 0.03 cd	331.75 ± 6.41 def	4.19 ± 0.10 c	0.32 ± 0.00 cd	0.51 ± 0.01 c	53.50 ± 2.63 d	4.77 ± 0.08 b
	N0D3	8.84 ± 0.28 g	0.18 ± 0.03 i	351.39 ± 6.18 ab	2.11 ± 0.19 hi	0.42 ± 0.03 b	0.24 ± 0.01 i	40.67 ± 1.99 f	3.26 ± 0.09 f
	N1D3	9.43 ± 1.10 fg	0.24 ± 0.02 gh	338.65 ± 3.50 cd	2.66 ± 0.48 fgh	0.36 ± 0.02 c	0.30 ± 0.02 h	40.67 ± 0.76 f	3.84 ± 0.18 d
	N2D3	13.26 ± 0.37 cd	0.41 ± 0.04 c	316.29 ± 5.39 i	4.34 ± 0.11 c	0.31 ± 0.00 cde	0.46 ± 0.01 d	52.60 ± 2.01 d	4.28 ± 0.15 c
	N3D3	11.22 ± 0.33 e	0.30 ± 0.03 e	326.16 ± 6.91 fgh	3.56 ± 0.30 d	0.32 ± 0.02 cd	0.40 ± 0.01 e	46.63 ± 1.95 e	4.77 ± 0.05 b
	N4D3	10.09 ± 0.41 f	0.29 ± 0.03 ef	332.76 ± 5.36 def	2.97 ± 0.22 ef	0.34 ± 0.02 cd	0.37 ± 0.01 f	41.90 ± 1.23 f	5.35 ± 0.10 a
LSD0.05 for
N		**	**	**	**	**	**	**	**
D		**	**	**	**	**	**	**	**
N × D		**	**	**	**	**	**	**	**
2017	N0D1	11.82 ± 0.16 e	0.23 ± 0.02 fg	336.71 ± 5.67 ef	3.50 ± 0.14 de	0.34 ± 0.01 c	0.38 ± 0.03 ef	42.73 ± 1.32 fgh	2.27 ± 0.05 j
	N1D1	12.08 ± 0.16 e	0.27 ± 0.03 ef	336.23 ± 3.76 ef	3.24 ± 0.35 e	0.38 ± 0.04 bc	0.39 ± 0.01 de	45.23 ± 2.21 ef	2.80 ± 0.21i
	N2D1	13.26 ± 0.15 d	0.34 ± 0.04 d	320.86 ± 2.09 g	3.73 ± 0.11 cd	0.36 ± 0.01 bc	0.41 ± 0.01 cd	51.90 ± 1.55 c	3.17 ± 0.08 h
	N3D1	9.33 ± 0.97 h	0.22 ± 0.02 gh	347.76 ± 1.15 c	2.35 ± 0.06 g	0.40 ± 0.04 b	0.34 ± 0.04 g	37.13 ± 1.88 ij	3.60 ± 0.05 g
	N4D1	8.71 ± 0.11 i	0.18 ± 0.02 hi	358.86 ± 2.25 a	1.83 ± 0.16 h	0.48 ± 0.05 a	0.26 ± 0.01 hi	34.67 ± 2.74 j	3.96 ± 0.10 ef
	N0D2	10.83 ± 0.1 f3	0.32 ± 0.02 d	341.41 ± 1.83 de	3.30 ± 0.13 e	0.33 ± 0.01 c	0.40 ± 0.01 de	43.83 ± 1.97 fg	3.92 ± 0.13 f
	N1D2	13.90 ± 0.16 c	0.43 ± 0.03 c	333.71 ± 4.50 f	4.03 ± 0.19 c	0.35 ± 0.01 c	0.44 ± 0.03 c	54.13 ± 1.29 bc	4.42 ± 0.13 d
	N2D2	15.66 ± 0.24 a	0.58 ± 0.03 a	310.56 ± 5.71 h	5.69 ± 0.07 a	0.28 ± 0.01 d	0.57 ± 0.02 a	62.73 ± 1.19 a	4.93 ± 0.08 c
	N3D2	14.97 ± 0.19 b	0.52 ± 0.02 b	324.65 ± 3.65 g	4.44 ± 0.36 b	0.34 ± 0.02 c	0.49 ± 0.01 b	56.77 ± 1.17 b	5.29 ± 0.26 b
	N4D2	13.31 ± 0.21 d	0.41 ± 0.01 c	340.58 ± 3.78 de	3.80 ± 0.11 cd	0.35 ± 0.01 bc	0.43 ± 0.03 cd	52.17 ± 1.08 c	6.00 ± 0.22 a
	N0D3	9.05 ± 0.27 hi	0.17 ± 0.01 i	357.20 ± 1.62 ab	2.06 ± 0.17 gh	0.44 ± 0.03 a	0.23 ± 0.01 i	39.83 ± 1.53 hi	3.48 ± 0.05 g
	N1D3	9.34 ± 0.45 h	0.18 ± 0.02 hi	351.71 ± 4.52 bc	1.95 ± 0.10 h	0.48 ± 0.01 a	0.28 ± 0.01 h	40.93 ± 1.07 gh	3.87 ± 0.16 f
	N2D3	12.75 ± 0.22 d	0.36 ± 0.02 d	336.11 ± 2.40 ef	3.45 ± 0.06 de	0.37 ± 0.01 bc	0.41 ± 0.01 cd	50.77 ± 0.51 cd	4.16 ± 0.05 e
	N3D3	10.67 ± 0.19 f	0.28 ± 0.03 e	341.30 ± 1.31 de	3.21 ± 0.10 e	0.33 ± 0.00 c	0.38 ± 0.00 ef	47.63 ± 1.60 de	4.65 ± 0.10 d
	N4D3	9.95 ± 0.27 g	0.25 ± 0.03 efg	346.17 ± 0.72 cd	2.68 ± 0.36 f	0.38 ± 0.04 bc	0.35 ± 0.01 fg	42.03 ± 0.91 fgh	5.41 ± 0.21 b
LSD0.05 for
N		**	**	**	**	**	**	**	**
D		**	**	**	**	**	**	**	**
N × D		**	**	**	**	**	**	**	ns

). The Pn increased from 10.32 to 14.11 μmol·m⁻²·s⁻¹ (mean of two growing seasons and three sowing densities) as the nitrogen fertilizer rate increased from 0 to 90 kg·ha⁻¹, and then decreased to 10.77 μmol·m⁻²·s⁻¹ (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 180 kg·ha⁻¹ resulted in an increase in Pn from 10.80 to 13.94 μmol·m⁻²·s⁻¹ (mean of two growing seasons and five nitrogen fertilizer rates), but it was decreased to 10.46 μmol·m⁻²·s⁻¹ (mean of two growing seasons and five nitrogen fertilizer rates) when increasing sowing density to 300 kg·ha⁻¹. The highest Pn occurred in the combination of 90 kg·ha⁻¹ N and 180 kg·ha⁻¹ sowing density (N2D2) in each growing season.

Nitrogen fertilizer, planting density, and their interaction had significant effects (p < 0.01) on the Gs during the two growing seasons (Table 2). The Gs increased from 0.24 to 0.44 mol·m⁻²·s⁻¹ (mean of two growing seasons and three sowing densities) as the nitrogen fertilizer rate increased from 0 to 90 kg·ha⁻¹, and then decreased to 0.28 mol·m⁻²·s⁻¹ (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 180 kg·ha⁻¹ resulted in an increase in Gs from 0.25 to 0.46 mol·m⁻²·s⁻¹ (mean of two growing seasons and five nitrogen fertilizer rates), but it was decreased to 0.27 mol·m⁻²·s⁻¹ (mean of two growing seasons and five nitrogen fertilizer rates) when increasing sowing density to 300 kg·ha⁻¹. The highest Gs occurred in the combination of 90 kg·ha⁻¹ N and 180 kg·ha⁻¹ sowing density (N2D2) in each growing season.

Nitrogen fertilizer, planting density, and their interaction had significant effects (p < 0.01) on the Ci during the two growing seasons (Table 2). Increasing the nitrogen fertilizer rate from 0 to 90 kg·ha⁻¹ resulted in a decrease in the Ci from 343.18 to 316.77 μmol·mol⁻¹ (mean of two growing seasons and three sowing densities), and then an increase to 344.14 μmol·mol⁻¹ (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 180 kg·ha⁻¹ resulted in a decrease in Ci from 340.05 to 324.66 μmol·mol⁻¹ (mean of two growing seasons and five nitrogen fertilizer rates), but it was increased to 339.77 μmol·mol⁻¹ (mean of two growing seasons and five nitrogen fertilizer rates) when increasing sowing density to 300 kg·ha⁻¹. The lowest Ci occurred in the combination of 90 kg·ha⁻¹ N and 180 kg·ha⁻¹ sowing density (N2D2) in each growing season.

Nitrogen fertilizer, planting density, and their interaction had significant effects (p < 0.01) on the Tr during the two growing seasons (Table 2). The Tr increased from 2.87 to 4.68 mmol·m⁻²·s⁻¹ (mean of two growing seasons and three sowing densities) as the nitrogen fertilizer rate increased from 0 to 90 kg·ha⁻¹, and then decreased to 2.88 mmol·m⁻²·s⁻¹ (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 180 kg·ha⁻¹ resulted in an increase in Tr from 2.91 to 4.52 mmol·m⁻²·s⁻¹ (mean of two growing seasons and five nitrogen fertilizer rates), but it was decreased to 2.90 mmol·m⁻²·s⁻¹ (mean of two growing seasons and five nitrogen fertilizer rates) when increasing sowing density to 300 kg·ha⁻¹. The highest Tr occurred in the combination of 90 kg·ha⁻¹ N and 180 kg·ha⁻¹ sowing density (N2D2) in each growing season.

The WUE was significantly (p < 0.01) affected by nitrogen fertilizer, planting density, and their interaction during the two growing seasons (Table 2). Increasing the nitrogen fertilizer rate from 0 to 90 kg·ha⁻¹ resulted in a decrease in the WUE from 0.37% to 0.31% (mean of two growing seasons and three sowing densities), and then increased to 0.39% (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 180 kg·ha⁻¹ resulted in a decrease in WUE from 0.39% to 0.31% (mean of two growing seasons and five nitrogen fertilizer rates), but it was increased to 0.37% (mean of two growing seasons and five nitrogen fertilizer rates) when increasing sowing density to 300 kg·ha⁻¹. The lowest WUE occurred in the combination of 90 kg·ha⁻¹ N and 180 kg·ha⁻¹ sowing density (N2D2) in each growing season.

The Ls was significantly (p < 0.01) affected by nitrogen fertilizer, planting density, and their interaction during the two growing seasons (Table 2). The Ls increased from 0.33 to 0.48 (mean of two growing seasons and three sowing densities) as the nitrogen fertilizer rate increased from 0 to 90 kg·ha⁻¹, and then decreased to 0.36 (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 180 kg·ha⁻¹ resulted in an increase in Ls from 0.34 to 0.49 (mean of two growing seasons and five nitrogen fertilizer rates), but it was decreased to 0.34 (mean of two growing seasons and five nitrogen fertilizer rates) when increasing sowing density to 300 kg·ha⁻¹. The highest Ls occurred in the combination of 90 kg·ha⁻¹ N and 180 kg·ha⁻¹ sowing density (N2D2) in each growing season.

3.1.2. Flag Leaf Chlorophyll Content

The flag leaf chlorophyll content (SPAD value) was significantly (p < 0.01) affected by nitrogen fertilizer, planting density, and their interaction during the two growing seasons (Table 2). The SPAD increased from 42.07 to 55.83 (mean of two growing seasons and three sowing densities) as the nitrogen fertilizer rate increased from 0 to 90 kg·ha⁻¹, and then decreased to 42.88 (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 180 kg·ha⁻¹ resulted in an increase in SPAD from 40.81 to 56.15 (mean of two growing seasons and five nitrogen fertilizer rates), but it was decreased to 44.37 (mean of two growing seasons and five nitrogen fertilizer rates) when increasing sowing density to 300 kg·ha⁻¹. The highest SPAD occurred in the combination of 90 kg·ha⁻¹ N and 180 kg·ha⁻¹ sowing density (N2D2) in each growing season.

3.1.3. Leaf Area Index

Nitrogen fertilizer and planting density had significant effects (p < 0.01) on the leaf area index (LAI) during the two growing seasons. Interaction between nitrogen fertilizer and planting density had a significant effect (p < 0.01) on LAI in 2016, but had no significant effect (p > 0.05) in 2017 (Table 2). The LAI increased from 2.95 to 3.89 (mean of two growing seasons and three sowing densities) as the nitrogen fertilizer rate increased from 0 to 90 kg·ha⁻¹, and then increased to 4.84 (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 180 kg·ha⁻¹ resulted in an increase in LAI from 2.96 to 4.33 (mean of two growing seasons and five nitrogen fertilizer rates), but it was decreased to 4.31 (mean of two growing seasons and five nitrogen fertilizer rates) when increasing sowing density to 300 kg·ha⁻¹. The highest LAI occurred in the combination of 180 kg·ha⁻¹ N and 300 kg·ha⁻¹ sowing density (N4D3) in 2016, and the combination of 180 kg·ha⁻¹ N and 180 kg·ha⁻¹ sowing density (N4D2) in 2017.

3.2. Grain Yield and Yield Components

3.2.1. Grain Yield

The grain yield was significantly (p < 0.01) affected by nitrogen fertilizer, planting density, and their interaction during the two growing seasons (

Table 3. Grain yield, plant height, panicle length, number of spikelets per panicle (NSP), number of grains per panicle (NGP), weight of grains per panicle (WGP), and 1000-kernel weight of oat in different nitrogen fertilizer and planting density treatments in two cropping seasons. Data are expressed as the mean of three replications (n = 3). Data within a column in the same year, sharing the same letter, are not significantly different at p < 0.05; ns, not significant; * and ** significant at p < 0.05 and p < 0.01, respectively. N, D, and N × D represent nitrogen fertilizer, planting density, and the interaction between nitrogen fertilizer and planting density, respectively. N0, N1, N2, N3, and N4 represent nitrogen fertilizer rates at five levels of 0, 45, 90, 135, and 180 kg·ha⁻¹, respectively. D1, D2, and D3 represent planting density at three sowing densities of 60, 180, and 300 kg·ha⁻¹, respectively.

Years	Treatments	Grain Yield (kg·ha−1)	Plant Height (cm)	Panicle Length (cm)	NSP	NGP	WGP (g)	1000-Kernel Weight (g)
2016	N0D1	2500.93 ± 60.63 j	97.94 ± 2.34 ef	17.77 ± 0.49 i	76.12 ± 1.54 c	174.73 ± 3.93 c	5.09 ± 0.11 d	32.38 ± 0.07 g
	N1D1	2668.60 ± 33.95 i	104.37 ± 0.93 cd	22.43 ± 0.22 e	82.98 ± 1.62 b	191.59 ± 2.36 b	5.50 ± 0.14 c	36.61 ± 0.27 c
	N2D1	2957.30 ± 39.93 g	108.88 ± 3.54 b	27.22 ± 0.91 a	88.14 ± 1.51 a	206.39 ± 4.53 a	6.18 ± 0.10 a	38.21 ± 0.53 a
	N3D1	2913.10 ± 51.42 gh	115.09 ± 2.38 a	24.95 ± 0.39 c	68.12 ± 0.65 de	140.38 ± 5.54 ef	5.32 ± 0.05 cd	34.59 ± 0.33 e
	N4D1	2846.17 ± 38.42 h	107.71 ± 1.77 bc	19.83 ± 0.29 g	62.33 ± 2.91 fg	100.32 ± 7.52 h	4.66 ± 0.30 e	30.48 ± 0.21 i
	N0D2	2919.00 ± 81.80 gh	87.66 ± 2.68 g	15.60 ± 0.42 j	63.57 ± 3.84 fg	138.08 ± 5.31 f	4.03 ± 0.18 h	30.86 ± 0.29 i
	N1D2	3423.97 ± 45.16 e	89.65 ± 0.90 g	21.17 ± 0.17 f	64.93 ± 0.53 ef	147.39 ± 3.93 de	4.32 ± 0.05 fg	35.59 ± 0.17 d
	N2D2	4002.00 ± 47.31 a	100.42 ± 2.86 de	26.08 ± 0.28 b	69.02 ± 0.47 d	150.91 ± 3.45 d	5.77 ± 0.01 b	37.23 ± 0.07 b
	N3D2	3868.63 ± 35.17 b	103.72 ± 2.28 cd	23.53 ± 0.27 d	59.94 ± 1.95 g	104.94 ± 1.12 h	4.42 ± 0.06 ef	33.26 ± 0.08 f
	N4D2	3757.43 ± 51.01 c	95.33 ± 3.04 f	18.67 ± 0.73 h	48.31 ± 1.16 j	63.56 ± 2.07 j	3.31 ± 0.04 i	29.55 ± 0.15 j
	N0D3	2623.53 ± 79.72 i	76.33 ± 2.08 i	13.01 ± 0.43 k	49.78 ± 4.67 ij	91.49 ± 5.09 i	2.55 ± 0.10 k	28.74 ± 0.07 k
	N1D3	3153.97 ± 28.82 f	83.48 ± 0.98 h	20.23 ± 0.22 g	52.77 ± 1.57 hi	105.56 ± 2.37 h	2.85 ± 0.33 j	32.43 ± 0.61 g
	N2D3	3624.37 ± 51.05 d	89.22 ± 1.23 g	25.03 ± 0.17 c	55.59 ± 1.04 h	118.11 ± 1.25 g	4.15 ± 0.03 gh	35.49 ± 0.34 d
	N3D3	3402.00 ± 66.55 e	95.54 ± 2.52 f	22.43 ± 0.35 e	49.02 ± 2.61 ij	69.22 ± 5.91 j	2.68 ± 0.12 jk	31.71 ± 0.18 h
	N4D3	3201.77 ± 33.25 f	89.22 ± 2.55 g	17.57 ± 0.23 i	26.50 ± 3.18 k	43.56 ± 2.92 k	2.54 ± 0.13 k	27.21 ± 0.15 l
LSD0.05 for
N		**	**	**	**	**	**	**
D		**	**	**	**	**	**	**
N × D		**	ns	**	**	**	**	**
2017	N0D1	2668.04 ± 49.88 i	103.42 ± 2.62 e	13.59 ± 0.12 f	56.89 ± 2.41 c	124.33 ± 4.26 c	3.43 ± 0.08 d	32.44 ± 0.92 def
	N1D1	2959.81 ± 76.43 fg	114.22 ± 1.02 c	17.62 ± 0.26 b	62.89 ± 2.41 b	135.17 ± 5.27 b	4.03 ± 0.10 b	35.13 ± 2.30 b
	N2D1	3057.08 ± 48.14 de	123.31 ± 2.22 b	18.33 ± 0.33 a	73.94 ± 0.82 a	173.19 ± 3.53 a	4.54 ± 0.16 a	36.53 ± 0.18 a
	N3D1	2918.13 ± 57.13 g	137.58 ± 3.74 a	17.53 ± 0.25 b	52.08 ± 0.68 d	124.42 ± 1.23 c	4.16 ± 0.10 b	33.88 ± 0.39 c
	N4D1	2793.06 ± 60.01 h	124.11 ± 2.36 b	16.47 ± 0.18 c	47.48 ± 1.03 e	67.42 ± 2.24 g	3.67 ± 0.09 c	32.13 ± 0.24 ef
	N0D2	2883.35 ± 49.59 g	92.67 ± 1.15 g	11.33 ± 0.54 g	44.42 ± 2.10 ef	82.24 ± 1.76 e	2.93 ± 0.09 ef	30.53 ± 0.35 h
	N1D2	3201.50 ± 16.63 bc	103.94 ± 1.23 e	15.64 ± 0.13 d	46.44 ± 1.07 e	113.56 ± 2.71 d	3.12 ± 0.09 e	34.74 ± 0.51 c
	N2D2	3653.92 ± 17.27 a	105.89 ± 2.34 de	17.54 ± 0.35 b	60.39 ± 2.68 b	122.72 ± 4.75 c	4.07 ± 0.09 b	35.68 ± 0.24 ab
	N3D2	3237.56 ± 41.19 b	122.08 ± 1.54 b	17.32 ± 0.35 b	41.69 ± 1.24 f	84.69 ± 2.89 e	2.99 ± 0.23 ef	32.62 ± 0.04 de
	N4D2	3138.89 ± 48.02 cd	109.67 ± 1.38 d	14.78 ± 0.55 e	34.08 ± 1.63 g	44.50 ± 2.18 i	2.66 ± 0.09 g	30.89 ± 0.68 gh
	N0D3	2763.95 ± 37.47 h	81.86 ± 0.55 h	9.43 ± 0.33 h	30.67 ± 1.76 h	51.19 ± 3.35 h	1.63 ± 0.14 i	29.54 ± 0.11 i
	N1D3	3043.19 ± 36.35 ef	89.78 ± 0.48 g	13.75 ± 0.09 f	32.63 ± 0.86 gh	56.17 ± 2.17 h	1.95 ± 0.06 h	30.14 ± 3.70 hi
	N2D3	3140.38 ± 71.85 cd	92.97 ± 3.89 g	16.42 ± 0.22 c	45.08 ± 3.06 e	73.11 ± 2.99 f	2.82 ± 0.07 fg	33.31 ± 0.08 cd
	N3D3	2973.63 ± 34.65 efg	108.50 ± 4.49 d	15.84 ± 0.46 ed	31.10 ± 1.40 gh	42.82 ± 1.05 i	1.88 ± 0.09 h	31.59 ± 0.14 fg
	N4D3	2885.04 ± 57.65 g	99.08 ± 2.07 f	13.49 ± 0.22 f	21.64 ± 1.54 i	24.81 ± 3.11 j	1.83 ± 0.11 h	29.41 ± 0.31 i
LSD0.05 for
N		**	**	**	**	**	**	**
D		**	**	**	**	**	**	**
N × D		**	ns	**	*	**	**	**

). The grain yield increased from 2726.47 to 3405.84 kg·ha⁻¹ (mean of two growing seasons and three sowing densities) as the nitrogen fertilizer rate increased from 0 to 90 kg·ha⁻¹, and then decreased to 3103.73 kg·ha⁻¹ (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 180 kg·ha⁻¹ resulted in an increase in grain yield from 2828.22 to 3408.63 kg·ha⁻¹ (mean of two growing seasons and five nitrogen fertilizer rates), but it was decreased to 3081.18 kg·ha⁻¹ (mean of two growing seasons and five nitrogen fertilizer rates) when increasing sowing density to 300 kg·ha⁻¹. The highest grain yield occurred in the combination of 90 kg·ha⁻¹ N and 180 kg·ha⁻¹ sowing density (N2D2) in each growing season, and grain yield in 2016 was higher than that in 2017.

3.2.2. Yield Components

Nitrogen fertilizer and planting density had significant effects (p < 0.01) on the plant height during the two growing seasons. The interaction between nitrogen fertilizer and planting density had no significant effect (p > 0.05) in both growing seasons (Table 3). Increasing the nitrogen fertilizer rate from 0 to 135 kg·ha⁻¹ resulted in an increase in plant height from 89.98 to 113.75 cm (mean of two growing seasons and three sowing densities), and then decreased to 104.19 cm (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 300 kg·ha⁻¹ resulted in a decrease in plant height from 113.66 to 90.60 cm (mean of two growing seasons and five nitrogen fertilizer rates). The highest plant height occurred in the combination of 135 kg·ha⁻¹ N and 60 kg·ha⁻¹ sowing density (N3D1) in each growing season.

Nitrogen fertilizer, planting density, and their interaction had significant effects (p < 0.01) on the panicle length during the two growing seasons (Table 3). Increasing the nitrogen fertilizer rate from 0 to 90 kg·ha⁻¹ resulted in an increase in panicle length from 13.46 to 21.77 cm (mean of two growing seasons and three sowing densities), and then a decrease to 16.80 cm (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 300 kg·ha⁻¹ resulted in a decrease in panicle length from 19.57 to 16.72 cm (mean of two growing seasons and five nitrogen fertilizer rates). The highest panicle length occurred in the combination of 90 kg·ha⁻¹ N and 60 kg·ha⁻¹ sowing density (N2D1) in each growing season.

The number of spikelets per panicle (NSP) was significantly (p < 0.05 or p < 0.01) affected by nitrogen fertilizer, planting density, and their interaction during the two growing seasons (Table 3). The NSP increased from 53.57 to 65.36 (mean of two growing seasons and three sowing densities) as the nitrogen fertilizer rate increased from 0 to 90 kg·ha⁻¹, and then decreased to 40.06 (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 300 kg·ha⁻¹ resulted in a decrease in NSP from 67.10 to 39.48 (mean of two growing seasons and five nitrogen fertilizer rates). The highest NSP occurred in the combination of 90 kg·ha⁻¹ N and 60 kg·ha⁻¹ sowing density (N2D1) in each growing season.

The number of grains per panicle (NGP) was significantly (p < 0.01) affected by nitrogen fertilizer, planting density, and their interaction during the two growing seasons (Table 3). Increasing the nitrogen fertilizer rate from 0 to 90 kg·ha⁻¹ resulted in an increase in panicle length from 110.34 to 140.74 (mean of two growing seasons and three sowing densities), and then a decrease to 57.36 (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 300 kg·ha⁻¹ resulted in a decrease in NGP from 143.79 to 67.60 (mean of two growing seasons and five nitrogen fertilizer rates). The highest NGP occurred in the combination of 90 kg·ha⁻¹ N and 60 kg·ha⁻¹ sowing density (N2D1) in each growing season.

Nitrogen fertilizer, planting density, and their interaction had significant effects (p < 0.01) on the weight of grains per panicle (WGP) during the two growing seasons (Table 3). Increasing the nitrogen fertilizer rate from 0 to 90 kg·ha⁻¹ resulted in an increase in WGP from 3.28 to 4.59 g (mean of two growing seasons and three sowing densities), and then a decrease to 3.11 g (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 300 kg·ha⁻¹ resulted in a decrease in WGP from 4.66 to 2.49 g (mean of two growing seasons and five nitrogen fertilizer rates). The highest WGP occurred in the combination of 90 kg·ha⁻¹ N and 60 kg·ha⁻¹ sowing density (N2D1) in each growing season.

Nitrogen fertilizer, planting density, and their interaction had significant effects (p < 0.01) on the 1000-kernel weight during the two growing seasons (Table 3). The 1000-kernel weight increased from 30.75 to 36.08 g (mean of two growing seasons and three sowing densities) as the nitrogen fertilizer rate increased from 0 to 90 kg·ha⁻¹, and then decreased to 29.95 g (mean of two growing seasons and three sowing densities) with further increasing N up to 180 kg·ha⁻¹. Increasing the sowing density from 60 to 300 kg·ha⁻¹ resulted in an increase in 1000-kernel weight from 34.24 to 30.95 g (mean of two growing seasons and five nitrogen fertilizer rates). The highest 1000-kernel weight occurred in the combination of 90 kg·ha⁻¹ N and 60 kg·ha⁻¹ sowing density (N2D1) in each growing season. However, it was not significantly (p > 0.05) higher than the nitrogen fertilization of 90 kg·ha⁻¹ and sowing density of 180 kg·ha⁻¹ treatment (N2D2) in 2017.

3.3. Relationship between Grain Yield and Flag Leaf Photosynthetic Characteristics, as Well as Yield Components

Regression analyses demonstrated that the grain yield was significantly (p < 0.01) correlated with Pn, Gs, Ci, Tr, WUE, Ls, SPAD, LAI, panicle length, NSP, and WGP (Figure 1A–H and Figure 2B,C,E). Remarkable (p < 0.05) correlations were found between the grain yield and NGP, as well as 1000-kernel weight (Figure 2D,F). However, grain yield was not significantly (p > 0.05) correlated with plant height (Figure 2A).

3.4. Relationship between Nitrogen Fertilizer or Planting Density and Flag Leaf Photosynthetic Characteristics, Grain Yield, and Yield Components

Regression analysis showed that the Pn, Gs, Ci, Tr, WUE, Ls, SPAD, LAI, grain yield, plant height, panicle length, number of NSP, NGP, WGP, and 1000-kernel weight were significantly (p < 0.01) correlated with nitrogen fertilizer (

Table 4. Regression of net photosynthetic rate (Pn), stomatal conductance (Gs), intercellular CO₂ concentration (Ci), transpiration rate (Tr), water-use efficiency (WUE), stomatal limitation value (Ls), leaf area index (LAI), grain yield, plant height, panicle length, number of spikelets per panicle (NSP), number of grains per panicle (NGP), weight of grains per panicle (WGP), and 1000-kernel weight of oat on nitrogen fertilizer or planting density; * and ** significant at p < 0.05 and p < 0.01, respectively.

	Nitrogen Fertilizer (x1)	Planting Density (x2)
Pn (y1)	y1 = −0.0003x12 + 0.0625x1 + 10.19 (R2 = 0.263 **)	y1 = −0.0002x22 + 0.0814x2 + 6.74 (R2 = 0.491 **)
Gs (y2)	y2 = −0.00002x12 + 0.0032x1 + 0.23 (R2 = 0.223 **)	y2 = −0.00001x22 + 0.0052x2 − 0.02 (R2 = 0.622 **)
Ci (y3)	y3 = 0.0025x12 − 0.4526x1 + 345.12 (R2 = 0.350 **)	y3 = 0.0011x22 − 0.3825x2 + 359.19 (R2 = 0.251 **)
Tr (y4)	y4 = −0.0002x12 + 0.0297x1 + 2.75 (R2 = 0.240 **)	y4 = −0.0001x22 + 0.0403x2 + 0.90 (R2 = 0.455 **)
WUE (y5)	y5 = 0.000006x12 − 0.0011x1 + 0.38 (R2 = 0.143 **)	y5 = 0.000005x22 − 0.0017x2 + 0.47 (R2 = 0.292 **)
LS (y6)	y6 = −0.00001x12 + 0.0025x1 + 0.32 (R2 = 0.217 **)	y6 = −0.00001x22 + 0.0037x2 + 0.16 (R2 = 0.538 **)
SPAD (y7)	y7 = −0.0013x12 + 0.2388x1 + 41.27 (R2 = 0.237 **)	y7 = −0.0009x22 + 0.3542x2 + 22.95 (R2 = 0.531 **)
LAI (y8)	y8 = 0.000005x12 + 0.0094x1 + 2.95 (R2 = 0.439 **)	y8 = −0.00005x22 + 0.0229x2 + 1.76 (R2 = 0.412 **)
Grain yield (y9)	y9 = −0.051x12 + 11.173x1 + 2719.90 (R2 = 0.333 **)	y9 = −0.0315x22 + 12.402x2 + 2197.60 (R2 = 0.409 **)
Plant height (y10)	y10 = −0.0011x12 + 0.2889x1 + 88.60 (R2 = 0.278 **)	y10 = 0.00007x22 − 0.1218x2 + 120.71 (R2 = 0.471 **)
Panicle length (y11)	y11 = −0.0008x12 + 0.157x1 + 13.35 (R2 = 0.433 **)	y11 = −0.000001x22 − 0.0114x2 + 20.26 (R2 = 0.072 *)
NSP (y12)	y12 = −0.0018x12 + 0.248x1 + 52.78 (R2 = 0.225 **)	y12 = 0.0000005x22 − 0.1153x2 + 74.01 (R2 = 0.478 **)
NGP (y13)	y13 = −0.0058x12 + 0.7468x1 + 109.22 (R2 = 0.351 **)	y13 = 0.00003x22 − 0.3285x2 + 163.39 (R2 = 0.444 **)
WGP (y14)	y14 = −0.0001x12 + 0.0221x1 + 3.20 (R2 = 0.131 **)	y14 = −0.00001x22 − 0.0043x2 + 4.96 (R2 = 0.546 **)
1000-kernel weight (y15)	y15 = −0.0006x12 + 0.1063x1 + 30.73 (R2 = 0.623 **)	y15 = −0.00003x22 − 0.0037x2 + 34.56 (R2 = 0.249 **)

). The Pn, Gs, Ci, Tr, WUE, Ls, SPAD, LAI, grain yield, plant height, number of NSP, NGP, WGP, and 1000-kernel weight were significantly (p < 0.01) correlated with planting density. Moreover, the panicle length was weakly (p < 0.05) correlated with planting density (Table 4).

4. Discussion

In this study, we demonstrated significant correlations of grain yield with Pn, Gs, Ci, Tr, WUE, Ls, SPAD, plant height, panicle length, NSP, NGP, WGP, and 1000-kernel weight. Our data showed that the flag leaf photosynthetic capacity and yield components are key factors to determine oat grain yield at the boot stage. Grain yield varied between years, with higher grain yield in a drier year (2016) than in 2017. This is interesting because earlier studies revealed that the period just before anthesis is a very sensitive period in crop and the photothermal quotient (radiation/temperature), and precipitation has a major influence on grain number and, hence, yield [10,26]. The reason may be that both nitrogen fertilizer and climatic conditions affected vegetative and reproductive development of the crop. The reproductive growth stage of oat in our study was mainly in July, while total rainfall in July 2017 was lower by 57 mm and mean temperature was higher by 2.8 °C than the respective long-term averages (30 year average). In contrast, total rainfall in 2016 was similar to (only 3.5 mm higher than) the long-term average. This pattern of the weather led to a more reproductive development in 2016 and more vegetative development in 2017, as evidenced by the change trend of LAI and plant height in 2 years. Similar conclusions were drawn by Anderson and McLean [27]; in their study, the yield response of oat to applied N in Western Australia depended on soil N status, seasonal rainfall, sowing date, seed rate, and cultivar.

Crop yield is determined by the efficiency of photosynthesis, as well as assimilate transport and distribution. Nitrogen nutrition plays a key role in these processes. Nitrogen fertilization contributes to the LAI, radiation interception, radiation use efficiency, dry matter partitioning to reproductive organs, and protein content of the plant and the seed [28]. Previous studies illustrated that nitrogen deficiency resulted in reductions in leaf photosynthetic capacity and grain yield in Triticum aestivum [29], Zea mays [30], and Sorghum bicolor [31]. In our experiment, the flag leaf photosynthetic capacity, grain yield, and yield components increased with the increase of nitrogen fertilizer from 0 to 90 kg·ha⁻¹, and then decreased with the increase of nitrogen fertilizer from 90 to 180 kg·ha⁻¹. Both the number of spikelets per panicle and number of grains per panicle can affect the critical yield components and number of grains per unit area, thereby influencing grain yield. Increased floret fertility within spikelets in response to nitrogen fertilization, resulting in an increased number of grains, was also reported by Peltonen-Sainio and Järvinen [8]. Reduced floret abortion and free kernel at greater nitrogen availability could increase grain yield [32]. The reason for the decrease in the flag leaf photosynthetic capacity, grain yield, and yield components may be that the weak light irradiance to the leaf as a result of shading due to abundant nitrogen fertilizer promoted leaf development, and oat grew excessively tall when nitrogen fertilizer was abundant and became susceptible to lodging, as proven by the change trend of plant height, especially in 2017, a phenomenon that was also reported by Vos, Van Der Putten, and Birch [14]. Previous studies on cereals reported that yield is primarily determined by grain numbers rather than by grain weight [33,34]. In the present study, regression analysis showed that the NGP and WGP were significantly correlated with nitrogen fertilizer, with a larger coefficient of determination (R²) value of nitrogen fertilizer with NGP than with WGP, in agreement with previous reports on oat [35].

Planting density is one of the most important factors coordinating the contradiction between crop group and individual [7]. In general, increasing plant population produces a greater grain yield for most crops, but high planting density could also cause a yield decline. In this study, the leaf photosynthetic capacity and grain yield increased to their maximums at a sowing density of 180 kg·ha⁻¹, while they subsequently decreased at sowing density of 300 kg·ha⁻¹. The reason for the decrease in leaf photosynthetic capacity and grain yield might be the weak light irradiance to the leaf by shading and strong competition at high planting density, as noted in the previous study [36]. While the yield components per plant decreased with the increase in planting density, maximum yield components appeared at the sowing density of 60 kg·ha⁻¹. Obviously, low competition at low planting density allowed oat a better use of available light for enhancing described yield component parameters. According to Peltonen-Sainio and Järvinen [8], in pure oat stands, an increase in tillering, grains per panicle, and harvest index with decreasing seeding rate was explained by more space per plant. More fertile spikelets per spike and grains per spike of wheat at lower crop density are attributed to less light competition [37]. Increased light availability enhances the differentiation rate of inflorescences during tillering stage, resulting in a higher number of spikelets (and florets) per inflorescence [16]. In the present study, we showed higher coefficient of determination (R²) values for the regression equation of planting density with the flag leaf photosynthetic capacity, grain yield, and yield components (except Ci, LAI, panicle length, and 1000-kernel weight) than of nitrogen fertilizer rates with these parameters. Although the nitrogen fertilizer and planting density had significant effects on the flag leaf photosynthetic capacity, grain yield, and yield components (except LAI and plant height), there was a larger effect of planting density than of nitrogen fertilizer, according to the coefficient of determination (R²) value. New paradigms and further research are needed to enable growers to achieve high grain yield and other characteristics required by millers, while recognizing that nitrogen does not have an overriding effect.

Variety, management, and environment (site and season) are the three groups of factors interacting to determine development and growth in crops and the structure of the crop and the panicles of oat [9]. Nitrogen fertilization and planting density are two of the most important crop management measures. In the present study, it is suggested that appropriate nitrogen fertilizer and planting density could effectively improve the leaf photosynthetic capacity and enhance the yield of oat. The combination which showed the greatest effect was nitrogen fertilization of 90 kg·ha⁻¹ and sowing density of 180 kg·ha⁻¹.

5. Conclusions

Nitrogen fertilizer and planting density had significant effects on flag leaf photosynthetic capacity, grain yield, and yield components of oat. The flag leaf photosynthetic capacity and grain yield increased and then decreased with the increase in nitrogen fertilizer and planting density. The yield components increased and then decreased with the increase in nitrogen fertilizer, while they decreased with the increase in planting density. Specifically, the combination of nitrogen fertilization of 90 kg·ha⁻¹ and sowing density of 180 kg·ha⁻¹ is recommended for oat production on a cool semiarid plateau or other agroecozones with similar environmental conditions.