Grain yield and yield-related traits
Grain yield and yield-related traits varied to some extent under different N application treatments (Table 1). The grain yield of MD treatment was remained the highest in both years. Mean grain yield of both rice cultivars under MD were 8.03–8.35 t ha− 1, which was 33.48–36.35% higher than the MB treatment. The MD treatment resulted the highest number of productive panicles ha− 1 and spikelet per panicle with 257.16 × 104 − 274.21 × 104 ha− 1 and 178.25–185.96, respectively, which was significantly higher than MB treatment. The MD and MDB treatment were not differed significantly (P > 0.05) in number of productive panicles ha− 1 and spikelet per panicle. Substantial difference was found in the grain yield and yield-related traits between nitrogen treatments (p < 0.05). Furthermore, the N × Y factor also had prominent effect on spikelet per panicle and harvest index (HI).
Table 1
Effects of mechanized deep placement of nitrogen fertilization on yield and its components of rice during 2019 and 2020
Treatment
|
Productive
|
Spikelet
|
Grain
|
1000-grain-
|
Harvest
|
Harvested
|
panicle
|
per
|
filling
|
weight
|
index
|
Yield
|
(104 ha− 1)
|
panicle
|
(%)
|
(g)
|
|
(t ha− 1)
|
2019
|
WFY615
|
|
|
|
|
|
|
|
CK
|
130.95 ± 5.92b
|
139.57 ± 6.65b
|
76.13 ± 0.08a
|
21.62 ± 0.31b
|
0.66 ± 0.03a
|
4.23 ± 0.33d
|
|
MDS
|
253.97 ± 3.53a
|
158.62 ± 3.07ab
|
80.22 ± 0.04a
|
23.34 ± 0.11a
|
0.66 ± 0.02a
|
7.17 ± 0.45b
|
|
MBS
|
277.78 ± 4.33a
|
180.19 ± 4.93a
|
83.15 ± 0.03a
|
23.74 ± 0.28a
|
0.60 ± 0.05a
|
8.09 ± 0.33a
|
|
MD
|
170.63 ± 7.92b
|
147.30 ± 6.92b
|
81.30 ± 0.04a
|
22.14 ± 0.61b
|
0.62 ± 0.06a
|
5.45 ± 0.49c
|
|
mean
|
208.33
|
156.42
|
80.20
|
22.71
|
0.64
|
6.24
|
|
YXYZ
|
|
|
|
|
|
|
|
CK
|
130.95 ± 4.54c
|
149.14 ± 5.77b
|
83.28 ± 0.02a
|
19.21 ± 0.12c
|
0.64 ± 0.02a
|
3.37 ± 0.33c
|
|
MDS
|
234.12 ± 7.92ab
|
160.17 ± 4.24ab
|
86.55 ± 0.02a
|
19.68 ± 0.43a
|
0.6 ± 0.05ab
|
7.09 ± 0.57a
|
|
MBS
|
250.24 ± 3.81a
|
176.41 ± 1.78a
|
86.42 ± 0.05a
|
21.71 ± 0.47a
|
0.50 ± 0.06c
|
7.67 ± 0.33a
|
|
MD
|
186.51 ± 6.01b
|
157.90 ± 2.99ab
|
82.11 ± 0.05a
|
20.75 ± 0.25b
|
0.55 ± 0.01b
|
5.80 ± 0.37b
|
|
mean
|
200.40
|
160.91
|
84.59
|
20.34
|
0.57
|
5.98
|
2020
|
WFY615
|
|
|
|
|
|
|
|
CK
|
184.11 ± 4.57c
|
140.99 ± 6.58b
|
78.81 ± 0.48b
|
22.84 ± 0.25ab
|
0.58 ± 0.03a
|
5.62 ± 0.15c
|
|
MDS
|
264.43 ± 2.98ab
|
173.32 ± 2.97ab
|
81.51 ± 0.41a
|
23.16 ± 0.23a
|
0.57 ± 0.01a
|
7.57 ± 0.74b
|
|
MBS
|
270.67 ± 3.86a
|
191.79 ± 3.77a
|
81.23 ± 0.17a
|
22.24 ± 0.11b
|
0.56 ± 0.01a
|
8.60 ± 0.11a
|
|
MD
|
220.37 ± 2.98bc
|
160.84 ± 5.81ab
|
80.12 ± 0.86ab
|
22.37 ± 0.64ab
|
0.57 ± 0.01a
|
6.92 ± 0.26b
|
|
mean
|
234.76
|
166.74
|
80.42
|
22.65
|
0.57
|
7.18
|
|
YXYZ
|
|
|
|
|
|
|
|
CK
|
150.09 ± 4.21c
|
131.05 ± 1.31c
|
77.83 ± 1.06c
|
19.75 ± 0.56a
|
0.48 ± 0.02c
|
4.76 ± 0.68c
|
|
MDS
|
244.25 ± 3.86a
|
176.59 ± 6.31ab
|
79.67 ± 0.13b
|
19.87 ± 0.19a
|
0.52 ± 0.03a
|
7.67 ± 0.07a
|
|
MBS
|
264.07 ± 5.60a
|
180.18 ± 2.55a
|
80.77 ± 0.28a
|
21.61 ± 0.92a
|
0.49 ± 0.01bc
|
8.38 ± 0.94a
|
|
MD
|
208.39 ± 6.97b
|
149.74 ± 3.34bc
|
80.15 ± 0.76ab
|
20.01 ± 1.65a
|
0.50 ± 0.01b
|
6.22 ± 0.95b
|
|
mean
|
216.50
|
159.39
|
79.61
|
20.31
|
0.50
|
6.76
|
Anova
|
Year(Y)
|
ns
|
ns
|
ns
|
*
|
ns
|
ns
|
|
Cultivar(C)
|
ns
|
ns
|
ns
|
**
|
**
|
ns
|
|
Nitrogen(N)
|
**
|
**
|
**
|
**
|
**
|
**
|
|
Y × C
|
ns
|
ns
|
ns
|
**
|
ns
|
ns
|
|
Y × N
|
ns
|
*
|
ns
|
ns
|
**
|
ns
|
|
N × C
|
ns
|
ns
|
ns
|
ns
|
ns
|
ns
|
|
Y × N × C
|
ns
|
ns
|
*
|
ns
|
ns
|
ns
|
Note: WFY615, Wufemgyou615; YXYZ, Yuxiangyouzhan. Within a column, means followed by the same letter are not significantly different at the 0.05 probability level according to least significant different test (LSD0.05). The same as below. |
Root morphological indexes
The root morphological indexes were different to some extent under different N application treatments at the HS stage (Table 2). The root morphological indexes for MD treatment were substantially improved for both rice cultivars. Mean total root volume (TRV) of both rice cultivars under MD were 17.81–21.59 cm3 per hill, and the total root length (TRL) was 69.79–77.46 m per hill, root superficial area (RSA) was 12.94–13.29 cm2 per hill, respectively. Moreover, during the two early seasons, difference in TRV, TRL, and RSA between the MD and MB treatments was significant, however, no significant difference was found between deep fertilizer treatment. The highest root dry weight (RDW) was recorded in MD treatment, which was significantly higher 17.21–24.11% than MB treatment. Significant differences were found in the root morphological traits and RDW between N treatments (P < 0.05), moreover, the Y × C factor also had obvious impact on TRV and RSA. The Y × N × C factor had also prominent impact on RDW and ARD.
Table 2
Effects of mechanized deep placement of nitrogen fertilization on root morphological indexes at heading stage during 2019 and 2020
Treatment
|
TRV
|
TRL
|
RDW
|
RSA
|
ARD
|
(cm3/hill)
|
(m/hill)
|
(g/hill)
|
(cm2/hill)
|
(mm/hill)
|
2019
|
WFY615
|
|
|
|
|
|
|
CK
|
11.33 ± 0.87c
|
56.63 ± 1.77c
|
6.22 ± 0.88b
|
9.38 ± 0.44c
|
9.55 ± 0.45b
|
|
MDB
|
17.04 ± 4.11b
|
68.22 ± 8.34b
|
10.11 ± 0.80a
|
12.27 ± 2.08ab
|
10.40 ± 0.03a
|
|
MD
|
22.47 ± 1.83a
|
79.25 ± 2.83a
|
10.35 ± 0.71a
|
13.02 ± 0.71a
|
10.61 ± 0.31a
|
|
MB
|
13.37 ± 1.39bc
|
60.79 ± 3.71bc
|
6.88 ± 1.02b
|
10.42 ± 0.92bc
|
8.95 ± 0.23b
|
|
mean
|
16.05
|
66.22
|
8.39
|
11.27
|
9.88
|
|
YXYZ
|
|
|
|
|
|
|
CK
|
11.69 ± 1.24c
|
56.36 ± 4.01c
|
5.97 ± 0.51c
|
9.57 ± 1.12b
|
9.62 ± 0.23ab
|
|
MDB
|
15.29 ± 1.36b
|
63.67 ± 2.51b
|
9.89 ± 0.95b
|
11.39 ± 0.62a
|
9.39 ± 0.47b
|
|
MD
|
19.39 ± 1.98a
|
72.99 ± 1.36a
|
11.36 ± 0.54a
|
11.46 ± 0.34a
|
10.36 ± 0.52a
|
|
MB
|
12.24 ± 0.67bc
|
56.48 ± 2.76c
|
9.38 ± 1.22b
|
10.15 ± 0.69ab
|
9.43 ± 0.47b
|
|
mean
|
14.65
|
62.38
|
9.15
|
10.64
|
9.70
|
2020
|
WFY615
|
|
|
|
|
|
|
CK
|
14.68 ± 1.09c
|
63.43 ± 2.21c
|
5.45 ± 1.04b
|
11.08 ± 0.55c
|
9.71 ± 0.15c
|
|
MDB
|
21.15 ± 1.96a
|
73.57 ± 1.68ab
|
9.33 ± 0.99a
|
14.35 ± 0.42a
|
12.57 ± 0.36a
|
|
MD
|
20.70 ± 1.70a
|
75.66 ± 3.98a
|
9.91 ± 1.70a
|
15.12 ± 0.99a
|
9.45 ± 0.26c
|
|
MB
|
17.81 ± 0.83b
|
69.77 ± 3.46b
|
8.17 ± 0.93ab
|
12.65 ± 0.86b
|
10.77 ± 0.03b
|
|
mean
|
18.58
|
70.61
|
8.22
|
13.3
|
10.63
|
|
YXYZ
|
|
|
|
|
|
|
CK
|
10.48 ± 0.55c
|
54.91 ± 1.12b
|
6.12 ± 1.48b
|
8.96 ± 0.28b
|
9.57 ± 0.21b
|
|
MDB
|
14.53 ± 2.96ab
|
66.13 ± 2.52a
|
9.13 ± 1.16a
|
12.00 ± 0.63a
|
10.50 ± 0.24a
|
|
MD
|
16.23 ± 1.38a
|
66.58 ± 6.02a
|
9.07 ± 0.27a
|
12.86 ± 1.5a
|
9.56 ± 0.33b
|
|
MB
|
11.72 ± 1.24bc
|
57.42 ± 2.81b
|
7.08 ± 1.58ab
|
9.58 ± 0.70b
|
9.39 ± 0.31b
|
|
mean
|
13.24
|
61.26
|
7.85
|
10.85
|
9.76
|
Anova
|
Year(Y)
|
ns
|
ns
|
ns
|
**
|
ns
|
|
Cultivar(C)
|
ns
|
ns
|
ns
|
ns
|
ns
|
|
Nitrogen(N)
|
**
|
**
|
**
|
**
|
**
|
|
Y × C
|
**
|
ns
|
ns
|
**
|
ns
|
|
N × C
|
ns
|
ns
|
ns
|
ns
|
**
|
|
Y × N
|
ns
|
**
|
ns
|
ns
|
ns
|
|
Y × N × C
|
ns
|
ns
|
*
|
ns
|
**
|
Note: TRV: total root volume, TRL: total root length, RDW: root dry weight, RSA: root superfcial area, ARD: average root diameter. |
Detection of chlorophyll contents and net photosynthetic rate (Pn)
The total Chl at whole growth stage and Pn at HS stage were differed to some extent under different N application treatments (Table 3). For example, the MD treatment resulted in the highest total Chl for both rice cultivars, and remained significantly higher than other treatments at MT stage. Total Chl content was remained lower at PI stage than MT stage, with marginal difference among treatments whereas total Chl was noted the lowest at the HS stage for each treatment. Mean total Chl at HS stage of both rice cultivars under MD were 2.47–3.06 mg g− 1, which was 25.15–55.12% higher than MB treatment, respectively. The highest Pn was obtained for MD treatment for both rice cultivars, and mean Pn under MD were 20.27–22.17 µmol m− 2 s− 1, which was 24.63–25.57% higher than MB treatment. Furthermore, the similar trend in Pn was observed i.e., MD > MDB > MB > CK.
Table 3
Effects of mechanized deep placement of nitrogen fertilization on chlorophyll contents and net photosynthetic rate during 2019 and 2020
Treatment
|
Total Chl (mg·g− 1)
|
Net photosynthetic
|
MT
|
PI
|
HS
|
rate (µmol m− 2s− 1)
|
2019
|
WFY615
|
|
|
|
|
|
CK
|
1.38 ± 0.07c
|
1.90 ± 0.05b
|
0.91 ± 0.07c
|
13.60 ± 0.76c
|
|
MDB
|
3.09 ± 0.16b
|
3.28 ± 0.06a
|
3.06 ± 0.29a
|
18.89 ± 0.19b
|
|
MD
|
3.94 ± 0.50a
|
3.36 ± 0.19a
|
3.00 ± 0.03a
|
22.45 ± 0.32a
|
|
MB
|
3.11 ± 0.18b
|
2.03 ± 0.07b
|
2.49 ± 0.27b
|
17.81 ± 1.65b
|
|
mean
|
2.88
|
2.64
|
2.37
|
18.19
|
|
YXYZ
|
|
|
|
|
|
CK
|
2.27 ± 0.12c
|
1.81 ± 0.09c
|
1.27 ± 0.05b
|
12.18 ± 0.96c
|
|
MDB
|
2.77 ± 0.07b
|
3.57 ± 0.14a
|
1.04 ± 0.05b
|
19.14 ± 0.18a
|
|
MD
|
3.36 ± 0.16a
|
3.84 ± 0.18a
|
2.27 ± 0.54a
|
20.29 ± 0.99a
|
|
MB
|
2.27 ± 0.02c
|
3.12 ± 0.08b
|
1.50 ± 0.14b
|
16.02 ± 0.41b
|
|
mean
|
2.67
|
3.09
|
1.52
|
16.91
|
2020
|
WFY615
|
|
|
|
|
|
CK
|
1.45 ± 0.25c
|
2.05 ± 0.09c
|
1.01 ± 0.06c
|
15.58 ± 1.36c
|
|
MDB
|
3.18 ± 0.11b
|
3.11 ± 0.03b
|
2.98 ± 0.11a
|
21.48 ± 1.72ab
|
|
MD
|
3.7 ± 0.07a
|
3.88 ± 0.06a
|
3.12 ± 0.15a
|
23.89 ± 1.92ab
|
|
MB
|
2.9 ± 0.15b
|
2.97 ± 0.1b
|
2.4 ± 0.07b
|
19.50 ± 0.25b
|
|
mean
|
2.81
|
3.01
|
2.38
|
18.11
|
|
YXYZ
|
|
|
|
|
|
CK
|
2.20 ± 0.18d
|
1.98 ± 0.31b
|
0.96 ± 0.11c
|
11.18 ± 0.8d
|
|
MDB
|
3.13 ± 0.06b
|
3.27 ± 0.18a
|
1.81 ± 0.2b
|
19.11 ± 0.8b
|
|
MD
|
3.35 ± 0.1a
|
3.42 ± 0.48a
|
2.66 ± 0.04a
|
20.24 ± 0.4a
|
|
MB
|
2.4 ± 0.08c
|
3.23 ± 0.06a
|
1.67 ± 0.14b
|
16.51 ± 0.34c
|
|
mean
|
2.77
|
2.98
|
1.78
|
16.76
|
Anova
|
Year(Y)
|
ns
|
ns
|
ns
|
ns
|
|
Cultivar(C)
|
ns
|
ns
|
**
|
**
|
|
Nitrogen(N)
|
**
|
**
|
**
|
**
|
|
Y × C
|
ns
|
**
|
*
|
**
|
|
Y × N
|
*
|
**
|
ns
|
ns
|
|
N × C
|
**
|
**
|
**
|
ns
|
|
Y × N × C
|
ns
|
ns
|
ns
|
ns
|
Note: MT: mid-tillering stage; PI: panicle initiation stage; HS: heading stage; MS: maturity stage; The same as below. |
Total above-ground biomass
Significant differences were noted in the total above-ground biomass at whole growth stage under different N application treatments (Fig. 2). For example, the TAB for deep fertilizer treatments (MD and MDB) were significantly higher than MB and CK treatment at PI and HS stage. However, no significant difference was found between MD and MDB. Mean TAB at MS stage of both rice cultivars under MD were 14.11–14.66 t ha− 1, which was 21.17–31.99% higher than MB treatment, respectively. Moreover, the MD and MDB treatment was remained statistically similar (P > 0.05) at MS stage. In the whole growth period, the similar trend for TAB of both rice cultivars was observed. It was also found that TAB was remained the highest at MS stage for each treatment than at other stages.
Nutrient accumulation
Significant differences were also observed in the total N accumulation (TNA) at whole growth stage under different N application treatments (Fig. 3). For instance, there were no significant differences in MD, MDB and MB treatments at MT stage, however, all treatments were substantially higher than CK treatment. The MD treatment had the highest TNA at PI stage, as compared to other treatments. Moreover, the MD treatment significantly increased the TNA in 2019, however, no significant difference was found between MD and MDB treatment in 2020. Mean TNA at MS stage of both rice cultivars under MD were 186.01-190.04 kg ha− 1, which was 26.38–44.15% higher than MB treatment, respectively. On the other hand, the MD and MDB treatments statistically similar (P > 0.05) at MS stage. In the whole growth period, similar trends were observed for TNA of both rice cultivars.
The total P accumulation (TPA) at whole growth stage were varied to some extent under different N applications treatments (Fig. 4). For example, no significant differences were noted in MD, MDB and MB treatments at MT stage for both rice cultivars, however, all treatments were remained significantly higher than CK. The MD treatment had the highest TNA at PI stage, as compared with other treatments. Moreover, the TPA as increased in MD treatment for WFY615, however, no significant difference was found between MD and MDB treatment in YXYZ regarding TPA. Furthermore, deep fertilizer treatments i.e., MDB and MD were remained statistically similar (P > 0.05) at HS and MS stage. Mean TPA at MS stage of both rice cultivars under MD were 32.82–35.5 kg ha− 1, which was remarkable higher 28.72–30.23% than MB treatment, respectively. In the whole growth period, the similar trend for TPA of both rice cultivars was observed.
The total K accumulation (TKA) at whole growth stage were differed to some extent under different N applications treatments (Fig. 5). For example, there were no significant differences in MD, MDB and MB treatments at MT stage for both rice cultivars, however, all treatments were remained higher than CK treatment. The MD treatment had the highest TKA at PI stage, as compared to other treatments, whereas the TKA in WFY615 were increased under MD treatment in 2019, however, no remarkable difference was found between MD and MDB treatment in YXYZ. The MD treatment led to the highest TKA at MS stage, Furthermore, mean TKA at MS stage of both rice cultivars under MD were 173.76-182.63 kg ha− 1, which was 25.61–37.33% higher than MB treatment, respectively. Similar trend for TKA of both rice cultivars was observed during whole growth period.
Glutamine synthetase (GS) and nitrate reductase (NR) activities
The nitrate reductase (NR) activity at all growth stages were remained variable to some extent under different N applications treatments (Fig. 6). For WFY615, the NR activity of MD treatment was remained the highest among all treatments at MT stage during both years, however, MB and MDB treatment was not differed significantly. For both rice cultivars, as compare to other treatments, the MD treatment led to the highest NR activity at HS stage, furthermore, mean NR activity at HS stage of both rice cultivars under MD were remained 111.59-118.24% higher than MB treatment.
The glutamine synthetase (GS) activity at all growth stages differed under different N applications treatments (Fig. 7). For both rice cultivars, the GS activity of MD treatment was the highest among all treatment at MT stage, while significant difference was found among MD and MDB treatment. However, MD and MDB treatment was not differed significantly (P > 0.05). Compared with other treatments, the MD treatment led to the highest GS activity at HS stage. Furthermore, mean GS activity at HS stage of both rice cultivars under MD were remained 16.97–23.36% higher than MB treatment. Moreover, significant differences in GS activity were found in all treatments of YXYZ in 2020.
Economic profitability
The total input, total return, net income, and BCR were differed under different N applications treatments (Table 5). Significant differences in total input were found among all treatments, while the highest total inputs in both seasons were observed in MB. However, MD generated maximum net income and total return, next to MDB and MB, while the lowest net income was observed in CK. The BCR was the highest in MD than all other treatments because of the highest total return and lower total input. Thus, present results demonstrated that mechanized deep N fertilization was a suitable and an economical fertilization technology for mechanically transplanted rice.
Correlation among root morphological indexes, nutrients accumulation, total above-ground biomass, and yield
The correlation coefficients amongst the root morphological indexes and nutrient accumulation i.e., TNA, TPA, and TKA, as well as TAB and yield (Table 4). The TNA had significant and positive correlations with root morphological indexes such as TRL, RSA, TRV, and RDW for both rice cultivars, correlations between TNA and ARD were remained non-significant. Moreover, the TPA had significant and positive correlations with TRV and RDW whereas the TKA has also significant associations with root morphological indexes except ARD. Moreover, significant and positive correlation among root morphological indexes, TAB and yield at was also noted at P<0.01.
Table 4
Correlation analysis between root morphological indexes and nutrients accumulation, total above-ground biomass, and yield
|
TRV
|
RSA
|
RDW
|
TRL
|
ARD
|
|
(cm3/hill)
|
(cm2/hill)
|
(g/hill)
|
(m/hill)
|
(mm/hill)
|
TNA (kg ha-1)
|
0.4877*
|
0.5071*
|
0.4177*
|
0.5136**
|
0.2812
|
TPA (kg ha-1)
|
0.4497*
|
0.2113
|
0.6232**
|
0.3496
|
0.1532
|
TKA (kg ha-1)
|
0.6157**
|
0.5874**
|
0.6473**
|
0.6299**
|
0.1783
|
TAB (t ha-1)
|
0.7476**
|
0.7375**
|
0.6359**
|
0.7426**
|
0.4305*
|
Yield (t ha-1)
|
0.7343**
|
0.7639**
|
0.6497**
|
0.7607**
|
0.5571**
|
Note: TRV: total root volume, TRL: total root length, RDW: root dry weight, RSA: root superficial area, ARD: average root diameter, TNA: total nitrogen accumulation, TPA: total phosphorus accumulation, TKA: total potassium accumulation, TAB: total above-ground biomass.
|
Table 5
Production cost, gross return, and beneficial cost ratio (BCR) in 2019 and 2020
Treatments
|
Production cost (US$ha−1)
|
Gross return (US$ha−1)
|
Net income
|
Fertilizer
|
Fertilizer
|
Machinery
|
Labor
|
Input
|
Total
|
Price
|
GY
|
Total
|
(US$ha−1)
|
BCR
|
application
|
costs
|
costs
|
costs
|
input
|
(US$ t-1)
|
(t ha−1)
|
return
|
2019
|
|
|
|
|
|
|
|
|
|
|
|
CK
|
0
|
0
|
314.93
|
354.16
|
177.08
|
846.17d
|
390.12
|
3.82d
|
1489.80d
|
643.63d
|
0.76b
|
MDB
|
165.27
|
47.03
|
314.93
|
354.16
|
177.08
|
893.20b
|
390.12
|
7.40b
|
2886.00b
|
1992.80b
|
2.23a
|
MD
|
165.27
|
34.58
|
314.93
|
354.16
|
177.08
|
880.75c
|
390.12
|
7.88a
|
3073.20a
|
2192.45a
|
2.49a
|
MB
|
165.27
|
76.09
|
314.93
|
354.16
|
177.08
|
922.26a
|
39012
|
5.625c
|
2193.75c
|
1271.49c
|
1.38b
|
2020
|
|
|
|
|
|
|
|
|
|
|
|
CK
|
0
|
0
|
314.93
|
354.16
|
177.08
|
846.17d
|
390.12
|
4.36d
|
1700.92d
|
854.75d
|
1.01d
|
MDB
|
165.27
|
47.03
|
314.93
|
354.16
|
177.08
|
893.20b
|
390.12
|
7.62b
|
2972.71b
|
2079.51b
|
2.33b
|
MD
|
165.27
|
34.58
|
314.93
|
354.16
|
177.08
|
880.75c
|
390.12
|
8.49a
|
3312.12a
|
2431.37a
|
2.76a
|
MB
|
165.27
|
76.09
|
314.93
|
354.16
|
177.08
|
922.26a
|
390.12
|
6.57c
|
2563.08c
|
1640.82c
|
1.78d
|
Note: The grain yields were represented the average of two rice cultivars
|