Comparative Yield, Fiber Quality and Dry Matter Production of Cotton Planted at Various Densities under Equidistant Row Arrangement

Abstract

The number of cotton plants grown per unit area has recently gained attention due to technology expense, high input, and seed cost. Yield consistency across a series of plant populations is an attractive cost-saving option. Field experiments were conducted to compare biomass accumulation, fiber quality, leaf area index, yield and yield components of cotton planted at various densities (D1, 1.5; D2, 3.3; D3, 5.1; D4, 6.9; D5, 8.7; and D6, 10.5 plants m⁻²). High planting density (D5) produced 21% and 28% more lint yield as compared to low planting density (D1) during both years, respectively. The highest seed cotton yield (4662 kg/ha) and lint yield (1763 kg/ha) were produced by high plant density (D5) while the further increase in the plant population (D6) decreased the yield. The increase in yield of D5 was due to more biomass accumulation in reproductive organs as compared to other treatments. The highest average (19.2 V_A gm m⁻² d⁻¹) and maximum (21.8 V_M gm m⁻² d⁻¹) rates of biomass were accumulated in reproductive structures. High boll load per leaf area and leaf area index were observed in high planting density as compared to low, while high dry matter partitioning was recorded in the lowest planting density as compared to other treatments. Plants with low density had 5% greater fiber length as compared to the highest plant density, while the fiber strength and micronaire value were 10% and 15% greater than the lowest plant density. Conclusively, plant density of 8.7 plants m⁻² is a promising option for enhanced yield, biomass, and uniform fiber quality of cotton.

1. Introduction

Cotton is an important cash crop grown worldwide as a major source of fiber [1]. Cotton is perennial but commercially grown as an annual crop and has indeterminate growth. China is the largest cotton-producing country in the world by contributing about 30% of the world’s cotton production [2]. Henan Province is one of the major cotton growing provinces of China, with more than 400 thousand hectares of land [3]. Plant density determination is one of the most important practices for increasing yield of cotton [4]. Plant density is the key factor for optimizing structures and increasing the photosynthetic capacity of the cotton canopy. High planting density has become common in cotton production systems. It has been reported that both too high and too low plant density reduces cotton yield by affecting light penetration and moisture availability, further influencing plant height, architecture, boll behavior, and crop maturity. An optimum plant density not only improves the yield and fiber quality of cotton but also reduces input costs by minimizing seed rate and fertilizer application without decreasing yield [5]. Low plant density produced a higher number of heavy bolls per plant, while both the number and weight of bolls reduced with increasing plant density [6,7]. Currently, suggested and practiced plant densities in China are 5.3 × 10⁴ to 7.5 × 10⁴ plants ha⁻¹ in the Yellow River Valley [8], 3.0 × 10⁴ plants ha⁻¹ in the Yangtze River Valley [9], and 22.7 × 10⁴ plants ha⁻¹ in the Northwest region. The difference between the plant densities among various locations is due to difference in climatic conditions which affect the yield and fiber quality of cotton.

Biomass accumulation in the cotton plant during the early growth period is an important factor for final yield determination. More biomass accumulation in early stages helps in better establishment of a crop while accumulation at late growth stages increases assimilation to the reproductive organs, resulting in a higher yield and quality of cotton [10]. Cotton plants accumulate more biomass in vegetative organs due to its indeterminate nature. More assimilate accumulation to vegetative and reproductive organs increases the shedding of fruit and leaves [11,12]. At maturity, the aboveground biomass becomes lower than the total due to the shedding of leaves and fruits [13]. Previous studies have confirmed that optimum plant density is the critical factor for establishing optimal canopy structure and leaf area index (LAI). Optimal LAI determines light penetration in the canopy [14,15,16]. Several researchers have examined the relationship between the plant density, LAI, and cotton production [17,18,19] and found that an increase in plant density results in higher LAI, while too-high LAI caused shading and reduced the yield [20,21]. Both LAI and yield increases slowly with an increase in plant density [22]. Fiber quality indicators including fiber strength, fitness, length, uniformity index, and fineness are negatively affected by environmental and genetic factors as well as poor management practices at flowering and boll formation stages [23,24]. Similarly, fiber quality is affected by plant density, irrigation, fertilization, and weather changes [23,25]. This study is conducted with the aim to assess the response of cotton yield and fiber quality, biomass accumulation, and partitioning of various plant densities to identify technological alternatives to make efficient use of land and increase yield and profitability of cotton.

2. Materials and Methods

2.1. Experimental Site

The study was conducted in 2016 and 2017 at the experimental station of the Institute of Cotton Research of Chinese Academy of Agricultural Sciences in Anyang, Henan, China (36°06′ N, 114°21′ E). The soil was medium loam in texture with a total N of 0.65 g kg⁻¹, P of 0.01 g kg⁻¹ and K of 0.11 g kg⁻¹. The monthly average temperature and relative humidity data of both years of cotton growing seasons are presented in Figure 1. The average temperature during the cotton growing season was 22 °C and 23 °C in 2016 and 2017, respectively. Annual rainfall was 713 mm in 2016 and 585 mm in 2017. Annual sunshine hours were 1737 h in 2016 and 1838 h in 2017. The average air temperature at the seedling and reproductive stages was cooler as compared to other growth stages. The overall cotton growing season in 2016 was cooler with more rainfall as compared to 2017.

2.2. Experimental Design

The experiment was conducted in a randomized complete block design (RCBD). Six plant densities (D1, 1.5; D2, 3.3; D3, 5.1; D4, 6.9; D5, 8.7; and D6, 10.5 plants m⁻²) were plotted randomly in three replications on clay loam soil. Each experimental plot area was 64 m² with 8 m in length and width. Each plot consisted of 10 rows, with a row spacing of 0.8 m, which was constant for all plant densities. Seeds of cotton mid maturity cultivar SCRC28 were sown by hand on flat beds with plastic mulching to conserve soil moisture from evaporation. Plastic mulch was removed after one month of full emergence. Seedlings were thinned to the required plant densities after three weeks of emergence. During both years, the land was prepared by ploughing, and irrigated in early spring before sowing. Sowing was done during the growing season on 22 April in 2016 and 2017.

A basal dose of 225 kg N ha⁻¹, 150 kg P₂O₅ ha⁻¹, and 225 kg K₂O ha⁻¹ was applied to the field before sowing. Irrigation was applied by flooding during the flowering stage at a total volume of approximately 45 m³. Crop management practices such as weeding, hoeing, pesticides, and irrigation were performed in a timely manner to enhance crop growth.

2.3. Data Collection

Data were recorded on cotton leaf area index, biomass accumulation at critical stages of crop growth, fiber quality, yield and yield components (boll m⁻² and boll weight) during 2016 and 2017 at different days after emergence.

2.3.1. Yield and Yield Components

Seed cotton yield (kg/ha) and lint yield (kg/ha) were recorded by hand-harvesting three times from each treatment. The boll moisture was reduced to less than 11% by air-drying and seed cotton of 100 bolls at first harvest were sampled for boll weight. Weight of single boll was calculated by dividing total seed cotton yield of 100 bolls by the total number of bolls. Lint percentage was calculated from lint yield of 100 bolls divided by seed cotton weight of 100 bolls.

2.3.2. Biomass Accumulation and Partitioning

The dry weight of cotton plants was recorded seven times during the growing season with an interval of 15 days at 42 days after emergence (DAE), 57 DAE, 72 DAE, 87 DAE, 102 DAE, 117 DAE, and 132 DAE. Three plants from each plot of three replications were uprooted randomly and dissected into the underground part (roots), leaves, stem, and reproductive structures. Samples were quickly placed for 30 min in an electric fan-assisted oven at 105 °C in order to stop metabolism. Samples were dried at 80 °C for 48 h to attain a constant weight. Dry matter partitioning was calculated by the ratio of the dry weight of reproductive organs (DWRO) (squares, flowers, green, and open bolls) to plant total biomass while boll load was calculated by dividing DWRO by leaf area. A logistic regression equation was used to describe biomass accumulation [26].

$$\mathrm{Y}=\frac{\mathrm{A}}{1+{\mathrm{be}}^{-\mathrm{kt}}}$$

(1)

In Equation (1) Y (kg) is the biomass, A (kg) the maximum biomass, t (d) is the number of days after emergence (DAE) while a and b are constants.

From Formula (1), the following equations were calculated:

$${\mathrm{t}}_{\mathrm{o}}=\frac{\mathrm{lnb}}{\mathrm{k}}({\mathrm{t}}_{\mathrm{o}}=\mathrm{t})$$

(2)

$${\mathrm{t}}_1=\frac{\mathrm{lnb}-\ln \left(2+\sqrt{3}\right)}{\mathrm{k}}$$

(3)

$${\mathrm{t}}_2=\frac{\mathrm{lnb}+\ln \left(2+\sqrt{3}\right)}{\mathrm{k}}$$

(4)

$${\mathrm{V}}_{\mathrm{M}}=\frac{\mathrm{Ak}}{4}$$

(5)

$$∆\mathrm{t}={\mathrm{t}}_2-{\mathrm{t}}_1$$

(6)

$${\mathrm{V}}_{\mathrm{A}}=\frac{{\mathrm{Y}}_2-{\mathrm{Y}}_1}{∆\mathrm{t}}$$

(7)

In the above equations, V_M (kg ha⁻¹ d⁻¹) is the highest rate of biomass accumulation, and t (d) is the maximum biomass fast accumulation period. Y₁ and Y₂ are the biomass at t₁ and t₂. V_A indicates the average biomass accumulation from t₁ to t₂ and ∆t (d) is the total period of average biomass accumulation.

2.3.3. Leaf Area Index

LAI of cotton plants were calculated by taking photos of leaves through a scanning machine (Phantom p800xl, MiCROTEK, Shanghai, China) and leaf area was calculated by using Image-Pro Plus 7.0 (Media Cybernetics, Rockville, MD, USA). The LAI was determined by dividing the total plant leaf area per unit ground area.

2.3.4. Fiber Quality

Fiber quality, including fiber length (mm), fiber uniformity, fiber strength (cN tex⁻¹), and fiber micronaire, were assessed by the Supervision, Inspection and Test Center of Cotton Quality, Ministry of Agriculture, in Anyang, Henan province of China using a high volume instrument (HVI-900) (Changing Technologies, Mainland, China) according to the internationally accepted ICC standard.

2.3.5. Statistical Analysis

Microsoft Excel 365 (Microsoft, Bothell, WA, USA) as used for the processing of data. SPSS 19.0 (SPSS Inc. Chicago, IL, USA) and Origin 2016 (OriginLab Corporation, Northampton, MA, USA) were used for the analysis of data. Figures were plotted by using Origin 2016. Duncan’s multiple range test at 5% probability level was used to test differences among mean values.

3. Results

3.1. Yield and Yield Components

Yield and yield components of cotton varied with plant density. Seed cotton yield and lint yield along with yield components were significantly affected by plant density except boll weight and lint percentage in both years (

Table 1. Comparison of boll m⁻² and boll weight (g) at various plant densities in 2016 and 2017 in the cotton growing season.

Treatment	Boll (m2)	Boll Weight (g)
Year 2016
Plant Density (PD)
D1	64.3f	6.2a
D2	72.7e	5.8a
D3	82.4d	5.7a
D4	90.1c	5.7a
D5	99b	5.7a
D6	104.4a	5.6a
Year 2017
Plant Density (PD)
D1	46.5e	6a
D2	51.3d	6a
D3	59.5c	5.9a
D4	64.4b	5.9a
D5	75.7a	5.7ab
D6	66.4b	5.6b
ANOVA
Y	0.1509	0.3616
D	0.0061	0.5045
Y × D	0.0001	<0.0001

Means followed by the same letters within the same category are statistically similar according to Duncan’s multiple range test at p < 0.05.

and

Table 2. Comparison of seed cotton and lint yield of various plant densities in 2016 and 2017 in the cotton growing season.

Treatment	Seed Cotton Yield (kg ha−1)	Lint Yield (kg ha−1)	Lint Percentage (%)
D1	3258e	1389e	42.8a
D2	3598d	1490d	41.6ab
D3	3989c	1574c	39.5bc
D4	4304b	1669b	38.8c
D5	4662a	1763a	37.9c
D6	4259b	1609bc	37.8c
ANOVA
Y	0.0008	<0.0001	0.8173
D	0.0005	0.0002	0.0700
Y × D	0.1827	0.7746	0.2969

Means followed by the same letters within the same category are statistically similar according to Duncan’s multiple range test at p < 0.05.

). During both years, D5 plant density (PD) produced the highest seed cotton and lint yield as compared to other plant densities. Highest seed cotton yield of 4662 kg ha⁻¹ and highest lint yield 1763 kg ha⁻¹ was produced by D5, which was followed by D4, D6, D3, D2, and D1. The highest lint percentage (43.5%) was recorded at D1, followed by D2, D3, D4, D5, and D6. The boll density per unit ground area generally increased with increasing plant density but the boll density of individual plants decreased with increasing plant density. More number of bolls m⁻² (105.4) was produced by D6 in 2016, while in 2017 more bolls m⁻² (75.7) was produced by D5. During both years, bigger bolls were produced by D1 as compared to other treatments.

3.2. Biomass Accumulation

Cotton plant biomass accumulation (CPB) was significantly affected by plant density and followed a normal logistic model by DAE (Figure 2). CPB increased as plant density increased and differences were found between the different densities. The D6 plant density had more CPB accumulation as compared to D1, D2, D3, D4, and D5 during both years. Vegetative organ biomass (VOB) during 2016 and 2017 was positively affected by plant density (

Table 3. Analysis of variance for the effect of year (Y) and plant density (PD) on biomass accumulation.

Source	42 DAE		57 DAE		72 DAE		87 DAE		102 DAE		117 DAE		132 DAE
	F	p-Value	F	p-value	F	p-Value	F	p-Value	F	p-Value	F	p-Value	F	p-Value
VOB
Y	9.808	0.026	52.945	0.001	1.485	0.277	0.126	0.737	8.859	0.031	13.983	0.013	21.422	0.006
D	18.048	0.003	10.969	0.010	29.071	0.001	80.460	0.000	24.908	0.002	16.879	0.004	14.033	0.006
Y × D	4.421	0.007	14.299	<0.0001	5.652	0.002	1.710	0.178	14.539	<0.0001	27.077	<0.0001	21.498	<0.0001
ROB
Y	39.283	0.002	14.100	0.013	15.211	0.011	31.227	0.003	19.709	0.007	1.360	0.296	27.305	0.003
D	5.136	0.048	3.505	0.097	54.730	0.000	20.054	0.003	140.225	<0.0001	61.440	0.000	25.162	0.001
Y × D	2.878	0.041	5.808	0.002	1.842	0.150	4.177	0.009	0.588	0.709	0.880	0.513	6.975	0.001
CPB
Y	10.214	0.024	50.051	0.001	6.404	0.052	108.192	0.000	5.092	0.074	5.802	0.061	23.851	0.005
D	17.733	0.003	10.436	0.011	67.943	0.000	602.531	<0.0001	103.482	<0.0001	26.573	0.001	17.484	0.003
Y × D	4.492	0.007	13.345	<0.0001	3.051	0.033	0.501	0.772	3.956	0.012	11.936	<0.0001	33.702	<0.0001

and

Table 4. Regression of cotton plant biomass accumulation at growing seasons 2016 and 2017.

Items	Treatment	Regression Equation	R2
Cotton plant biomass (2016)	D1	Y = 697.81055/(1 + 348.84488e−0.06161t)	0.9925 ***
	D2	Y = 1011.12157/(1 + 310.52754e−0.06365t)	0.9909 ***
	D3	Y = 1302.14044/(1 + 174.95064e−0.05567t)	0.9878 ***
	D4	Y = 1435.81604/(1 + 134.01637e−0.05418t)	0.9783 ***
	D5	Y = 1575.89801/(1 + 138.95428e−0.05546t)	0.9844 ***
	D6	Y = 1594.84306/(1 + 111.53485e−0.05384t)	0.9741 ***
Vegetative organ biomass	D1	Y = 295.33909/(1 + 298.23496e−0.06989t)	0.9879 ***
	D2	Y = 484.50669/(1 + 743.27663e−0.08699t)	0.9796 ***
	D3	Y = 651.12066/(1 + 143.02059e−0.06253t)	0.9782 ***
	D4	Y = 701.83792/(1 + 131.19935e−0.06514t)	0.9569 ***
	D5	Y = 773.54645/(1 + 108.47096e−0.06237t)	0.9738 ***
	D6	Y = 807.8694/(1 + 122.95547e−0.06699t)	0.9549 ***
Reproductive organ biomass	D1	Y = 403.53388/(1 + 1096.07473e−0.06673t)	0.9918 ***
	D2	Y = 494.27754/(1 + 2706.97502e−0.07808t)	0.9956 ***
	D3	Y = 591.54947/(1 + 1594.68086e−0.07241t)	0.9943 ***
	D4	Y = 668.08731/(1 + 1195.61028e−0.07009t)	0.9946 ***
	D5	Y = 741.0549/(1 + 1078.67833e−0.07091t)	0.9915 ***
	D6	Y = 734.4757/(1 + 978.71309e−0.06826t)	0.9929 ***
Cotton plant biomass (2017)	D1	Y = 648.96068/(1 + 258.79492e−0.06442t)	0.9936 ***
	D2	Y = 838.50226/(1 + 150.01915e−0.06353t)	0.9886 ***
	D3	Y = 951.32332/(1 + 103.58331e−0.06163t)	0.9857 ***
	D4	Y = 1070.21853/(1 + 83.38347e−0.05897t)	0.9918 ***
	D5	Y = 1144.95086/(1 + 95.96115e−0.06213t)	0.9936 ***
	D6	Y = 1165.41818/(1 + 74.52055e−0.05987t)	0.9938 ***
Vegetative organ biomass	D1	Y = 268.83729/(1 + 327.895e−0.08594t)	0.9911 ***
	D2	Y = 404.49245/(1 + 299.20271e−0.0908t)	0.9883 ***
	D3	Y = 460.00848/(1 + 292.09854e−0.09749t)	0.9855 ***
	D4	Y = 505.9789/(1 + 101.20589e−0.08053t)	0.9898 ***
	D5	Y = 534.94333/(1 + 88.38951e−0.07954t)	0.9911 ***
	D6	Y = 586.38098/(1 + 97.73996e−0.08321t)	0.9890 ***
Reproductive organ biomass	D1	Y = 354.07086/(1 + 8856.53541e−0.09421t)	0.9938 ***
	D2	Y = 415.62152/(1 + 2397.24183e−0.08411t)	0.9860 ***
	D3	Y = 473.372/(1 + 1661.97624e−0.08182t)	0.9825 ***
	D4	Y = 539.06687/(1 + 1185.99181e−0.07956t)	0.9877 ***
	D5	Y = 582.96854/(1 + 1712.18326e−0.08643t)	0.9869 ***
	D6	Y = 548.10109/(1 + 1929.68342e−0.08696t)	0.9867 ***

***, significant at the 0.001 probability level.

). The VOB increased linearly with the increase in plant density. The highest PD, D6, produced more VOB as compared to other plant densities while individual plant VOB decreased as density increased due to resource competition among plants. Reproductive growth of cotton started from the appearance of the first square. Less biomass accumulated to reproductive organs of cotton which increases linearly with further growth. Treatment D5 produced more ROB in 2016 and 2017, followed by D6, D4, D3, D2, and D1.

3.3. Simulation of Biomass Accumulation

Simulation of biomass accumulation based on Equation (1) followed the logistic function and all the biomass accumulation were significant. Calculation from Equations (2)–(7) based on Table 2 illustrates the day of starting and termination of cotton biomass fast accumulation period (FAP) during 2016 and 2017. The averaged highest speed for CPB in all plant densities were 68 and 114 DAE in 2016, and 56 and 98 in 2017, with the highest average (V_A = 16 and 14 gm m⁻² d⁻¹) and maximum rate (V_M = 18 and 15 gm m⁻² d⁻¹) (

Table 5. Eigen values of cotton biomass accumulation at growing season 2016.

Items	Treatment	Fast Accumulation Period				Fastest Accumulation Point
		t1 (DAE)	t2 (DAE)	∆t (d)	VA (gm m−2 d−1)	VM (gm m−2 d−1)	at DAE
Cotton plant biomass	D1	73.7	116.4	42.8	9.4	10.7	95.0
	D2	69.5	110.8	41.4	14.1	16.1	90.2
	D3	69.1	116.4	47.3	15.9	18.1	92.8
	D4	66.1	114.7	48.6	17.1	19.4	90.4
	D5	65.2	112.7	47.5	19.2	21.8	89.0
	D6	63.1	112.0	48.9	18.8	21.5	87.6
	Average	67.8	113.9	46.1	15.7	18.0	90.8
Vegetative organ biomass	D1	62.7	100.4	37.7	4.5	5.2	81.5
	D2	60.9	91.1	30.3	9.2	10.5	76.0
	D3	58.3	100.4	42.1	8.9	10.2	79.4
	D4	54.6	95.1	40.4	10.0	11.4	74.9
	D5	54.0	96.3	42.2	10.6	12.1	75.1
	D6	52.2	91.5	39.3	11.9	13.5	71.8
	Average	57.1	95.8	38.7	9.2	10.5	76.5
Reproductive organ biomass	D1	85.2	124.6	39.5	5.9	6.7	104.9
	D2	84.4	118.1	33.7	8.5	9.6	101.2
	D3	83.7	120.0	36.4	9.4	10.7	101.8
	D4	82.3	119.9	37.6	10.3	11.7	101.1
	D5	79.9	117.1	37.1	11.5	13.1	98.5
	D6	81.6	120.2	38.6	11.0	12.5	100.9
	Average	82.8	120.0	37.1	9.4	10.7	101.4

t₁ is the starting and t₂ is the termination point of the fast accumulation period (FAP). ∆t is the total duration of FAP. V_A is the average and V_M is the maximum rate of biomass accumulation during FAP. DAE represents days after emergence.

and

Table 6. Eigen values of cotton biomass accumulation at growing season 2017.

Items	Treatment	Fast Accumulation Period				Fastest Accumulation Point
		t1 (DAE)	t2 (DAE)	∆t (d)	VA (gm m−2 d−1)	VM (gm m−2 d−1)	at DAE
Cotton plant biomass	D1	65.8	106.7	40.9	9.2	10.5	86.2
	D2	58.1	99.6	41.5	11.7	13.3	78.9
	D3	53.9	96.7	42.7	12.9	14.7	75.3
	D4	52.7	97.3	44.7	13.8	15.8	75.0
	D5	52.3	94.7	42.4	15.6	17.8	73.5
	D6	50.0	94.0	44.0	15.3	17.4	72.0
	Average	55.5	98.2	42.7	13.1	14.9	76.8
Vegetative organ biomass	D1	52.1	82.7	30.6	5.1	5.8	67.4
	D2	48.3	77.3	29.0	8.1	9.2	62.8
	D3	44.7	71.7	27.0	9.8	11.2	58.2
	D4	41.0	73.7	32.7	8.9	10.2	57.3
	D5	39.8	72.9	33.1	9.3	10.6	56.3
	D6	39.2	70.9	31.7	10.7	12.2	55.1
	Average	44.2	74.9	30.7	8.6	9.9	59.5
Reproductive organ biomass	D1	82.5	110.5	28.0	7.3	8.3	96.5
	D2	76.9	108.2	31.3	7.7	8.7	92.5
	D3	74.5	106.7	32.2	8.5	9.7	90.6
	D4	72.4	105.5	33.1	9.4	10.7	89.0
	D5	70.9	101.4	30.5	11.0	12.6	86.1
	D6	71.9	102.1	30.3	10.4	11.9	87.0
	Average	74.8	105.7	30.9	9.1	10.3	90.3

t₁ is the starting and t₂ is the termination point of fast accumulation period (FAP). ∆t is the total duration of FAP. V_A is the average and V_M is the maximum rate of biomass accumulation during FAP. DAE represents days after emergence.

).

Cotton plant biomass accumulation was found significant among plant densities. In 2016, a fast accumulation period in D5 started at 65 DAE and terminated at 113 DAE, which lasts for 48 DAE with the ighest average (19.2 V_A gm m⁻² d⁻¹) and maximum rate (22 V_M gm m⁻² d⁻¹) at 89 DAE. The lengthiest fast accumulation period for CPB was noted in D6, which lasts for 49 DAE with the average rate of 18.8 V_A gm m⁻² d⁻¹ (Table 5).

The fast accumulation period of CPB in 2017 for D6 started earlier at 50 DAE and terminated at 94 DAE, while D1 FAP terminated last at 107 DAE. The highest average (15.6 V_A gm m⁻² d⁻¹) and maximum rate (17.8 V_M gm m⁻² d⁻¹) were noted in D5, followed by D4, D6, D3, D2, and D1 (Table 6).

Vegetative organ biomass responded positively to plant density. The earliest and highest FAP of VOB in both years was observed at D6 with the average rate (12 and 10.7 V_A gm m⁻² d⁻¹), which lasts for 39 and 32 DAE, and maximum rate (13.5 and 12.2 V_M gm m⁻² d⁻¹), which lasts for 72 and 55 DAE in 2016 and 2017, respectively. Both average and maximum VOB accumulation rates of D6 were 62%, 23%, 25%, 16%, and 11% higher than D1, D2, D3, D4, D5 in 2016 and 52%, 24%, 8%, 17%, and 13% higher than D1–D5 in 2017 (Table 5 and Table 6).

The highest average rate (11.5 V_A gm m⁻² d⁻¹) of reproductive structures biomass was observed in D5, which started at 80 DAE and terminated at 117 DAE and lasted for 37 DAE, with a maximum rate (13 V_M gm m⁻² d⁻¹) at 99 DAE in 2016 (Table 4). Both average and maximum ROB accumulation rates of D5 were observed to be higher as compared to D1, D2, D3, D4, and D5. The earliest FAP in D5 began at 80 DAE, while the last terminated FAP was observed in D1, which ended at 125 DAE. In 2017, initial FAP of ROB began in D5 which lasted for 31 DAE and terminated at 101 DAE, with the highest average rate (11 V_A gm m⁻² d⁻¹) and maximum rate (12.6 V_M gm m⁻² d⁻¹) at 86 DAE, followed by D6, D4, D3, D2, and D1 (Table 5 and Table 6).

3.4. Dry Matter Partitioning (DWRO/PB)

Dry matter partitioning, as indicated by the ratio of the dry weight of reproductive organs to plant biomass (DWRO/PB), increase slowly as the plant changes from one growth stage to another and peak stage of dry matter partitioning was observed at 120 DAE during 2016 and 2017 (Figure 3). During different growth stages, significant differences were observed between treatments (

Table 7. Analysis of variance for the effect of year (Y) and plant density (PD) on dry matter partitioning.

Source	42 DAE		57 DAE		72 DAE		87 DAE		102 DAE		117 DAE		132 DAE
	F	p-Value	F	p-Value	F	p-Value	F	p-Value	F	p-Value	F	p-Value	F	p-Value
Y	9.371	0.028	3.769	0.110	0.178	0.691	8.596	0.033	15.602	0.011	205.764	0.0000	13.113	0.015
D	64.200	0.0000	1.009	0.496	1.112	0.455	1.415	0.356	3.231	0.112	23.395	0.002	18.230	0.003
Y × D	0.591	0.707	26.356	0.000	5.720	0.002	6.189	0.001	1.904	0.139	0.249	0.935	1.498	0.235

). The DWRO/PB of D1 was observed to be higher as compared to other treatments.

3.5. Boll Load (DWRO/LA)

Boll load, as indicated by the ratio of the dry weight of reproductive organs by leaf area (DWRO/LA), was found to be significantly higher in D6 as compared to other treatments and significant differences were observed between different treatments (

Table 8. Analysis of variance for the effect of year (Y) and plant density (PD) on boll load.

Source	42 DAE		57 DAE		72 DAE		87 DAE		102 DAE		117 DAE		132 DAE
	F	p-Value	F	p-Value	F	p-Value	F	p-Value	F	p-Value	F	p-Value	F	p-Value
Y	59.210	0.001	1.984	0.218	3.392	0.125	9.491	0.027	9.909	0.025	4.792	0.080	4.161	0.097
D	1.902	0.249	2.038	0.227	8.939	0.016	9.628	0.013	8.725	0.016	20.358	0.002	32.806	0.001
Y × D	5.208	0.003	1.565	0.215	6.946	0.001	7.186	0.001	5.631	0.002	4.300	0.008	1.138	0.373

). The DWRO/LA increased gradually with an increase in plant density and changing from one growth stage to another (Figure 4). At 132 DAE, DWRO/LA of D6 was 14%–82% and 4%–76% higher than treatment D1–D6 during 2016 and 2917, respectively.

3.6. Leaf Area Index

The leaf area index (LAI) at different days after emergence is shown in Figure 5. The LAI of D6 was higher during both years as compared to other treatments and increased linearly as plant density increased (

Table 9. Analysis of variance for the effect of year (Y) and plant density (PD) on leaf area index.

Source	42 DAE		57 DAE		72 DAE		87 DAE		102 DAE		117 DAE		132 DAE
	F	p-Value	F	p-Value	F	p-Value	F	p-Value	F	p-Value	F	p-Value	F	p-Value
Y	25.001	0.004	15.040	0.012	0.413	0.549	0.030	0.869	2.068	0.210	7.016	0.045	12.421	0.017
D	5.591	0.041	16.289	0.004	20.785	0.002	37.873	0.001	76.869	0.000	20.444	0.002	28.831	0.001
Y × D	7.678	0.000	1.573	0.213	1.237	0.329	0.896	0.503	1.032	0.426	3.175	0.029	0.635	0.676

). The LAI increased with the growth of the cotton plants and reached a peak at 102 DAE and then decreased linearly. LAI of high plant density reached 4.3 in 2016 and 4 in 2017, while in the case of lower plant density, it reached up to 1.3 in 2016 and 1.5 in 2017. In the last growth stages, no significant differences were observed in plant densities D4–D6.

3.7. Fiber Quality

Fiber quality parameters were significantly influenced by plant density in both years (

Table 10. Comparison of fiber quality parameters of various densities of growing seasons 2016 and 2017.

Treatment	Length	Strength	Elongation	Uniformity Index	Microniare
	(mm)	(cN/tex)	(%)	(%)
D1	28.86c	31.26a	6.40a	84.77b	5.33a
D2	29.40bc	30.97ab	6.28a	84.83b	4.93b
D3	29.52bc	30.10b	6.25a	85.17ab	4.83bc
D4	29.82b	29.10c	6.24a	85.47ab	4.80bc
D5	30.63a	28.80c	6.20a	86.07a	4.70c
D6	30.60a	28.33c	6.22a	85.53ab	4.63c
ANOVA
Y	0.1170	0.0798	0.0000	<0.0001	0.0001
D	0.8030	0.0001	0.5164	0.0109	0.0007
Y × D	0.0614	0.9776	0.4447	0.9268	0.8639

Means followed by the same letters within the same category are statistically similar according to Duncan’s multiple range test at p < 0.05.

). An increase in plant density led to longer fiber length while decreasing strength and micronaire value. Low planting density had low length and greater strength and micronaire value as compared to high treatments. The fiber length of high plant density D5 and D6 were statistically similar while the length of D5 was 5% and 7% longer than the lowest planting density in both years, respectively. Fiber strength of lowest plant density was recorded to be 9% and 10% higher, while micronaire was observed to be 15% and 9% higher as compared to the highest planting density in 2016 and 2017, respectively. Planting density had no significant effect on fiber elongation and uniformity index during both growing seasons.

3.8. Economic Analysis

The net returns were affected by different plant density. Net returns were determined on the basis of production cost and returns from the cotton crop. The highest net returns were obtained from D5 (1750 USD ha⁻¹ and 1393 USD ha⁻¹) during 2016 and 2017, respectively, while the lowest was obtained from D1 (

Table 11. Cotton yield, production cost, and net returns of 2016 and 2017.

Treatment	Lint Yield		Returns from Cotton		Production Cost		Net Returns
	Kg ha−1		USD ha−1		USD ha−1		USD ha−1
	2016	2017	2016	2017	2016	2017	2016	2017
D1	1459	1320	3040	2705	2429	2389	611	316
D2	1588	1391	3310	2850	2325	2257	985	593
D3	1685	1463	3510	2998	2297	2259	1213	739
D4	1762	1576	3671	3230	2269	2232	1402	998
D5	1843	1683	3840	3449	2090	2056	1750	1393
D6	1726	1493	3589	3061	2137	2102	1452	959

Production cost includes fertilizer, seed, and labor cost. Labor cost includes labor for planting, management, and harvesting. One labor unit per day cost was 6.02 USD in 2016 and 5.92 USD in 2017. Wholesale lint price of 1 kg was 2.08 USD in 2016 and 2.05 USD in 2017. Values were converted from Chinese Yuan to USD according to the official rate (USD 1 = 6.65 yuan in 2016 and 6.76 yuan in 2017).

). Seed and labor cost mostly affected net returns. More labor was required for low plant densities due to more vegetative branches as compared to high plant densities.

4. Discussion

The main purpose of this study is to explore and compare different plant densities in response to cotton yield, leaf area index, dry matter partitioning, and fiber quality at different growth stages. Higher plant density is the key management practice for obtaining greater numbers of bolls per unit area, but in most cases, the yield enhances up to an optimum density, after which further increase in plant population decreases yield. Different regions of China have different optimum densities and lint production, which depends on climatic conditions along with other management practices. The Xinjiang autonomous region has the recommended PD of 21.0 × 10⁴ to 24.0 × 10⁴ plants ha⁻¹ [27]; followed by Yellow River Valley with a PPD of 3.0 × 10⁴, 4.5 × 10⁴, and 6.0 × 10⁴ plants ha⁻¹ for hybrid Bt cotton, indigenous Bt cotton, and Bt cotton, respectively [28,29], and for late sowing, PD is 7.5 × 10⁴ ha⁻¹ [30]; while in the Yangtze River Valley where hybrid seeds are commonly used, it has the PD of 3.0 × 10⁴ plants ha⁻¹ [31]. Our results are consistent with previous studies that have shown that cotton yield increases up to a certain limit with increasing PD, while too low and too high plant density cause a reduction in yield [32]. In this study, yield and yield components were significantly affected by plant density, excluding boll weight and lint percentage. High yield and yield components were noted in plant density D5. Yield and number of bolls produced by a single plant of the treatment D5 was lower as compared to other treatments but was more based on per unit area. These results are consistent with Mao et al. [33], who reported that high plant population increase bolls m⁻² while the weight of individual bolls decreases.

More biomass production is the foundation of high yield [34,35,36]. In this study, biomass accumulation was higher in 2016 as compared to 2017, which might be due to differences in environmental conditions. Total plant biomass and vegetative organ biomass accumulation were high in higher plant density while higher reproductive organ biomass was accumulated in D5 as compared to other treatments. In early growth stages of the cotton plant, plant density did not affect reproductive structure biomass accumulation, while after 87 DAE ROB, accumulation was influenced significantly. High biomass accumulation in high plant density was due to a greater number of plants per unit ground area with more vegetative growth. Our results are in line with other researchers who also reported that high plant density resulted in high biomass production [37,38]. Both high and low plant density lead to reductions in reproductive organ biomass. The less ROB production in high population might be due to less light penetration to the lower parts of plants, followed by a reduction in temperature and increased relative humidity in the cotton canopy, which enhanced fruit shedding as compared to other plant densities [39].

The ratio of dry weight of reproductive organs to plant biomass (DWRO/PB) also affected the yield of cotton [40]. In this study, the highest ratio was obtained in the lowest plant density D1 as compared to other treatments which showed less differences. Similar results were previously obtained by Dai et al. [37], who also reported high DWRO/PB in the lowest density. Boll load is also an important indicator of lint yield. In this study, a high and significant ratio of dry weight of reproductive organs per leaf area (DWRO/LA) was observed in high plant densities (D4–D6), mostly in the late growth stages. Our results are supported by Dong et al. [30], where high boll load led to an increase in leaf senescence and a decrease in cotton yield and quality. The high DWRO/LA in late growth stages is due to high competition for nutrients and assimilates between vegetative and reproductive growth after the bloom stage [40].

Leaf area index is an important factor that affects biomass production of cotton [41]. LAI is also one of the physiological parameters which determine crop yield and predict crop production up to some extent. For obtaining high yield, it is necessary to maintain optimum LAI for more light penetration and high light use efficiency, mostly at late growth stages: that is, the flowering and boll setting stages [22].

Cotton fiber is the extension of seed epidermal cells. Fiber quality indicators are affected by plant density and environmental factors [10]. In the present study, cotton fiber indicators were significantly affected by plant density. Low plant density had high strength and micronaire value as compared to high plant density, while the length of low plant density is shorter as compared to high and moderate plant density. Our results are in agreement with previous research that have reported high strength and micronaire and short fiber length at low planting density [5,33]. The lower fiber quality at high planting density may be due to less photosynthesis, which reduces carbohydrate supply for fiber formation. For obtaining good quality fiber, cultivar selection is of great importance, while managing plant populations to maintain genetic potential is the secondary part [42,43].

Economic benefit plays an important role in the success of agriculture business. In the Yellow River Valley, due to fast urbanization, high labor costs and a shortage of labor have become a challenge to traditional intensive cotton production [2]. Labor cost specifically affects the profitability of the cotton crop. High density has less vegetative branches as compared to high plant density, which needs less labor for vegetative branch removal and other field management.

5. Conclusions

In the present study, planting density positively affected yield, fiber quality, and dry matter accumulation and partitioning of cotton crop under equidistant row arrangement. Optimum or moderate plant density (8.7 plants m⁻²) resulted in high reproductive organ biomass accumulation at later growth stages as compared to other treatments. More reproductive organ biomass accumulation in this density increased the yield of cotton. Good quality fiber was obtained at low and moderate plant densities as compared to higher ones. In conclusion, 8.7 plants m⁻² is regarded to be an optimum plant density in term of high yield, uniform fiber quality, and dry matter accumulation. The finding of this research offers an alternative to cotton growers who use conventionally wider rows and lower plant population ha⁻¹.